INTERIM  REPORT 

O.R.A.  155-421 

Contract  #  702930013 

A  HYDROLOGIC  STUDY  OF  THE  SULPHUR  ARTESIAN  GROUNDWATER  SYSTEM 
AND  ASSOCIATED  WATERS  AT  CHICKASAW  NATIONAL  RECREATIONAL  AREA, 
SULPHUR,  OKLAHOMA 


Prepared  By: 

Jim  F.  Harp,  Ph.D.,  Professor 
Joakim  G.  Laguros,  Ph.D.,  Professor 
Stephen  S.  McLin,  Ph.D.,  Assoc.  Professor 
Leonard  B.  West,  Ph.D.,  Assoc.  Professor 
Charles  John  Barthel,  Jr. 


Civil  Engineering  and  Environmental  Science  Dept. 
202  West  Boyd,  Room  334,  Carson  Engineering  Center 
The  University  of  Oklahoma,  Norman,  Oklahoma  7  3019 


Submitted  To: 

The  United  States  Department  of  the  Interior 
National  Park  Service 


Submitted  By: 

The  Bureau  of  Water  Resources  Research, 
Oklahoma  Research  Administration 
The  University  of  Oklahoma 
Norman,  Oklahoma 


June  19  85 


TABLE  OF  CONTENTS 

Page 

LIST  OF  TABLES vii 

LIST  OF  FIGURES x 

LIST  OF  PLATES xiv 

ABSTRACT XV 

Chapter 

I  INTRODUCTION  1 

1.1  The  Meaning  of  Hydrogeology 1 

1 . 2  Man  and  Springs 1 

1.3  Scope  of  the  Project 3 

I I  BACKGROUND 5 

2 . 1  Geography  and  Hydrology 5 

2.2  Geology  of  the  Area 13 

2.3  Groundwater  Flow  Through 

Geologic  Materials 27 

2.4  Spring  Flow  Assessment  and 

Related  Methodology 4  8 

2.5  Previous  Hydrogeologic  Findings  of 

Artesian  Spring  Flow, 

Sulphur,  Oklahoma,  Area 7  3 

III  METHODOLOGY 98 

3.1   Geologic  Control  of  Groundwater  Flow....  98 


TABLE  OF  CONTENTS  (continued) 


Page 


3.2  Determination  of  the  Physical  Properties 

of  the  Simpson-Arbuckle  Aquifer 101 

3.3  Wells  and  Springs 110 

3.4  Determination  of  Recharge  Areas 113 

IV  FINDINGS 116 

4  . 1   Geology 116 

4  .  2   Petrophysical  Results 130 

4.3  Assessment  of  Aquifer  Properties 137 

4.4  Wells  and  Springs 162 

4.5  Recharge  Areas  Responsible  for  Spring 

Flow  in  the  Travertine  District, 

Chickasaw  National  Recreational  Area.  184 

4.6  Concise  Nomenclature  and  Origin  of  the 

Fresh  Water  Springs,  Buffalo  and 

Ante  lope 191 

4.7  Hydrologic  Budget  of  the  Sulphur  Area...  192 

V  DISCUSSIONS  OF  FINDINGS 193 

5.1   Geologic  Control  of  Groundwater 193 

5 .  2   Aquifer  Properties 195 

5.3  Data  Analysis  of  Wells  and  Springs 198 

5.4  Analysis  of  Spring  Flow  Occurrence  and 

the  Recharge  Areas 201 

VI  SUMMARY  AND  CONCLUSIONS 2  05 

NOMENCLATURE 209 

BIBLIOGRAPHY 213 

APPENDIX  I 218 


TABLE  OF  CONTENTS  (continued) 

Page 

APPENDIX  II 223 

APPENDIX  III 229 

APPENDIX  IV 2  3  2 


LIST  OF  TABLES 

Table  Page 

1.  Surface  Elevations  of  Springs 10 

2 .  Rainfall  Records 14 

3.  Typical  Porosities  for  Porous  Materials 36 

4.  Calculated  Discharge  Parameters  for  Selected 

Karst  Springs 53 

5.  Meinzer's  Classification  of  Spring 

Discharge 55 

6.  Water  Quality  Analysis  of  Selected 

Wells  and  Springs  in  Project  Area 83 

7.  Relative  Effects  of  Production  Wells 

in  the  Vicinity  of  Antelope  and  Buffalo...     87 

8.  Effect  of  Hydrogeologic  Setting  on 

Carbonate  Aquifers 99 

9.  Fracture  Measurements  of  the  Pontotoc 

Group,  Antelope  Springs 119 

10.  Fracture  Measurements  of  the  Pontotoc 

Group 121 

11.  Fracture  Measurements  of  an  Indurated 

Sandstone  of  the  Simpson  Group,  Glass 

Sand  Quarry 122 

12.  Fracture  Measurements  of  the  Arbuckle 

Group,  Mill  Creek  and  State  Highway  7 124 

13.  Fracture  Measurements  of  the  Arbuckle  Group, 

Intersection  of  State  Highways  7  and  12...    125 

14.  Fracture  Measurements  of  the  Arbuckle 

Group,  Collins1  Ranch 126 

vii 


Vlll 


LIST  OF  TABLES  (continued) 
Table  Page 

15.  Laboratory  Data  and  Effective  Porosity 

Measurements  of  Selected  Rock  Specimens...    131 

16.  Slug  Test  Data  of  Crawford's  Well 138 

2 

17.  Input  Parameters  (T  =  .1  ft  /hour)  and 

Resultant  Output 145 

2 

18.  Input  Parameters  (T  =  1.0  ft  /hour)  and 

Resultant  Output 149 

19.  Water  Level  Fluctuations  and  Barometric 

Pressure  Changes 158 

20.  Determination  of  Average  Aquifer  Trans- 

missivity  as  a  Function  of  Spring 

Discharge 161 

21.  Ground  Water  Levels,  East  Observation  Well...  163 

22.  Ground  Water  Levels,  West  Observation  Well...  164 

23.  Stage  Data,  Buffalo  Springs 165 

24.  Stage  Data,  Antelope  Springs 166 

25.  Spring  Flow  in  Cubic  Feet  Per  Second, 

Buffalo  Springs 167 

26.  Spring  Flow  in  Cubic  Feet  Per  Second, 

Antelope  Springs 168 

27.  Piezometric  Levels  of  Selected  Wells  in 

the  Arbuckle  Group 172 

28.  Well  Pumpage  Data,  Area  Well  Fields 176 

29.  Mean  Annual  Flow  of  Antelope  and  Buffalo 

Springs,  and  Annual  Precipitation 179 

30.  Mean  Precipitation  and  Mean  Flow  of 

the  Springs  on  a  Monthly  Basis 181 

31.  Input  Variables  and  Resultant  Values  of 

Spring  Variance  for  Antelope  and 

Buffalo  Springs 182 


IX 

LIST  OF  TABLES  (continued) 
Table  Page 

32.  Data  Input  and  Resultant  Values  of 

Exhaustion  Coefficient  and  Response 

Time  of  Antelope  Spring 183 

33.  Results  of  the  Potentiometric  Survey 185 

A.  1   ISTEP  Function  Values 227 


LIST  OF  FIGURES 


Figure  Page 

1.  Location  of  project  area,  Murray  County, 

Oklahoma 6 

2.  Location  of  Sulphur  in  relation  to 

the  Arbuckle  Mountains 7 

3 .  Map  showing  the  eastern  and  western  regions 

of  the  Arbuckle  Mountains 8 

4.  Map  depicting  location  of  springs  in  the 

Travertine  District  of  Chickasaw  National 
Recreational  Area 9 

5.  Regional  stratigraphic  section  across  the 

Southern  Oklahoma  Geosyncline 16 

6.  Stratigraphic  columns  (Pre-Pennsylvanian)  of 

the  Arbuckle  Mountains 18 

7.  Geologic  map  of  the  Arbuckle  Mountains 20 

8 .  Surface  contour  map  of  the  basement  rock  of 

the  Arbuckle  Mountain  area 2  5 

9.  Hydrologic  cycle 28 

10.  A  sketch  of  an  apparatus  to  demonstrate 

Darcy '  s  Law 31 

11.  Approximate  range  of  hydraulic  conductivity 

in  soil  and  rock 31 

12.  A.  Water-table  aquifer 

B  .  Artesian  aquifer 40 

13.  Elemental  volume  of  flow  through  a  confined 

aquifer 40 


x. 


XI 


LIST  OF  FIGURES  (continued) 


Figure  Page 

14.  Equipotential  lines,  h.  through  hg ,  in  a 

two-dimensional  flow  field  and  the 

gradient  of  h 49 

15.  A  schematic  representation  of  the  scale 

on  karst  permeability  or  hydraulic 

conductivity 4  9 

16.  Schematic  hydrograph  for  high  discharge  from 

a  carbonate  spring 51 

17.  Two-dimensional  fracture  pattern  showing 

the  relation  between  true  flow  and 

direction  inferred  by  drawing  orthogonal 

lines  to  the  potentiometric  contours 57 

18.  Histogram  depicting  frequencies  of  fractures 

for  various  lithologies 57 

19.  Four  stages  of  karstic  landform  development..     59 

20.  Permeability,  open  fissure  density,  and 

aggregate  fissure  density  of  limestone 

rocks  versus  the  geologic  time  scale 60 

21.  Parts  A  and  B  depict  flow  systems  as  a  func- 

tion of  basin  depth,  topography,  and 

potential  distributions 62 

22.  A.  Hydrograph  of  observation  well  depicting 

changes  in  piezometric  level  in  relation 
to  spring  discharge 
B.  Hydrograph  of  observation  well  showing 

long-term  deprivation  of  irrigation 64 

23.  Drawdown  and  build-up  curves  for  observation 

wells  in  a  limestone  karst  aquifer 69 

24.  A  sketch  depicting  water  table  depth 

measurement  utilizing  a  steel  surveyor's 

tape  that  can  be  read  with  an  accuracy 

of  .005  feet 69 

25.  The  Stiff  graphical  method 71 


Xll 


LIST  OF  FIGURES  (continued) 


Figure  Page 

26.  Water  quality  analysis  of  the  Sulphur, 

Oklahoma,  area,  by  modified  Stiff 

diagrams  in  milliequivalents  per  liter....      79 

27.  Location  of  selected  wells  and  springs; 

Sulphur,  Oklahoma,  area 8  2 

28.  Well  layout 84 

29.  Structural  contour  map  of  the  surface 

of  the  Arbuckle  Group 128 

30.  Geologic  and  hydrogeologic  map  of  the 

Sulphur,  Oklahoma,  area 129 

31.  Magnification  is  39X;  mosaic  of  aphanitic 

dolomite  anhedra  of  the  West  Spring 

Formation,  Arbuckle  Group 133 

32.  Magnification  is  39X;  poorly-sorted, 

calcite-cemented  calcarenite  of  the 

Oil  Creek  Formation,  Simpson  Group 134 

33.  Magnification  is  39X;  calcite-cemented, 

fairly  well  sorted,  subangular  to  sub- 
rounded  arenite  of  the  Basal  Oil  Creek 
Formation,  Simpson  Group 135 

34.  Magnification  is  39X;  conglomerate  of  the 

Pontotoc  Group  consisting  of  some  fine 

to  medium  grained  quartz  particles,  and 

fine  grained  to  pebble-size  lime  clasts, 

embedded  in  a  calcite  and  argillaceous 

cement 136 

35.  Results  of  a  slug  test  performed  on  a 

selected  well  in  the  Arbuckle  Group 

outcrop 139 

36.  Prickett-Lonnquist  model  (1971)  that 

was  used  to  simulate  well  drawdown  at 

a  sink  node 141 

37.  Multi-purpose  water  well  report  on 

Municipal  Well  No.  9,  city  of  Sulphur 144 


Xlll 


LIST  OF  FIGURES  (continued) 


Figure  Page 

38.  Strip  chart  from  a  Stevens'  Water-Level 

Recorder  positioned  at  the  East  Observa- 
tion Well,  depicting  continuous  moni- 
toring of  water  level  fluctuations, 
h  (feet) 154 

39.  Strip  chart  from  a  Stevens'  Water-Level 

Recorder  positioned  at  the  West  Observa- 
tion Well,  depicting  continuous  moni- 
toring of  water  level  fluctuations, 
h  (feet) 155 

40.  Strip  charts  from  a  microbarograph  posi- 

tioned at  the  Travertine  Nature  Center, 
depicting  continuous  readings  of  baro- 
metric pressure,  P   (inches  of  Hg) 156 

3. 

41.  Graphical  solution  to  obtain  the  baro- 

metric efficiency,  BE 159 


42 .  Well  and  spring  hydrographs 169 

43.  Long-term  precipitation  histogram  versus 

periods  of  no  flow  of  the  fresh  water 

springs 173 

44.  Inch-year,  time  bar  chart  that  depicts 

above-  or  below-normal  precipitation 

amounts  for  a  specific  period,  in 

relation  to  no-flow  periods  of  Antelope 

and  Buffalo 175 

45.  Potentiometric  map  (Jan.  and  Feb.,  1985) 

of  the  Simpson-Arbuckle  Aquifer,  Sulphur, 
Oklahoma,  area 190 

46.  Nomograph  that  correlates  recharge  area 

to  spring  discharge  by  annual  recharge 

amount 203 

A.l   Prickett-Lonnquist  Model  for  an  isotropic, 

non-leaky,  confined  aquifer 225 

A. 2   Simplified  flowchart  of  the  Prickett- 
Lonnquist  Model 228 


XIV 


LIST  OF  PLATES 


Plate  Page 

Plate  I      Geology  of  the  Sulphur,  Oklahoma, 

Area in  folder 

Plate  II    Potentiometric  Map  (Winter,  1976- 

1977)  of  the  Eastern  Arbuckle 
Mountains,  Sulphur,  Oklahoma, 
Area in  folder 

Plate  III    Topographic  Map  of  the  Sulphur, 

Oklahoma,  Area in  folder 

Plate  IV    Lineament  Map  of  Sulphur,  Oklahoma, 

Area in  folder 


ABSTRACT 

Increasing  cessation  of  fresh  water  artesian  spring 
flow  in  the  Travertine  District,  Chickasaw  National  Rec- 
reation Area,  prompted  this  investigation  to  address  the 
parameters  that  define  subsurface  flow  in  the  Sulphur, 
Oklahoma  area. 

The  Sulphur  Fault  is  responsible  for  the  origin  of 
the  fresh  water  artesian  springs,  Buffalo  and  Antelope, 
and  controls  groundwater  movement  from  the  recharge  area 
downdip  to  the  springs.  Buffalo  and  Antelope  Springs  de- 
rive their  water  from  the  confined  Arbuckle  Aquifer.  The 
intensely  fractured  dolomites  of  the  Arbuckle  group  out- 
crop east  of  Sulphur,  constituting  the  recharge  area  of 

2 
approximately  5  mi   for  the  fresh  water  springs. 

Observation  wells  and  spring  hydrograph  analysis 

help   to   establish   a   relationship   between   discharge, 

piezometric  elevation  of  the  confined  Arbuckle  Aquifer, 

and  seasonal  precipitation  events.   The  long-term  trend 

of  the  hydrographs  reveal  that  the  piezometric  elevation 

of  the  Sulphur  artesian  basin  is  declining  at  a  rate  of 

.65ft/year.   If  this  trend  continues,  permanent  cessation 

of  spring  flows  will  occur  in  about  40  years. 

xv 


XVI 

Statistical  and  spring  hydrograph  analysis  have 
provided  the  basis  for  the  classification  scheme  of 
Antelope  and  Buffalo  Springs.  As  a  result,  Antelope  and 
Buffalo  are  classified  as  third-magnitude,  intermediate- 
response,  diffuse-flow,  fresh  water  artesian  springs. 

Finally,  a  hydrologic  budget  was  developed  for  the 

Sulphur,  Oklahoma  area.   The  recharge  area  for  the  City 

2 

of  Sulphur  is  approximately  15  mi   with  an  average  annual 

recharge  rate  of  about  4.6  inches  per  year;  the  total 

amount  of  water  available  on  a  yearly  basis  is  approxi- 

9 
mately  1.2  x  10   gallons  per  year.   However,  current 

9 
analysis  reveals  that  in  excess  of  1.4  x  10   gallons  per 

year  is  withdrawn  from  the  Sulphur  artesian  basin. 


CHAPTER  I 


INTRODUCTION 


1.1  The  Meaning  of  Hydrogeology 

Hydrogeology  is  a  discipline  that  is  primarily  con- 
cerned with  the  flow  of  water  that  is  influenced  by  geo- 
logic processes  (Narasimhan,  1982) .  Specifically,  hydro- 
geology entails  surface  water  and  groundwater  interac- 
tions governed  by  the  geometry,  material  distributions, 
and  boundary  and  initial  conditions  of  the  porous  media. 

Occasionally,  geohydrology  and  hydrogeology  are  used 
synonymously  in  the  literature.  Interestingly  enough, 
the  term,  hydrogeology,  does  not  appear  in  the  1981  edi- 
tion of  the  Websters  New  Collegiate  Dictionary.  Yet, 
geohydrology  does,  and  the  definition  is  as  follows: 
"Science  that  deals  with  the  character,  source,  and  mode 
of  occurrence  of  underground  water".  However,  for  the 
purpose  of  this  study,  hydrogeology  connotes  an  emphasis 
on  groundwater  flow  governed  by  geological  parameters. 

1.2  Man  and  Springs 

Mankind,  throughout  recorded  history,  has  been  fas- 
cinated and  attracted  to  springs.   Ancient  civilizations 
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settled  in  close  proximity  to  springs  in  arid  regions 
(Davis  and  DeWiest,  1966).  Mineral  springs  and  thermal 
springs  of  reputed  medicinal  and  therapeutic  value  were 
exploited  by  vigorous  commercial  ventures  that  resulted 
in  the  establishment  of  permanent  settlements  in  the 
United  States.  The  city  of  Sulphur  in  south-central 
Oklahoma  owes  its  origins  to  the  commercialism  of  the 
mineral  springs  and  the  life-sustaining,  freshwater 
artesian  springs. 

The  land  encompassing  these  springs  were  given  to 
the  Choctaw  and  Chickasaw  Indian  Tribes  in  1835  by  the 
U.S.  Government  (Barker  and  Jameson,  1975).  By  1880 
white  ranchers  had  moved  in  by  intermarriage  and  leases, 
and  began  to  take  over  the  economic  and  political  aspects 
of  the  region.  The  Chickasaws  were  fearful  that  they 
would  lose  these  springs  to  private  interests  pending 
statehood.  So,  they  relinquished  control  of  the  springs 
back  to  the  U.S.  Government  in  1902.  This  was  accom- 
plished so  that  all  people  could  enjoy  the  use  of  these 
springs  for  years  to  come. 

The  small  tract  of  land,  approximately  1.5  square 
miles,  along  the  southern  periphery  of  the  city  of 
Sulphur,  became  a  national  park  in  1906.  The  federal 
land  was  named  Piatt  National  Park  in  commemoration  of  a 
United  States  Senator,  Orville  H.  Piatt,  who  served  on 


the  Committee  of  Indian  Affairs. 

The  park  was  enlarged  to  15  square  miles  in  1976  and 
subsequently  renamed,  Chickasaw  National  Recreation  Area. 
Additions  to  the  park  system  included  Arbuckle  Reservoir, 
Veterans  Lake,  and  the  Vendome  Well. 

Charles  Gould,  a  geologist,  forwarded  a  report  to 
the  Secretary  of  the  Interior  in  1906  that  stated  that 
there  were  33  springs  in  the  park.  To  date,  there  are 
only  5  springs  that  flow  with  any  degree  of  regularity. 
A  strong  possibility  exists  that  all  springs  will  cease 
flowing  in  the  next  40  years  unless  groundwater  conserva- 
tion measures  are  taken.  The  increasing  frequency  of 
cessation  of  springflow  can  be  attributed  to  depriva- 
tional  effects  of  production  wells  in  close  vicinity  to 
the  park . 

1.3   Scope  of  the  Project 

The  primary  thrust  of  the  project  will  be  confined 
to  the  hydrogeologic  investigation  of  the  two  freshwater 
artesian  springs  in  the  park,  Antelope  and  Buffalo.  The 
mineralized  springs  will  be  included  as  a  supplementary 
analysis  of  the  structural  and  lithological  parameters 
that  govern  groundwater  movement  in  the  Sulphur, 
Oklahoma,  area. 

The  intentions  of  this  study  are  to  accurately  ad- 
dress the  two  artesian  springs  in  the  form  of  data  analy- 


sis  and  provide  a  concise,  descriptive  narrative  of  the 
same. 

The  principal  objectives  of  this  investigation  are: 
(1)  review  previous  findings  and  literature  pertinent  to 
the  study;  (2)  deduce  the  origin  of  the  two  artesian 
springs;  (3)  assign  concise  nomenclature  to  the  springs; 
(4)  assess  the  geologic  characteristics  of  the  aquifer 
responsible  for  spring  discharge;  (5)  determine  the  re- 
charge area;  and  (6)  utilize  the  Prickett-Lonnquist 
groundwater  model  to  approximate  aquifer  properties,  at 
selected  points. 


CHAPTER  II 


BACKGROUND 


2.1   Geography  and  Hydrology  at  the  Project  Area 

The  area  of  investigation  is  located  in  and  adjacent 
to  Sulphur,  Murray  County,  Oklahoma  (Figure  1).  South 
and  east  of  the  city  of  Sulphur  is  the  Arbuckle  Mountains 
(Figure  2) .  Reference  to  the  Arbuckle  Mountains  in 
south-central  Oklahoma  is  somewhat  misleading,  since  8  0 
percent  of  the  roughly  triangular  area  of  1,000  square 
miles  is  rolling  plains  (Ham,  1978)  .  The  Arbuckle 
Mountains  can  be  further  subdivided  into  western  and 
eastern  regions  (Figure  3) .  The  study  area  is  located  in 
the  western  portion  of  the  eastern  part  of  the  Arbuckle 
Mountains . 

Buffalo  and  Antelope  Springs  are  located  at  the  east 
end  of  Travertine  District,  Chicasaw  National  Recreation- 
al Area  (Figure  4) .  The  remainder  of  the  flowing  springs 
are  located  in  the  central  portion  of  the  park  and  are 
somewhat  mineralized.  Surface  elevations  of  Buffalo  and 
Antelope  are  1078  feet  and  1080  feet,  respectively  (Table 
1)  .   The  maximum  surface  elevation  in  the  park  is  1160 
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Figure  2.   Location  of  Sulphur  in  relation 
to  the  Arbuckle  Mountains  (After 
Johnson,  1984) . 
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Figure  3.   Map  showing  the  eastern  and  western  regions 
of  the  Arbuckle  Mountains  (Modified  after 
Fairchild  et  al. ,  1983). 
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Figure    4.      Map   depicting   location   of    springs    in    the 
Travertine    District   of   Chickasaw  National 
Recreation   Area    (After   National   Park    Service, 
1984)  . 
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Table  1.   Surface  Elevations  of  Springs 

Spring  Surface  Elevation  (feet) 

Antelope  1080 

Buffalo  1078 

Black  Sulphur  94  2 

Pavillion  938 

+Hillside  960 

Bromide/Medicine  925 


Hillside  Spring  is  now  considered  to  be  one  of  the 
Pavillion  Springs  by  the  National  Park  Service.   For 
the  purpose  of  this  study,  the  original  nomenclature  of 
the  spring  will  be  retained. 

Data  obtained  from  the  1984  United  States  Department  of 
the  Interior,  National  Park  Service,  Topographic  Map. 
Contour  interval  is  2  feet. 
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feet,  at  the  northeast  corner.  The  minimum  surface  ele- 
vation is  915  feet,  along  Rock  Creek  at  the  southwest 
corner  of  the  park. 

The  Travertine  District  is  comprised  primarily  of  a 
vibrant  growth  of  hardwoods.  This  particular  ecotone  is 
sustained  by  the  springs,  and  by  the  perennial  streams, 
Travertine  and  Rock  Creek.  Travertine  Creek  traverses 
the  park  in  a  westerly  direction  until  its  confluence 
with  Rock  Creek  at  the  center  of  the  park.  Then,  Rock 
Creek  flows  in  a  southwesterly  direction  and  reaches  its 
base  level  at  Arbuckle  Lake,  surface  elevation  of  872 
feet. 

Outside  the  park,  to  the  north,  the  terrain  consists 
of  rolling  hills  heavily  dissected  with  numerous  ravines 
and  intermittent  streams.  Riparian  vegetation  abounds  in 
the  stream  bottoms,  and  grasses  grow  mainly  on  the  hills. 
The  maximum  surface  elevation  is  approximately  1130  feet. 
Rock  Creek  is  the  major  stream,  and  flows  in  a  southerly 
direction. 

To  the  northeast  and  east  of  the  park,  within  a 
7-mile  radius,  the  terrain  consists  of  a  plains  topog- 
raphy, moderately  dissected  with  intermittent  streams. 
The  maximum  surface  elevation  is  1309  feet.  Numerous 
low-lying  rock  outcrops  of  the  Arbuckle  Group  appear  in 
this  area.  Mill  Creek  is  the  largest  intermittent  stream 
in  this  area  and  flows  in  a  southerly  direction. 
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South  of  the  park  within  a  2-mile  radius,  the  ter- 
rain is  hilly,  heavily  dissected  with  numerous  ravines 
and  intermittent  streams.  The  maximum  topographic  relief 
in  the  entire  study  area  outside  the  park  boundary  occurs 
along  U.S.  Highway  177,  with  a  differential  elevation  of 
156  feet  for  a  distance  of  one  mile.  Riparian  vegetation 
exists  along  the  stream  bottoms,  cacti  and  grasses  grow 
on  the  rocky  ridges  and  hilltops.  A  major  intermittent 
stream,  Wilson  Creek,  flows  in  a  westerly  direction  be- 
fore reaching  its  base  level,  Veterans  Lake.  The  maximum 
surface  elevation  is  1238  feet. 

Southeast  of  the  park  within  a  5-mile  radius,  the 
terrain  is  characterized  by  comparatively  wider  valleys 
and  gently  rolling  hills.  The  south-flowing  Mill  Creek 
is  the  major  intermittent  stream.  The  maximum  surface 
elevation  is  1160  feet. 

Generally,  the  terrain  in  the  study  area  is  charac- 
terized by  a  gently  rolling  topography  overlying  rela- 
tively flat-lying,  fractured  and  faulted  carbonates 
(Fairchild  et  al.,  1983).  However,  the  topographic  re- 
lief in  the  park  is  comparitively  greater  than  that  of 
the  surrounding  study  area.  Also,  just  south  of  the  park 
are  northwest-trending  ridges  that  are  relatively  pro- 
nounced. Both  of  these  topographic  expressions  can  be 
attributed  to  synclinal  deformation  and  faulting.  In 
addition,  karst  topography  is  not  well  developed  in  the 
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project  area  (Fairchild  et  al.,  1983). 

The  study  area  is  located  in  a  moist  subhumid  zone. 
The  climate  is  generally  characterized  by  long  summers 
and  relatively  mild  winters  of  short  duration.  The  mean 
annual  air  temperature  for  Sulphur  is  62.5°F  for  the  time 
period  of  1941  to  1970  (Fairchild  et  al.,  1983).  The 
mean  annual  precipitation  is  37.96  inches  for  the  time 
period  of  1919  to  1984  (Table  2) . 

Maximum  temperatures  commonly  occur  during  the 
months  of  July  and  August  (Fairchild  et  al.,  1983).  Min- 
imum temperatures  occur  during  the  months  of  December  and 
January.  Maximum  precipitation  generally  occurs  during 
the  spring  months,  from  April  through  June  (Table  2). 
Another  wet  season  of  lesser  magnitude  takes  place  during 
September  and  October.  The  dry  season  normally  occurs 
from  November  through  February. 

2.2   Geology  of  the  Area 

The  Arbuckle  Mountains  consist  mostly  of  Paleozoic 
carbonates  and  Late  Paleozoic  elastics  deposited  partly 
upon  a  craton  of  Precambrian  granite  and  partly  in  a 
geosyncline  where  the  basement  rocks  are  Middle  Cambrian 
Volcanics  (Figure  5)  .  The  western  part  of  the  Arbuckle 
Mountains,  named  the  Arbuckle  Anticline,  represents  an 
intense  orogenic  event  of  the  geosyncline,  composed  of 
34,000  feet  of  Paleozoic  sediments.   The  eastern  portion 
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Table  2.   Rainfall  Record 

Chickasaw  National  Recreation  Area 
Sulphur,  Okla. 

For  the  Years  1917-1984 


Year 

Jan 

Feb 

Mar 

Apr 

May 

June 

July 

Aug 

Sept 

Oct 

Nov 

Dec 

TOTAL 

1917 

.62 

.17 

4.68 

2.28 

3.04 

3.12 

2.15 

3.09 

T 

1.72 

T 

1918 

.62 

4.15 

4.37 

2.56 

3.92 

1919 

1.21 

2.39 

2.18 

3.40 

5.35 

3.67 

3.81 

7.25 

.79 

8.04 

3.86 

1.10 

43.04 

1920 

2.65 

.14 

1.45 

3.22 

6.08 

4.88 

1.64 

2.15 

5.38 

5.64 

1.18 

1.48 

35.89 

1921 

1.66 

2.17 

4.91 

4.27 

2.80 

8.00 

3.90 

2.83 

.78 

.05 

1.11 

1.06 

33.54 

1922 

1.40 

1.36 

3.09 

8.31 

6.08 

1.11 

2.58 

.79 

.04 

2.33 

4.50 

.40 

31.99 

1923 

2.51 

1.28 

1.80 

5.30 

8.88 

3.87 

T 

2.51 

6.54 

10.75 

2.81 

4.26 

50.51 

1924 

.89 

.79 

3.42 

5.88 

3.10 

3.63 

.35 

.72 

1.98 

.14 

.49 

2.85 

24.24 

1925 

.99 

.78 

.08 

2.19 

4.83 

.02 

2.95 

1.32 

5.08 

3.63 

1.48 

.32 

23.67 

1926 

3.47 

.72 

2.56 

2.56 

4.40 

3.17 

6.19 

8.28 

5.84 

5.39 

.66 

3.31 

46.55 

1927 

4.78 

2.97 

1.95 

9.99 

.76 

6.44 

7.70 

2.93 

5.00 

3.19 

.77 

3.65 

50.13 

1928 

3.30 

2.99 

1.10 

7.86 

3.20 

9.77 

4.97 

3.06 

.76 

3.55 

2.68 

1.74 

44.98 

1929 

2.04 

1.22 

3.52 

1.87 

10.46 

2.49 

1.84 

.14 

5.05 

6.76 

2.36 

2.78 

40.53 

1930 

2.10 

2.76 

1.23 

3.73 

7.99 

2.57 

1.18 

3.42 

1.73 

2.47 

4.02 

2.34 

35.54 

1931 

.35 

7.20 

4.30 

2.02 

1.83 

1.63 

2.95 

.72 

.14 

6.22 

6.03 

1.44 

34.83 

1932 

7.33 

3.41 

1.73 

3.71 

2.11 

4.06 

1.48 

1.76 

.44 

2.54 

.22 

9.30 

38.09 

1933 

2.50 

2.02 

4.40 

2.25 

13.78 

.54 

2.22 

5.95 

4.96 

.   .46 

1.82 

1.44 

42.34 

1934 

2.40 

2.72 

3.80 

2.97 

3.00 

2.81 

.00 

.87 

6.38 

1.06 

5.98 

.48 

32.47 

1935 

1.89 

.85 

5.38 

4.31 

12.54 

7.48 

2.34 

4.14 

5.45 

3.84 

3.09 

3.05 

54.36 

1936 

.43 

.40 

1.61 

1.07 

8.38 

.60 

2.52 

T 

12.68 

3.56 

.31 

1.62 

33.17 

1937 

2.47 

.18 

4.06 

3.68 

2.22 

2.79 

4.63 

4.38 

.05 

3.27 

2.10 

2.93 

32.76 

1938 

2.60 

9.55 

4.61 

2.35 

6.05 

4.83 

2.57 

1.43 

2.44 

.50 

2.80 

.68 

40.41 

1939 

3.07 

1.73 

1.96 

2.78 

2.70 

4.35 

1.26 

2.67 

.59 

2.91 

2.28 

1.20 

27.50 

1940 

.45 

2.87 

T 

6.82 

9.33 

6.58 

6.83 

2.47 

.22 

3.28 

6.67 

2.87 

48.39 

1941 

2.88 

3.18 

.47 

6.89 

4.62 

7.71 

1.06 

6.63 

4.02 

14.89 

1.59 

1.58 

55.52 

1942 

.46 

2.30 

1.70 

10.60 

4.28 

8.18 

.77 

4.76 

3.68 

5.41 

2.24 

2.53 

46.93 

1943 

.15 

.74 

3.14 

4.23 

7.81 

2.73 

1.21 

T 

1.21 

1.63 

.12 

3.87 

26.84 

1944 

2.78 

4.76 

4.63 

2.75 

5.86 

2.54 

3.12 

2.77 

.85 

4.44 

3.82 

2.65 

38.97 

1945 

1.70 

4.44 

9.76 

8.61 

1.41 

10.99 

4.58 

6.13 

11.13 

.78 

.72 

.08 

60.33 

1946 

6.42 

3.66 

4.01 

2.83 

5.66 

3.43 

.84 

6.17 

2.72 

.14 

6.80 

8.42 

51.64 

1947 

.32 

.34 

1.21 

9.04 

8.04 

4.28 

1.57 

1.13 

3.33 

1.61 

2.87 

2.70 

39.34 

1948 

.98 

4.19 

2.51 

.70 

7.59 

6.64 

3.62 

1.32 

.12 

.78 

.48 

1.27 

30.20 

1949 

5.97 

2.69 

3.38 

2.18 

6.14 

4.90 

.58 

2.75 

6.28 

4.54 

.00 

1.48 

40.84 

1950 

3.02 

1.86 

.35 

2.03 

6.44 

3.39 

5.91 

7.94 

2.48 

.71 

-  .22 

.10 

34.45 

(continued) 
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Table  2.   Rainfall  Record  (continued) 

Chickasaw  National  Recreation  Area 
Sulphur,  Okla. 

For  the  Years  1917-1984 


Year 

Jan 

Feb 

Mar 

Apr 

May 

June 

July 

Aug 

Sept 

Oct 

Nov 

Dec 

TOTAL 

1951 

1.14 

3.61 

1.35 

1.43 

5.88 

6.87 

5.54 

2.92 

1.78 

3.15 

1.84 

.31 

35.82 

1952 

.31 

1.23 

3.33 

4.72 

4.86 

.50 

2.85 

.90 

.10 

.03 

4.47 

1.45 

24.75 

1953 

.55 

.97 

3.14 

4.36 

4.51 

2.28 

7.72 

1.44 

1.16 

5.15 

2.22 

1.05 

34.55 

1954 

1.49 

.55 

.93 

5.45 

8.21 

4.34 

.15 

.83 

.36 

5.72 

.36 

3.08 

31.47 

1955 

1.47 

2.18 

2.45 

1.81 

4.77 

1.79 

2.09 

1.08 

7.37 

.21 

.02 

.40 

25.64 

1956 

1.47 

2.14 

.39 

2.89 

4.81 

1.55 

1.43 

1.39 

.04 

3.73 

3.93 

3.15 

26.92 

1957 

2.13 

2.09 

4.37 

9.14 

13.24 

6.09 

2.34 

1.68 

11.41 

2.25 

4.86 

1.46 

61.06 

1958 

2.44 

.51 

4.15 

2.06 

4.85 

3.52 

4.30 

1.91 

.78 

.33 

1.94 

.86 

27.65 

1959 

.45 

.79 

3.36 

2.74 

4.45 

2.87 

3.04 

2.07 

1.81 

8.96 

2.50 

2.77 

35.81 

1960 

2.74 

3.08 

2.10 

2.76 

6.89 

1.84 

3.59 

3.16 

5.33 

4.38 

.21 

4.95 

41.11 

1961 

.10 

2.00 

3.62 

.92 

3.75 

3.51 

3.70 

.64 

6.87 

3.72 

3.20 

1.74 

33.67 

1962 

.34 

.35 

2.55 

2.97 

1.91 

8.91 

4.13 

1.19 

4.69 

4.85 

2.75 

1.79 

36.43 

1963 

.07 

.04 

5.02 

3.22 

.81 

.56 

2.28 

1.99 

1.77 

.04 

2.29 

1.45 

19.54 

1964 

.96 

1.59 

2.89 

1.79 

8.52 

.98 

.70 

4.98 

6.33 

.47 

6.70 

.72 

36.63 

1965 

2.14 

1.55 

1.65 

2.60 

5.90 

1.69 

1.58 

2.48 

2.04 

1.95 

.83 

1.32 

25.73 

1966 

1.26 

1.87 

1.26 

5.09 

.71 

2.18 

4.65 

5.36 

2.04 

.1.86 

.94 

.69 

27.91 

1967 

.34 

.80 

1.05 

8.17 

6.68 

6.73 

3.89 

.20 

5.81 

5.37 

.79 

2.13 

41.96 

1968 

4.13 

1.27 

3.82 

2.35 

9.47 

4.02 

5.75 

.97 

4.43 

2.87 

5.73 

2.04 

46.85 

1969 

1.25 

2.97 

2.55 

5.21 

5.05 

3.30 

1.74 

1.61 

2.65 

7.92 

.51 

3.47 

38.23 

1970 

.18 

1.74 

3.90 

4.97 

2.26 

4.82 

1.68 

.80 

11.73 

16.43 

.79 

.77 

50.07 

1971 

1.89 

1.15 

.68 

2.76 

2.83 

2.63 

4.05 

4.92 

2.79 

6.62 

.57 

5.46 

36.35 

1972 

.36 

.56 

.61 

3.95 

4.05 

2.39 

1.41 

3.56 

2.52 

7.69 

3.87 

.88 

31.85 

1973 

3.01 

3.00 

6.44 

6.04 

2.72 

8.12 

3.94 

.99 

10.90 

4.89 

8.02 

.86 

58.93 

1974 

.90 

1.76 

1.45 

4.91 

3.64 

3.47 

2.10 

5.49 

8.65 

7.10 

2.48 

1.53 

43.48 

1975 

1.98 

4.26 

4.71 

3. OS 

6.85 

2.43 

4.73 

1.98 

3.55 

1.24 

2.05 

1.50 

38.33 

1976 

.07 

.29 

3.70 

3.88 

7.78 

1.33 

1.90 

1.18 

1.36 

4.71 

.90 

1.57 

28.67 

1977 

1.91 

1.39 

6.39 

3.51 

5.29 

2.48 

2.04 

2.83 

1.89 

1.73 

1.29 

.34 

31.09 

1978 

1.15 

2.82 

3.69 

3.00 

8.12 

3.33 

2.13 

1.95 

2.37 

1.31 

4.46 

.58 

34.91 

1979 

2.33 

1.31 

4.56 

3.68 

4.75 

6.83 

.80 

2.72 

2.63 

3.18 

.95 

1.69 

35.43 

1980 

1.44 

1.26 

.91 

.56 

9.58 

2.50 

.18 

.00 

10.73 

1.74 

1.79 

2.29 

32.98 

1981 

.29 

2.61 

2.52 

2.29 

8.51 

3.46 

3.95 

1.81 

2.10 

16.26 

1.71 

.26 

45.77 

1982 

1.75 

2.29 

1.60 

1.04 

15.71 

6.21 

2.95 

1.04 

1.02 

2.84 

3.43 

1.77 

41.65 

1983 

2.57 

3.70 

2.95 

2.95 

10.67 

3.64 

.90 

3.28 

2.11 

4.53 

2.61 

.57 

40.48 

1984 

.98 

1.98 

3.97 

3.86 

2.83 

4.29 
3.95 

3.51 

2.52 

1.62 

7.55 

4.58 

5.10 

42.59 

Avg. 

2.08 

2.10 

2.80 

3.94 

5.62 

2.79 

2.61 

3.62 

4.02 

2.43 

2.02 

37.96 

(modified  after  Harp  et  al.,  1976) 
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of  the  Arbuckle  Mountains  consists  of  Paleozoic  sediments 
of  a  thickness  of  10,000  feet  upon  the  relatively  more 
stable  craton  (Ham,  1978)  .  Besides  the  marked  differ- 
ences in  thickness  of  sediments  between  the  two  provin- 
ces, major  facies  changes  also  occur  (Figure  6)  .  The 
most  notable  facies  change  occurs  in  the  Arbuckle  Group, 
predominantly  limestone  in  the  geosyncline  and  dolomite 
on  the  craton. 

The   structural   and   depositional   history   of   the 
Arbuckles  can  be  summarized  chronologically  as  follows 
(Ham,  1978)  . 

1.  The  geosyncline  was  filled  with  volcanics  dur- 
ing the  Cambrian. 

2.  Marine  sedimentation  in  the  form  of  carbonates, 
sandstones,  and  shales  continued  from  Late 
Cambrian  through  Mississippian. 

3.  The  early  Pennsylvanian  marked  the  inception  of 
the  emergence  of  the  craton,  sedimentation  of 
elastics  continued  on  the  geosyncline,  and  con- 
glomerates were  deposited  along  the  margins  of 
the  craton. 

4.  The  uplifted  craton  remained  distinct  from  the 
geosyncline  through  Middle  Pennsylvanian. 

5 .  The  comparitively  thinner  sediments  of  pre 
Pennsylvanian  continued  to  erode  on  the  craton 
during   Middle   and   Late   Pennsylvanian   (the 
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basement  rock  of  the  craton  was  6  miles  above 
the  floor  of  the  geosyncline) . 

6.  Compressional  forces  during  the  Upper 
Pennsylvanian  produced  close  folding  and  high- 
angle  thrust  faulting  of  the  thick  prism  of 
sediments  in  the  geosyncline  against  the  rela- 
tively immovable  craton.  The  orogeny  slightly 
affected  elements  of  the  craton. 

7.  Two  orogenic  pulses  yielded  the  older  Ada  Con- 
glomerate and  the  younger  Vanoss  Conglomerate 
(Hart,  1974  and  Ham,  1978)  . 

8.  Additional,  less  pronounced  orogenic  pulses 
continued  through  Lower  Permian,  marking  the 
close  of  the  Arbuckle  Orogeny. 

Erosional  processes  during  and  after  the  orogeny  have 
morphologically  produced  the  surface  geology  of  the 
present-day  Arbuckle  Mountains  (Figure  7)  .  In  addition, 
Figure  7  depicts  the  areal  extent  of  this  hydrogeologic 
investigation.  The  surficial  and  structural  geology  of 
the  project  area  is  depicted  in  greater  detail  in  Plate  I 
(Ham,  1954) .  Included  is  a  geologic  section  of  the  sub- 
surface geology;  based  on  known  stratigraphic  thickness 
and  sequence  of  geologic  units,  and  observed  surface 
structure. 

The  following  geologic  units  in  chronological  order 
appear  as  surface  outcrops  in  the  project  area:   (1)  the 
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Arbuckle  Group  (Upper  Cambrian  to  Lower  Ordovician) ,  (2) 
the  Simpson  Group  (Middle  Ordovician)  ,  (3)  the  Viola 
Limestone  (Upper  Ordovician) ,  (4)  the  Sylvan  Shale  (Upper 
Ordovician) ,  (5)  the  Ada  Formation  (Upper  Pennsylvania) , 
and  (6)  the  Vanoss  Formation  of  the  Pontotoc  Group  (Upper 
Pennsylvanian) .  The  remainder  of  the  geologic  units  in 
the  study  area,  the  Tishomingo  Granite  (Precambrian)  and 
the  Timbered  Hills  Group  (Upper  Cambrian) ,  are  situated 
at  considerable  depth  below  the  surface. 

The  Tishomingo  Granite  is  the  local  basement  rock  of 
the  study  area,  estimated  thickness  of  10  miles.  The 
Tishomingo  Granite  is  generally  coarse  crystalline,  pink, 
with  much  microcline  and  biotite  (Ham  et  al.,  1964)  . 

The  Timbered  Hills  Group  is  comprised  of  the  Reagan 
Sandstone  and  the  Honey  Creek  Formations  (Col.  B,  Figure 
6).  The  Reagan  Sandstone  is  a  feldspathic,  glauconitic 
sequence  deposited  on  the  Tishomingo  Granite  (craton) 
during  Late  Cambrian,  75  to  450  feet  thick  (Ham,  1978) . 
The  Honey  Creek  formation  is  a  fossiliferous ,  sandy 
dolomite,  about  225  feet  thick,  that  lies  conformably 
over  the  Reagan  Sandstone.  The  Timbered  Hills  Group  was 
determined  to  be  deposited  in  a  shallow  marine  environ- 
ment on  a  slowly  subsiding  craton. 

The  Arbuckle  Group  consists  of  the  following  forma- 
tions in  chronological  order:  (1)  Butterly-Royer 
Dolomite  (Late  Cambrian) ,  750  feet  thick;   (2)  McKenzie 
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Hill  Formation  (Early  Ordovician) ,  750  feet  thick;   (3) 
Cool  Creek  Formation  (Early  Ordovician) ,  750  feet  thick; 
(4)   Kinblade  Formation   (Early  Ordovician),   1,000  feet 
thick;   and   (5)   West   Spring   Creek   Formation   (Early 
Ordovician),  750  feet  thick  (Col.  B,  Figure  6). 

The  highly  resistant  Arbuckle  Group,  of  shallow 
water  marine  deposition,  approximate  thickness  of  4,000 
feet  is  primarily  dense,  crystalline  dolomite  with  some 
thin  beds  of  white  friable  sandstones  (up  to  10  feet 
thick)  and  shales  (Col.  B,  Figure  6).  The  Arbuckle  rocks 
are  in  part,  richly  fossiliferous  in  the  Arbuckle  anti- 
cline. Yet,  there  is  poor  faunal  preservation  on  the 
craton  (shelf)  because  of  dolomitization. 

The  Simpson  Group  (Middle  Ordovician),  600  to  1,600 
feet  thick,  is  divided  into  the  following  five  formations 
in  chronological  order:  Joins  Formation,  Oil  Creek 
Formation,  McLish  Formation,  Tulip  Creek  Formation,  and 
Bromide  Formation  (Col.  B,  Figure  6)  .  The  five  forma- 
tions represent  five  complete  sedimentary  cycles  in  a 
marine  environment  with  a  basal  sand  in  each  of  the  four 
upper  ones  (Decker  et  al.,  1931).  The  Tulip  Creek  and 
Joins  are  very  thin  or  lacking  in  the  study  area  and  have 
been  incorporated  into  the  Bromide  and  Oil  Creek  forma- 
tions, respectively,  by  previous  investigators  (Ham  et 
al.,  1954).  Each  formation  of  the  Simpson  Group  in  the 
study   area   basically   contains   greenish-gray   shale, 
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skeletal  calcarenites  or  algal-mat  carbonates;  and  a 
friable,  basal  sand  of  varying  thickness,  50  to  350  feet. 
The  basal  sand  is  composed  of  well  sorted,  rounded,  fine- 
grained frosted  quartz. 

The  Viola  Limestone  (Upper  Ordovician)  of  marine  en- 
vironment is  gray,  fine  to  coarse  grained  limestone,  350 
feet  thick.  The  Viola  lies  unconformably  on  the  Bromide 
Formation  of  the  Simpson  Group  (Col.  B,  Figure  6)  .  The 
Sylvan  shale  (Upper  Ordovician),  150  feet  thick,  lies 
disconformably  on  the  Viola.  The  Sylvan  is  a  dark  green- 
ish-gray shale,  well  laminated,  with  numerous  bedding 
surfaces  that  are  crowded  with  graptolites. 

In  ascending  chronological  order,  the  Ada  Formation 
and  the  Vanoss  Group  (both  units,  Upper  Pennsylvanian) 
represent  two  distinct,  depositional  and  deformational 
episodes  during  the  Arbuckle  Orogeny  (Upper 
Pennsylvanian) .  The  Ada  Formation  primarily  consists  of 
a  well-cemented  conglomerate  of  limestone  cobbles,  red 
shales,  and  bituminous  sandstone.  Thickness  of  the  Ada 
Formation  varies  from  0  to  500  feet  in  the  project  area. 
The  Ada  lies  upon  older  rocks  as  an  angular  unconformity, 
deformation  of  this  strata  occurs  locally.  The  Vanoss 
Group  consists  of  maroon  shales,  arkoses,  and  limestone 
conglomerates.  The  Vanoss,  0-500  feet  thick,  is  unde- 
formed,  and  lies  upon  older  rocks  as  an  angular  unconfor- 
mity.  For  the  purpose  of  this  study  only  and  for  simpli- 
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fication,  the  Ada  and  the  Vanoss  will  be  summarily 
joined;  to  be  called  the  following  hydrogeologic  unit, 
the  Pontotoc  Group. 

The  central  part  of  the  Arbuckle  Mountains,  covered 
by  this  investigation,  encompasses  the  following  major 
structural  features:  (1)  the  Hunton  Anticline,  (2)  the 
Hickory  Syncline,  and  (3)  the  Sulphur  Syncline  (Figure  7 
and  Plate  I) .  The  Hickory  Syncline  is  a  localized  struc- 
tural unit  and  subsidiary  to  the  Hunton  Anticline,  vicin- 
ity of  Hickory,  Oklahoma.  In  actuality,  the  Hickory 
Syncline  is  composed  of  two  synclinal  structures,  one 
north  of  Hickory  and  one  south  of  Hickory  (Section,  Plate 
I)  .  A  major  fault,  trending  in  a  west,  southwest  direc- 
tion, subtends  both  synclines.  The  Simpson  Group  has 
been  preserved  in  both  synclines  from  erosion  by  down 
faulting  and  accompanying  drag  folds  (Ham,  1945)  .  East 
of  Sulphur  at  the  eastern  fringe  of  the  project  area,  is 
the  northwest-trending  Hunton  Anticline  (Figure  7)  .  The 
basement  rock  is  approximately  3,200  feet  below  the  sur- 
face in  this  area  (Figure  8) .  The  relatively  flat-lying, 
highly-resistant,  Upper  Arbuckle  Group  outcrops  1  mile 
east  of  the  Travertine  District  boundary  and  continues 
eastward. 

The  west,  northwest-trending  Sulphur  Fault,  south- 
east of  Sulphur,  is  the  major  fault  in  the  project  area. 
The   fault   is   the   northern  boundary   of   the   Sulphur 
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EXPLANATION 


Contour  interval  2000  feet 
(1000-foot  contours   shown  by  dotted  lines  in  a  few  areas  of  low  relief) 


Control  wells  showing  elevation,  sea-ami  datum,  at 
top  of  basement  rocks 
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DEPICTS   EXTENT  OF 
THE   PROJECT  AREA 


Confirmed  by  examination  of  samples 
and/or  electric  log 


Reported  well,  not  confirmed  by  exami- 
nation of  samples  or  electric  log 


Major  fault,  normal  or  reverse 
Hachures  on  downthrown  aide 


Top  of  basement  rocks  projected  from 
stratigraphic  and  structural  data  in  deep 
veil* 


ea 


Outcrop  Oof 


t  rocks 


Figure    8.       Surface     contour   map   of    the   basement   rock 

of   the  Arbuckle  Mountain   area.         The   solid 
lines,    rectangular   in   shape,    outline   the 
general   extent  of   the  project   area,    vicintiy 
Sulphur,    Oklahoma    (Modified  after  Ham,    1964)  , 
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Syncline,  and  delineates  the  Arbuckle  Group  to  the  north 
(upthrown  side)  ,  from  the  Simpson  Group  preserved  in  the 
synclinal  graben  to  the  south  (Section,  Plate  I)  .  The 
amount  of  high-angle  throw  along  the  Sulphur  Fault  in- 
creases in  the  same  direction  as  the  west,  northwest 
plunging  Sulphur  Syncline. 

The  Sulphur  Syncline  begins  plunging  to  the  north- 
west in  Sections  8  and  17  of  Township  IS,  Range  3E;  with 
moderate  dips  of  10  degrees  northwest  (Plate  I)  .  The 
south  limit  of  the  graben  is  bounded  by  an  unnamed  fault. 
The  fault  decreases  in  throw  in  the  same  direction  as  the 
plunging  syncline.  This  would  account  for  the  steep  dips 
of  the  Simpson  on  the  down  thrown  side  to  the  north  and 
northeast  in  Sections  11  and  12,  Township  IS,  Range  3E . 
According  to  Ham  (1945)  ,  strictly  synclinal  conditions 
prevail  just  south  and  southeast  of  Sulphur.  The  syn- 
cline plunges  under  the  Upper  Pennsylvanian  cap  rock  to- 
wards Sulphur,  completely  encompassing  the  Travertine 
District.  The  graben  is  flanked  to  the  north  by  the 
west-northwest  trending  Sulphur  Fault,  exact  location 
unknown. 

In  general,  the  major  structural  features  of  the 
project  area  trend  in  a  west,  northwest  direction.  The 
major  faults  parallel  the  synclinal  and  anticlinal  axes. 
The  resistant  Arbuckle  Group,  relatively  flat  lying,  out- 
crops to  the  northeast  and  east  of  Sulphur.   The  Simpson 
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Group  is  preserved  in  grabens,  the  Sulphur  Syncline  and 
the  Hickory  Syncline.  The  Pontotoc  Group,  outcrops 
throughout  the  entire  park,  capping  the  northwest- 
plunging  Sulphur  Syncline.  The  Pontotoc  covers  the  re- 
mainder of  the  project  area,  gently  dipping  to  the  north- 
west. 

2.3   Groundwater  Flow  Through  Geologic  Materials 

Introduction  of  the  hydrologic  cycle  is  instrumental 
in  understanding  the  processes  to  the  mode,  occurrence, 
and  origin  of  water  (Figure  9) .  The  word,  cycle,  explic- 
itly connotes  no  end  or  start  point.  The  hydrologic 
equation,  a  succinct  statement  of  the  law  of  mass  conser- 
vation, for  a  specific  period  of  time,  quantifies  the 
hydrologic  cycle  in  the  form  of  the  following  equation 
(Fetter,  1980) . 

Inflow  =  Outflow  ±  Changes  in  Storage 

I'    =   .  0'    ±      S'  (2.1) 

For  example,  inflow,  I',  could  be  represented  by  precipi- 
tation, surface  flow,  groundwater,  seepage,  and  overland 
flow.  Storage,  S1,  could  be  in  the  form  of  a  lake  or 
reservoir.  Outflow,  0' ,  could  be  evapotranspiration,  an 
outlet  stream,  and  seepage  from  the  reservoir  or  lake. 
If  outflow  exceeds  inflow,  storage  will  decrease. 
Conversely,  if  inflow  exceeds  outflow,  storage  will 
increase.   The  hydrologic  equation  can  be  applied  to  any 
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Volcanoes 


Figure  9.   Hydrologic  cycle  (After  Fetter,  1980) 
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size  system  and  is  time  dependent  (Fetter,  1980). 

For  the  purpose  of  this  investigation,  the  hydro- 
logic  budget  will  be  defined  as  inflow  equals  outflow 
(Fairchild  et  al.,  1983).  Inflow  is  the  amount  of  annual 
rainfall  that  falls  on  a  specific  area.  Outflow  entails 
evapotranspiration,  surface  flow,  and  groundwater;  all  of 
which  result  from  precipitation  events  for  the  calendar 
year  (a  set  period  of  time) . 

The  outflow  element  of  the  hydrologic  budget, 
groundwater,  becomes  the  inflow  parameter  of  the  subsur- 
face hydrologic  equation.    The  groundwater   (units   in 

3 
length  )  introduced  to  the  aquifer  is  called  recharge. 

Recharge  is  the  fractional  amount  of  the  annual  rainfall 

(units  in  length)  that  infiltrates  to  a  water  table,  dis- 

2 
tributed  over  an  identified  area  (units  in  length  ) .   The 

aquifer,  which  stores  and  transmits  groundwater,  is  the 
storage  element  of  the  hydrologic  equation.  Changes  in 
the  potentiometric  level  of  the  aquifer  for  a  specified 
time  and  area  are  functions  of  inflow  and  outflow.  Out- 
flow entails  natural  and  artifical  discharge.  Natural 
discharge  includes  springs  and  seeps,  and  subsurface  flow 
from  the  local  system  into  a  regional  one.  Artifical 
discharge  consists  of  domestic,  municipal,  and  industrial 
production  wells. 

Darcy's  Law  and  the  continuity  principle  (law  of 
mass  conservation)  are  fundamental  in  the  development  of 
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the  physics  of  groundwater  flow  (Wang  et  al.,  1982).  The 
continuity  principle  is  a  mathematical  formulation  (a 
partial  differential  equation)  of  the  hydrologic  equation 
for  aquifers,  confined  or  unconfined.  Darcy's  Law  re- 
lates the  specific  discharge  vector  to  the  gradient  of 
the  potential,  and  is  normally  valid  for  a  Newtonian  flu- 
id (such  as  water)  and  nonturbulent  flow  (Narasimhan, 
1982)  .  Laminar  or  nonturbulent  flow  for  a  porous  media 
is  generally  considered  to  be  less  than  a  Reynolds  Number 
of  1  (Collins,  1961) . 

Henry  Darcy,  a  French  engineer,  in  1856,  conducted 
an  experiment  (Figure  10)  of  water  flow  through  a  verti- 
cal column  of  sand  and  arrived  at  the  following  equation: 

h   -  h 
Q  =  -KA  — (2.2) 

L2  -  Ll 

Where  volume  discharge  rate,  Q,  is  proportional  to  the 

head  drop  (h2  -  h- )  to  a  cross-sectional  area,  A,  and 
inversely  proportional  to  the  length  difference  (L^  - 
Lx)  . 

The  negative  sign  indicates  that  groundwater  flows 
in  the  direction  of  head  loss  (Wang  et  al.,  1982). 
Finally,  a  proportionality  constant  of  a  geologic  materi- 
al, K,  known  as  hydraulic  conductivity,  is  included  in 
the  right  side  of  Equation  1.1.  Dimensional  analysis  of 
K  reveals  units  of  LT  .  The  general  form  of  Darcy's  Law 
is  given  by  the  following  differential  equation: 
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Figure  10.   A  sketch  of  an  apparatus  to  demonstrate 

Darcy's  Law  (After  Davis  and  Dewiest,  1966) 
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Figure  11.  Approximate  range  of  hydraulic  conduc- 
tivity in  soil  and  rock  (After  Houstin 
et  al. ,  1982) . 
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Q       dh 
q=-=-K  —  (2.3) 

A       dl 

Where  specific  discharge  or  apparent  velocity,  q,   is 

equal  to  the  hydraulic  conductivity,  K,  times  the  gradi- 

ent,  -=y  (Davis  and  Dewiest,  1966)  .  Also,  Figure  10  pic- 
torially  shows  that  hydraulic  head,  h,  is  a  function  of 
elevation  ,  z,  from  some  arbitrary  datum,  and  pressure 
head,  P/y/  where  the  pressure,  P,  is  the  pore  pressure, 
and  y,  is  the  specific  weight.  Hydraulic  head  in  these 
equations  is  energy  per  unit  weight  of  water. 

The  total  energy  of  an  incremental  volume  of  fluid 
in  motion  is  equal  to  the  sum  of  kinetic,  gravitational, 
and  fluid-pressure  energy: 

Etv  =  1/2  pv2  +  pgz  +  P  (2.4) 

where  E  is  total  energy  per  incremental  volume,  p  is 
density,  g  is  the  gravitational  acceleration  constant, 
and  P  is  pressure.   If  equation  (1.3)  is  divided  through 

by  density,  the  result  is  total  energy  per  unit  mass: 

2 
E    =  v  /2  +  gz  +  p/p  =  constant  (2.5) 

which  is  the  Bernoulli  equation. 

The  kinetic  energy  term,  velocity  squared,  is 
generally  considered  to  be  negligible  in  flow  through  po- 
rous media  and  is  dropped  from  the  above  equation.  The 
force  potential  is  expressed  in  terms  of  energy  per  unit 
mass  (Fetter,  1980) : 
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pgz 
$  =  gz  +  P/p  =  gz  +  (2.6) 


and  h  =  z  +  h  ,  the  hydraulic  head,  and  h_  is  feet  of 
p        J  P 

pressure,  then 

$  =  gh  (2.6) 

Both  force  potential  and  hydraulic  head  are  poten- 
tials. Darcy's  Law  can  be  expressed  in  hydraulic  head 
and  also  in  force  potential  for  one  dimensional  flow: 

q  -  A  -    K  dL  U.I) 

Q  _    K   d$  n    R. 

*  "  A  ~  "  g  dL  (2*8) 

Hydraulic  conductivity  is  not  only  a  function  of  charac- 
teristics of  the  medium,  it  depends  on  the  properties  of 
the  fluid  as  well.  The  hydraulic  conductivity  is  given 
by  the  following  equation: 

K  =  K1  pg/y  (2.9) 

where  K  is  intrinsic  permeability  and  is  characteristic 
of  the  medium  only,  and  y  is  the  dynamic  viscosity  of  the 
fluid  (Davis  and  Dewiest,  1966).  Intrinsic  permeability 
is  proportional  to  the  average  pore  diameter: 

K.  =  cd2  (2.10) 

where  c  (dimensionless  units)  is  a  shape  factor  which  ac- 
counts for  stratification,  packing,  arrangement  of 
grains,  size  distribution,  and  porosity. 

Porosity  of  a  porous  media,  n,  is  a  fraction  of  the 
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bulk   volume,    V    ,    occupied   by   voids    (Collins,    1961) 
B 

V  Volume  of  Pores        1-  V        1-  Volume  of  Solids 

n  =  —  =     =  —  =       (2.11) 

V  Bulk  Volume        V         Bulk  Volume 
B  B 

where  bulk  volume  is  equal  to  the  sum  of  the  volume  of 
pores,  V  ,  and  the  volume  of  solids,  V  . 

P  o 

Porosity  of  a  porous  material  is  a  function  of  grain 
size  distribution,  grain  configuration,  and  diagenetic 
processes.  Primary  porosity  is  related  to  pore  spaces 
between  grains.  Secondary  porosity  is  related  to 
fractured  media  and  subsequent  enlargement  by  solution 
activities  of  groundwater.  Effective  porosity  is  the 
fractional  amount  of  interconnected  pores  for  a  given 
bulk  volume  through  which  groundwater  can  be  transmitted. 
Absolute  porosity  is  the  total  fractional  amount  of  voids 
for  a  given  bulk  volume,  regardless  of  its  ability  to 
transmit  fluids  (Chilingarian  et  al.,  1975). 

In  general,  hydraulic  conductivity  is  a  measure  of 
the  ability  of  porous  media  to  transmit  groundwater  and 
is  related  to  effective  porosity.  As  interconnected  pore 
diameters  become  smaller,  resistance  to  flow  increases, 
resulting  in  a  decrease  in  hydraulic  conductivity  and  a 
corresponding  increase  in  hydraulic  gradient.  Converse- 
ly, hydraulic  conductivity  approaches  infinity  in  large 
solution  cavities  and  fractures  (Chilingarian  et  al., 
1975)  .  A  wide  spectrum  of  values  for  hydraulic  conduc- 
tivity of  soil  and  rock  are  shown  in  Figure  11  (Houstin 
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et  al. ,  1982)  . 

Values  of  porosity  vary  considerably  for  sedimentary 
rocks.  elastics  generally  range  from  .03  to  .30.  Carbon- 
ates range  from  less  than  .01  to  about  .30  (Fetter, 
1980) .  Table  3  lists  typical  porosities  of  porous  mate- 
rials (Houstin  et  al.,  1982).  Indurated  sediments  at 
depths  of  several  thousand  feet  might  have  porosities 
that  are  only  half  to  one-third  the  porosities  of  the 
same  lithology  of  the  surface  or  at  shallow  depths 
(Goodman,  19  79)  . 

The  ensuing  discussion  of  aquifer  characteristics  is 
necessary  prior  to  the  development  of  the  continuity 
equation,  and  subsequent  methodology  of  hydrogeologic  ex- 
ploration. 

An  aquifer  is  a  geologic  unit  that  stores  and  trans- 
mits water  rapidly  enough  to  supply  adequate  amounts  to 

wells.   Also,  a  porous  media  with  at  least  an  intrinsic 

-2 
permeability   of   10     darcy   constitutes   an   aquifer 

(Fetter,  1980)  . 

The  intrinsic  permeability  of  a  confining  layer  is 

-2 
generally  less  than  10    darcy.   The  confining  layer  can 

be  further  subdivided  into  aquifuges  (impermeable) , 
aquitards  (low  permeability)  ,  and  aquicludes  (upper 
boundary  or  a  lower  boundary  to  a  flow  system,  low  per- 
meability) .  Confining  layers  of  large  aerial  extent, 
hundreds  of  square  miles,  can  provide  significant  amounts 
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Table  3.   Typical  Porosities  for  Porous  Materials 


Soft  skeleton 


Material     Range  in  n    Commonly  en- 

(%)  countered 

values  of  n  (%) 


Stiff  skeleton 


Sand 

25 

- 

50 

40 

+ 

8 

Silt 

40 

- 

63 

52 

+ 

8 

Clay 

40 

- 

94 

62 

+ 

12 

Marble 

0.3 

- 

1.1 

0.7 

+ 

0.4 

Granite 

0.1 

- 

20 

2 

+ 

1 

Limestone 

0.4 

- 

43 

10 

+ 

20 

Sandstone 

0.7 

- 

34 

15 

+ 

79 

Shale 

1.6 

- 

25 

15 

+ 

7 

Tuff 

14 

- 

40 

25 

+ 

6 

Chalk 

25 

— 

50 

35 

+ 

7 

(after  Houstin,  1982;  and  originally  from  Clark,  1966; 
Riley,  1972;  and  Goodman,  1979) 
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of  recharge  to  the  underlying  aquifer. 

Water-table  aquifer  (also  referred  to  as  phreatic  or 
unconfined  aquifer)  is  defined  as  a  continuous  geologic 
material  of  high  intrinsic  permeability  from  the  land 
surface  to  the  base  (Figure  12)  .  Recharge  is  normally 
downward  seepage  from  the  unsaturated  zone. 

Confined  aquifers  or  artesian  aquifers  are  overlain 
by  confining  layers.  Recharge  occurs  from  confining  lay- 
ers or  where  the  aquifer  outcrops  (Figure  12) .  The 
piezometric  level  of  an  artesian  aquifer  is  above  the 
base  level  of  the  overlying,  confining  layer.  The  dis- 
tance above  the  base  level  of  the  confining  layer  to 
which  water  will  rise  is  the  pressure  head  of  the  force 
potential. 

Transmissivity  is  the  measure  of  the  amount  of  water 
that  is  transmitted  laterally  by  the  full  saturated 
thickness  of  an  aquifer.  Transmissivity,  T,  is  the  prod- 
uct of  hydraulic  conductivity,  K,  and  saturated  thick- 
ness, h1  : 

T  =  Kh'  (2.12) 

2 
For  this  study,  units  for  transmissivity  are  ft  /day. 

Storage  coefficient,  S  ,  is  the  volume  of  water  than 

an  aquifer  can  expel  or  absorb  from  storage  per  unit  area 

per  unit  change  in  head  (Fetter,  1980) .   Water  is  stored 

or  withdrawn  from  a  confined  aquifer  by  elastic  behavior 

of  the  media.   Water  extracted  or  absorbed  in  a  phreatic 
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aquifer  is  a  function  of  gravity  and  elastic  mechanisms. 
Specific  storage,  S  ,  is  the  amount  of  water  released  per 
unit  volume  of  aquifer  due  to  compressibility  of  the 
geologic  unit  and  compressibility  of  the  water: 

Ss  =  y(<*  +  nB)  (2.13a) 

Y   =  pg  (2.13b) 

and 

Ss  =  pg  (a  +  nS)  (2.14) 

where      p  =  density  of  water, 

g  =  acceleration  of  gravity, 

Y  =  specific  weight  of  water, 

a  =  matrix  compressibility  of  the  aquifer, 

n  =  porosity, 

6  =  compressibility  of  the  water. 

The  storage  coefficient  of  a  confined  aquifer  is  the 

product  of  specific  storage,  S  ,  and  the  aquifer  for 

thickness ,  b' : 

S   =  b'  Sc  (2.15) 

c        S 

The  storage  coefficient  of  a  phreatic  aquifer  is  sum  of 

the  specific  yield,  S  ,  and  the  specific  storage,  S  : 

S   =  S   +  h'S^  (2.16) 

c    y      S 

where  h1  is  the  saturated  thickness  of  a  phreatic 
aquifer.  It  is  customary  to  set  the  storage  coefficient 
equivalent  to  the  specific  yield,  because  specific 
storage  is  at  least  two  orders  of  magnitude  less  than 
specific  yield  (Fetter,  1980).   Values  of  specific  yield 
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for  an  unconfined  aquifer  range  from  .02  to  .30.  Storage 
coefficients  of  a  confined  aquifer  are  usually  less  than 
.005. 

A  homogeneous  aquifer  has  the  same  rock  properties 
throughout  the  porous  media.  Narasimhan  (1982)  stated 
that  homogeneity  is  a  frequency  distribution  of  rock 
properties  that  does  not  vary  spatially  for  a  geologic 
unit.  Conversely,  a  heterogeneous  media  is  characterized 
by  frequency  distributions  of  rock  properties  that  vary 
spatially. 

If  all  measurable  properties  of  a  media  at  some  fi- 
nite element  or  point  are  spatially  equal  the  unit  is 
isotropic.  If  a  certain  property  of  the  media  at  some 
finite  element  has  a  preferred  orientation  and  is  not 
spatially  equal,  the  unit  is  anisotropic.  Examples  of 
anisotropic  behavior  of  intrinsic  permeability  of  a  media 
are:  (1)  fractured  carbonates,  and  (2)  preferred  grain 
orientation  in  elastics. 

The  derivation  of  the  law  of  mass  conservation  (the 
continuity  principle)  for  an  isotropic,  homogeneous,  ele- 
mental volume  (Figure  13)  of  a  confined  aquifer  is  given 
in  Appendix  I.  The  resultant  continuity  equation  of  two- 
dimensional  transient  flow  is: 

82h      82h      S    3h 

2   +   n   =   -°   —  (2.17) 

8x       8y       T    8t 

If  the  confined  aquifer  is  anisotropic,  (where  T   ^ 
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•  •  Water.'tab'le.'j'f'i^j^ Lake ^ 


Recharge 
area 


PotentJom_etric_s^^^ 


Figure  12.   A.  Water-table  aquifer 
B.  Artesian  aquifer 
(After  Fetter,  19  80) . 


Figure  13.   Elemental  volume  of  flow  through 

a  confined  aquifer  (After   Fetter, 
1980) . 
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T  )  ,  and  heterogeneous  (where  T  and  T  vary  spatially 
within  the  aquifer)  then,  the  following  two-dimensional, 
transient  flow  equation  applies: 

3      3h      3      3h         3h 

—  T   —   +   —  T   —   =   S   —  +  Q  (x,y,t)      (2.18) 

3x   X  3x      3y   Y  3y       C  3t     n 

where  +Q   per  unit  area  denotes  a  sink,  and  -Q   repre- 
n  r  n    r 

sents  recharge  per  unit  area  (Wang  et  al.,  1982).  There 
is  no  known  analytic  solution  to  Equation  (2.18).  As  a 
result,  numerical  methods  are  necessary  to  solve  the 
preceding  equation.  A  well  documented  and  widely  used 
groundwater  model,  Prickett  and  Lonnquist  (1971),  can  be 
utilized  to  solve  Equation  2.18.  In  addition,  the 
Prickett-Lonnquist  model  can  solve  phreatic,  two- 
dimensional  transient  flow: 

3     3h      3     3  h       3  h 

—  h  —  +  —  h  —  =   S   —  (2.19) 

3x    3x      3y    3y      y  3t 

which  is  the  Boussinesq  equation  (Fetter,  1980)  .  Appen- 
dix 2  contains  the  Prickett-Lonnquist  model,  and  an  ac- 
companying flow  chart  and  discussion. 

The  remainder  of  this  section  will  be  devoted  to 
anisotropic  flow  behavior.  Before  proceeding  into  this 
topic,  an  understanding  of  isotropic  flow  is  necessary. 

For  isotropic  porous  media  and  steady  state 
conditions,  the  specific  discharge,  q,  flows  parallel  to 
the  gradient,   and  normal   to  the  equipotential   lines 
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* 
(Figure   14  ) .   A  potentiometric  map  constructed   from 

field  data,  consists  of  isobars  of  equipotential  points 

that   depict  water   surface   contour   elevations   (force 

potential,  $)  .   If  hydraulic  conductivity  varies  in  the 

three  principal  axes,  then  specific  discharge  will  flow 

obliquely  to  the  equipotential  lines  (Davis  and  Dewiest, 

1966)  . 

For  this  study,  anisotropic  flow  behavior  will  be 

specifically   focused   on   fractured  media.    J.E.   Gale 

(1982)   states   that   fracture   hydrogeology   is   poorly 

nurtured  and  not  well  understood  when  compared  to  fluid 

movement  through  porous  media. 

Gale  (1982)  maintained  that  the  hydraulic  charac- 
teristics of  fractured  media  should  be  a  function  of:  the 
permeability  tensor,  K.  .,  which  reflects  the  nature  of 
the  aperture  distribution;  the  effects  of  the  stress 
tensor,  c.  .,  on  fracture  permeability;  and  the  geometry 
of  the  different  fracture  sets  and  the  degree  of  inter- 
connection. He  also  agreed  with  Castillo  (1972)  that  the 
concept  of  flow  through  porous  media  is  inadequate  to  de- 
scribe flow  behavior  in  jointed  rock. 

Hydraulic  conductivity,  K-,  of  a  single  fracture 
through  parallel  plates  for  a  single-phase,  laminar  flow 


* 

Figure  14  is  a  graphical  solution  to  the  Laplace  equa- 

2 
tion,  V  h  =  0. 


43 


of  a  viscous  incompressible  fluid  is: 

pg      o 

K   =  (2b)z  (2.20) 

1         12y 

where  2b  is  the  fracture  aperture  (Gale,  1982).  The  spe- 
cific discharge,  q,  and  continuity  equations  for  inter- 
connected, continuous  fractures  are,  respectively: 

3h 
q  =  -  K,  —  (2.21) 


and 


f  » 


3        8h 


—   K   —   =0  (2.22) 

where  h  is  the  hydraulic  head  and  steady  state  conditions 
prevail . 

In  addition,  empirical  models  were  developed  to  de- 
scribe variations  of  fracture  permeability  with  normal 
stress.  One  such  model  by  Snow  (1968)  ,  that  accounts  for 
normal  stress  phenomena  of  hydraulic  conductivity  of 
parallel  plate  flow  is  (Gale,  1982)  : 

Kh  =  Kq  +  [Kn  (2b)2/Sf]  (P  -  Pq)  (2.23) 

where 

K,  =  horizontal  hydraulic  conductivity  of  fractures, 

K   =  hydraulic  conductivity  at  initial  pressure  P  , 

K   =  normal  stiffness  factor  of  the  fracture, 
n 

S-  =  fracture  spacing, 
P   =  pressure. 
Another  empirical  model  that  accounts  for  pressure,  by  P. 
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Jones  (1975)  ,  for  carbonate  rocks  is: 
K  =  Kq  [log  (Ph/P)]3 

where  K   is  a  constant,  and  K,  hydraulic  conductivity,  is 
o 

equal  to  zero  at  P,  ,  the  "healing"  effective  pressure 
(Gale,  1982) . 

A  basic  flow  equation  for  multiple  fracture  systems 
was  developed  by  Snow  (1965) .  The  flow  of  a  single  frac- 
ture is: 

(2b)3 
q   = pg  (6.  .  -  M.  .)  I.  (2.24) 

3       12y        1D     1D 

where     q .   =   fracture-flux  vector  per  unit  width  of 

fracture  [L3/T] , 

6.  .  =  Kronecker  delta, 
ID 

M.  .  =   n.  n.,  a  matrix  formed  by  direction 
13      i   D 

cosines  of  the  normal  to  the  fracture, 
I.  =  potential  gradient  vector. 
Snow  maintains  that  for  each  fracture  that  intersects  a 
sample  line  has  it  image  at  a  distance  L  equal  to  the 
length  of  the  sample  line  and  in  the  same  direction  as 
that  of  the  sample  line.  Hence,  the  spacing  between  the 
fracture  and  its  repetition  is  as  follows: 

w  =  L  |n.D. |  (2.25) 

where  w  =  spacing  between  the  fractures, 
L   =  length  of  the  sample  line, 

n.  =  direction  cosines  of  normal  to  fracture 

l 

plane , 
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D.  =  direction  cosines  of  sampling  line. 
1 

Thus,  the  permeability  of  a  medium  for  a  single 

fracture  is  given  by: 

1      b3 

K    =—  ~ 7    («ij  "  M±j)  (2.26) 

1J    3L   |n.D.|     J      J 

2 
where  K. .  is  the  intrinsic  permeability  tensor  (L  ) .   The 

permeability  of  a  given  station  is  the  sum  total  of  each 

individual  fracture  (Gale,  1982). 

Other  important  factors  put  forth  by  Gale  (1982) 
are:  (1)  a  strong  correlation  exists  between  surface 
fracture  geometry  and  subsurface  fracture  geometry;  (2) 
fracture  porosity  decreases  logarithmically  with  depth; 
(3)  fractures  have  highly  nonlinear  deformation  charac- 
teristics; and  (4)  fracture  interconnection  with  other 
fractures  determines  to  a  great  extent  the  rock  mass  per- 
meability. 

Gale  (1982)  intentionally  restricted  his  discussion 
on  fracture-flow  phenomenon  to  crystalline  rock  (igneous 
and  metamorphic  rocks) .  He  cited  at  the  beginning  of  his 
article  that  fractured  carbonates  are  hydrologically  dis- 
tinct from  other  fractured  rocks,  because  fluid  flow  is 
dominated  by  solution  channels. 

Duguid  (1973)  developed  a  system  of  equations  that 
is  coupled  by  the  interaction  of  the  fluid  out  of  the 
primary  pores  and  into  the  fractures.  The  interaction 
term  describes  the  mass  flux  from  the  pores  into  frac- 


46 


tures  by  a  pressure  differential.  Gale  (1982)  asserts 
that  Duguid's  theoretical  approach  to  non-steady, 
fractured  porous  media  flow  is  the  most  promising  devel- 
opment as  of  late.  The  system  of  equations  derived  by 
Duguid  (1973)  are: 


3P-  ap9   r   k 

(1  -  n2)  n1  6 +  (1  -  n2)  n2  6  — -  + VZP1  =  0 

3t  3t     p     y 


8P1  3P~   r 

(1  -  nx)  n1  6 +  (1  -  n1)  n2  6  — +  V  •  (q2g)  =  0 

3 1  3 1     p 


(q2s)  =  -(K2/y)  ■  [p  (3(q2s)/3t)  +  VP2]  (2.27) 

where 


n1  =  primary  porosity  or  matrix  porosity, 

n_  =  fracture  porosity, 

6  =  coefficient  of  compressibility  of  water, 

P..  =  pore  pressure  in  the  matrix, 

P~  =  pore  pressure  in  the  fracture, 

r  =  mass  flow  rate  from  block  to  fracture, 

K.  =  intrinsic  permeability  of  the  matrix, 

K_  =  intrinsic  permeability  along  the  fracture, 

q2  =  space-averaged  velocity  over  the  fractures 
relative  to  the  matrix, 

Duguid  developed  a  numerical  solution  to  the  preced- 
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ing  equations  by  using  the  Galerkin  finite  element  method 
(Duguid  et  al.,  1974).  Duguid  et  al .  (1977),  used  the 
preceding  numerical  method  to  solve  a  hypothetical, 
leaky,  confined  aquifer.  A  very  interesting  result  re- 
ported by  Duguid  et  al.  (1977) ,  was  that  discharge  (con- 
stant discharge  problem)  is  divided  between  the  primary 
pores  and  fractures  in  the  ratio  of  matrix  permeability 
to  fracture  permeability. 

Gringarten  (1982)  summed  it  up  quite  well  as  to  when 
usage  of  deterministic  models  and  homogeneous  models  are 
appropriate.  He  states  that  deterministic  models  are 
based  on  duplication  of  the  geometry,  and  discontinuous 
behavior  of  the  fracture  system;  and  are  useful  for 
small-scale  geotechnical  investigations.  Whereas,  ho- 
mogeneous models  that  incorporate  anisotropic  behavior, 
are  useful  for  a  continuum  of  fractured  media  on  a  re- 
gional basis,  such  that  fracture  distributions  are  not 
too  dispersed. 

Karst  hydrogeology  is  somewhat  more  complicated  than 
fracture  flow  because  conduits  are  introduced  into  the 
overall  flow  regime.  Basmaci  (1977)  maintains  that  so- 
lution cavities  and  the  like,  constitute  another  class  of 
porosity,  and  should  be  differentiated  from  primary  and 
fracture  porosities.  In  addition,  he  states  that  solu- 
tion cavities  behave  as  pressure  conduits  or  seepage 
faces,  and  should  be  treated  as  boundaries  to  flow  of  wa- 
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ter  in  karst  regions. 

Figure  15  depicts  the  mean  hydraulic  conductivity 
and  corresponding  distribution  for  a  given  volume  (after 
Basmaci,  1977).  Mean  values  of  hydraulic  conductivity  of 
laboratory  samples  are  the  smallest.  Mean  values  ob- 
tained from  drilled  holes  are  2  orders  of  magnitude  larg- 
er than  laboratory  samples.  Mean  values  for  a  basin  are 
3  orders  of  magnitude  larger  than  drilled-hole  values. 
The  high  values  for  a  basin  are  a  function  of  the  karsti- 
fication  (development  of  solution  cavities) . 

In  summary,  an  understanding  of  the  physics  of 
groundwater  flow  is  crucial  in  properly  assessing 
groundwater  regimes.  In  addition,  the  rate  of  water  flow 
and  skeletal  deformation  of  porous  materials  are  contin- 
gent on  the  mechanical  behavior  of  the  following  main  pa- 
rameters: (1)  porosity,  (2)  permeability,  and  (3)  com- 
pressibility (Houstin  et  al.,  1982). 

2.4   Spring  Flow  Assessment  and  Related  Methodology 

Natural  discharge  from  groundwater  occurs  in  the 
form  of  springs  and  seeps.  A  classification  based  on 
spring  discharge  is  given  in  Table  3  (after  Davis  and 
Dewiest,  1966)  .  A  spring  can  be  further  defined  by  the 
following  parameters:  type  of  aquifer,  chemical  charac- 
teristics, water  temperature,  direction  of  water  mi- 
gration, relation  to  topography,  and  geologic  structure. 
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Figure  14. 


Equipotential  lines,  h-,  through  hg  , 
in  a  two-dimensional  flow  field  and 
the  gradient  of  h  (After  Fetter,  19  80) 
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Figure    15. 


A  schematic  representation  of  the  scale 
on  karst  permeability  or  hydraulic  con- 
ductivity (After  Basmaci  et  al. ,  1977) . 
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A  spring  is  any  natural  discharge  of  groundwater  on 
the  surface  that  is  large  enough  to  flow  into  a  rivulet 
(Davis  and  Dewiest,  1966)  .  Seeps  discharge  in  the  form 
of  surface  flow,  and  as  a  result  do  not  form  rivulets. 

Spring  discharge  is  a  function  of  the  aquifer  per- 
meability, area  contributing  recharge  to  the  aquifer,  and 
the  quantity  of  recharge.  A  useful  analytical  tool  in 
aquifer  assessment,  related  to  spring  discharge,  is  pre- 
sented by  White  et  al.  (1974).  The  authors  assert  that 
spring  hydrographs  should  provide  information  on  response 
of  aquifers  to  precipitation  events.  The  spring 
hydrograph,  Figure  16,  depicts  two  key  parameters  than 
can  be  obtained.  The  first  parameter,  the  discharge 
ratio,  is  defined  as  the  ratio  of  maximum  discharge,  Q  , 
to  base  flow  discharge,  Q_ , 

Discharge  Ratio  =  QM/Q„  (2.28) 

The  discharge  ratio  is  a  measure  of  the  "flashiness"  of 
an  aquifer.  Another  parameter  that  can  be  extracted  from 
a  spring  hydrograph  is  the  characteristic  response  time, 
t  ,  which  is  defined  as  the  length  of  time  necessary  for 
the  discharge  to  fall  to  1/e  (approximately  -37%)  of  Q  . 
The  exhaustion  coefficient,  e,  is  also  related  to 
Burdon's  (1963)  exhaustion  coefficient  by 


t   =  -  =  exp(l)  =  0.368  (2.29) 

e 

The  recession  limb  of  the  hydrograph  provides  the  follow- 
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Figure  16.   Schematic  hydrograph  for  high  discharge 
from  a  carbonate  spring  (After  White 
et  al.  ,  1974) . 
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ing  exponential  decay  function: 

Q  =  Q  exp(-et)  (2.30) 

where  Q  is  discharged  at  time  t,  and  Q  is  discharge  at 
time  equal  to  zero.  White  (1974)  maintains  that  dis- 
charges near  the  hydrograph  crest  don't  fit  the  simple 
exponential  form  very  well. 

White  et  al.  (1974),  identified  at  least  four 
hydrogeologic  factors  that  control  spring  response:  (1) 
aquifer  hydrogeology ,  (2)  stage  of  development  of  the 
conduit  system,  (3)  catchment  area  for  the  springs,  and 
(4)  piezometric  level  of  the  groundwater.  Table  4  lists 
typical  values  of  discharge  ratio  and  response  time  for 
fast,  intermediate,  and  slow  response  springs.  In 
addition,  the  rate  of  recovery,  t_ ,  for  springs  is  a 
function  of  the  remaining  storage  of  groundwater  (Figure 
16)  . 

The  variability,  V*  ,  of  a  spring  can  be  determined 

cl 

by  the  following  formula: 

Q     -  Q  • 
V   =  -^ SIS  (2b31) 

Qmd 

where  Q     is  the  maximum  discharge,  Q  .   is  the  minimum 

max  3     min 

discharge,  and  Q  ,  is  the  median  discharge  (Davis  and 
Dewiest,  1966)  .  The  time  period  for  calculation  of 
variability  must  be  specified.  Variability  is  a  key  in- 
dicator of  the  storage  capacity  of  an  aquifer.  For  exam- 
ple, aquifers  of  large  storage  capacity  may  have  only  a 
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Table  4.   Calculated  Discharge  Parameters 
for  Selected  Karst  Springs 


Springs 


max  base 


R 
(days) 


Source  of  Data 


Fast-Response  Springs 

Rock  Spring,  Pa.  4  2 

Penns  Cave,  Pa.  96 

Davis  Spring,  W.Va.  91 


7.3     Jacobson  (1973)  , 

19.0    White  and 

mack  (1968) 

4.1     USGS 


Intermediate-Response  Springs 

Thompson  Spring,  PA     9.5  70 

Tuscumbia  Spring,  AL     7.0  65 

Aghia  Eleousa  (Greece)   7.5  86 


Jacobson  (1973 

USGS 

Aronis  et  al . , 
(1961) 


Slow-Response  Springs 


San  Marcos  Springs,  TX   1.6 
Silver  Spring,  FL       1.5 


Ras-el-Ain  Spring, 
Syria 


1.16 


2,070 
(5.67 
years) 


Garza  (1962) 

Ferguson  et  al., 
(1947) 


Burdon  and 

Safadi  (1963) 


(after  White,  1974) 
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slight  fluctuation  in  discharge,  and  as  a  result,  a  low 
variability  (less  than  one) . 

Almost  all  first  magnitude  springs  eminate  from 
basalt,  carbonate,  boulder  or  gravel  aquifers  (Table  5). 
Other  aquifers  such  as  sandstones,  conglomerates,  and 
alluvial  sand  lack  sufficient  transmissivity  to  form 
first-  or  second-magnitude  springs. 

The  most  common  cause  of  spring  occurrence  is  later- 
al or  vertical  variation  of  permeability.  Vertical  vari- 
ation would  include  a  permeable  unit  overlying  an  imperm- 
eable one.  Lateral  variation  could  be  structurally  con- 
trolled by  faulting  where  a  permeable  unit  is  offset  next 
to  an  impermeable  one.  Other  examples  of  spring  occur- 
rence are:  (1)  intensely  fractured  zones  of  brittle  rock 
caused  by  faulting;  (2)  granular  or  unconsolidated  rocks 
crossed  by  faults  that  become  an  impediment  to 
groundwater  flow;  and  (3)  folded  and  tilted  geologic 
units,  permeable  or  impermeable  or  both,  exposed  at  the 
surface  (Davis  and  Dewiest,  1966)  . 

The  remainder  of  this  section  will  discuss  method- 
ology to  properly  assess  spring  flow  parameters.  The 
primary  emphasis  is  on  carbonate  rocks,  since  this  par- 
ticular lithology  is  the  dominant  rock  type  in  the  study 
area. 

To  determine  the  recharge  area  for  a  spring  requires 
an  extensive  amount  of  data  collection.   Preliminary  in- 
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Table  5.   Meinzer's  Classification 
of  Spring  Discharge 


Magnitude 


English  Units 


Metric  Units 


First        Greater  than  100  ft  /sec 

3 
Second       10  to  100  ft  /sec 

3 
Third        1  to  10  ft  /sec 

3 

Fourth       100  gal/min  to  1  ft  /sec 

Fifth  10  to  100  gal/min 

Sixth  1  to  10  gal/min 

Seventh  1  pt/min  to  1  gal/min 

Eighth  less  than  1  pt/min 


Greater  than  2.83  m  /sec 

3 
0.283  to  2.83  m  /sec 

28.3   to   283   liters/sec 

6.31   to   28.3    liters/sec 

0.631   to  6.31   liters/sec 

63.1   to  631  ml/sec 

7.9   to   63.1   ml/sec 

less  than   7.9  ml/sec 


(after   Davis    and  Dewiest,    1966) 
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formation  gathering  includes  acquisition  of:  black  and 
white  aerial  photographs  of  the  study  area,  geologic  maps 
and  hydrogeologic  maps  (if  available) ,  and  previous 
hydrogeologic  findings.  Information  obtained  from  aerial 
photographs  include:  geomorphology ;  lineaments  traces, 
such  as  fractures  and  faults;  formation  contacts;  and 
attitudes  of  formations  (direction  of  dip  and  strike) . 

The  next  phase  of  the  investigation  entails  on-site 
reconnaissance  that  substantiates  or  refutes  previous 
geologic  investigations.  Fracture  measurements  are  taken 
at  certain  locations  to  determine  orientation  and  inten- 
sity. In  addition,  observation  of  soil  cover,  vegetation 
cover,  and  landforms  are  recorded. 

Fracture  measurements  could  conceivably  lend  insight 
to  flow  direction  and  deformational  behavior  of  the  rock. 
If  a  strong  fracture  orientation  exists  in  a  given 
locale,  combined  with  known  potentiometric  levels,  direc- 
tion of  flow  can  be  reasonably  ascertained  (Figure  17)  . 
Also,  strong  correlation  exists  between  orientations  of 
lineations  observed  from  aerial  photographs  versus  mea- 
surements of  fractures  in  the  field  (Hanson,  1973). 
Hanson  (1973)  concluded  that  96%  of  the  springs  in  his 
study  area  are  associated  with  mapped  linear  trends. 

Stearns  (1967)  demonstrated  that  the  frequency  of 
fractures  is  a  function  of  lithology  (Figure  18) .  Frac- 
ture frequency  is  approximately  four  times  higher  for 
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Figure  17.   Two-dimensional  fracture  pattern  showing 
the  relation  between  true  flow  and  direc- 
tion inferred  by  drawing  orthogonal  lines 
to  the  potentiometric  contours  (After 
Davis  and  DeWiest,  1966) . 
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Figure  18. 


Histogram  depicting  frequencies  of 
fractures  for  various  lithologies 
(After  Stearns,  1967). 
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dolomites  than  it  is  for  limestones  (Figure  18)  .  The 
fracture  space  index  on  the  vertical  scale  of  the 
histogram,  Figure  18,  is  the  average  number  of  fracture 
planes  normal  to  a  traverse  of  100  feet.  In  addition, 
Stearns  (1967)  maintains  that  fractures  can  be  used  to 
predict  attitude  of  faults  and  fault  types. 

Morphogenetic  development  of  karst  topography  is  de- 
pendent on  solubility  of  the  carbonate  and  the  amount  of 
rainfall  (Tuttle,  1975) .  Figure  19  depicts  the  four 
stages  of  karst  development  of  a  soluble  carbonate  rock. 
Usually,  dolomites  are  more  resistant  to  dissolution 
activities  than  are  limestones  (Jakucs,  1977).  Jakucs 
(1977) ,  asserts  that  dolomites  very  rarely  duplicate 
textbook  karst  forms,  and  in  actuality  are  semi-karstic 
only.  Another  interesting  item  put  forth  by  Jakucs 
(1977)  ,  open  fissure  (fracture)  density  and  permeability 
of  Lower  Paleozoic  limestone  rocks  are  quite  low,  even 
though  aggregate  fracture  density  is  high  (Figure  20) . 

According  to  Dreiss  (1974)  ,  the  extent  of  conduit 
development  in  dolomites  is  a  function  of  increasing 
amounts  of  calcium  carbonate,  opaque  minerals,  and  poros- 
ity. In  addition,  Dreiss  (1974)  asserts  that  fine- 
grained carbonates  are  more  susceptible  to  karstif ication 
than  coarse-grained  carbonates;  attributing  this  to  in- 
creased surface  area  in  fine  grained  carbonates.  Also, 
she  maintains  that  bedding  planes  rather  than  fractures 
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Figure  19.   Four  stages  of  karstic  landform 
development  (After  Tuttle,  197  5) 
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Figure  20. 


Permeability,  open  fissure  density, 
and  aggregate  fissure  density  of 
limestone  rocks  versus  the  geologic 
time  scale  (After  Jakucs ,  1977). 
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in  carbonates,  control  cavern  occurrence. 

It  is  important  during  the  field  phase  that  an  in- 
vestigator take  note  of  geomorphological  development  of 
carbonates  in  the  outcrop  area.  Davis  and  Dewiest 
(1966)  state  that  a  carbonate  unit  is  permeable  at  depths 
greater  than  1,000  feet  if  solution  openings  or  swallow 
shafts  and  closely  spaced  joints  are  present,  and  shale 
partings  are  absent.  Increased  permeability  in  carbon- 
ates occurs  along  the  crests  of  anticlines  and  the 
troughs  of  synclines  (Davis  and  DeWiest,  1966) .  In  addi- 
tion, Davis  (1966)  maintains  that  in  relatively  homoge- 
neous, carbonate  rocks,  the  greatest  amount  of  solution 
activity  occurs  at  natural  discharge  points. 

The  field  phase  should  include  data  collection  and 
analysis  of  faults  as  a  structural  parameter  that  governs 
groundwater  flow.  Kastning  (1977)  asserts  that  no  gener- 
al statement  of  fault  influence  on  groundwater  flow  in 
carbonates  cannot  be  made  without  some  detailed  examina- 
tion. 

Tracers  such  as  tritium,  and  helium  can  give  much 
needed  insight  to  the  groundwater  flow  system  (Davis  and 
DeWiest,  1966)  .  Tracers  can  be  used  to  delineate  the 
recharge  area,  determine  seepage  velocity  (Marine,  1967), 
and  determine  local  or  regional  flow  systems. 

Local  or  regional  flow  systems  (Figure  21)  can  be 
deduced  from  temperature  logs,  variation  of  spring  dis- 
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charge,  and  the  degree  of  mineralization  (Fetter,  1980)  . 
Temperature  gradients  within  the  earth's  upper  crust  in- 
creases at  about  1°C  per  100  meters  of  depth.  Regional 
systems  are  generally  characterized  by  a  higher  degree  of 
mineralization.  In  addition,  water  temperatures  of  local 
flow  systems  are  close  to  the  mean  annual  air  tempera- 
ture. 

The  use  of  black  and  white  aerial  photographs,  a 
form  of  remote  sensing  analysis,  constitutes  a  viable 
preliminary  reconnoiter  of  a  project  area.  Radar  imagery 
is  another  useful  remote  sensing  technique  for  geological 
investigations.  SLAR  (Side  Looking  Airborne  Radar),  is 
extremely  useful  in  detecting  lineations  in  heavily  for- 
ested areas  (Sabins,  1978). 

Observation  well  hydrographs  combined  with  spring 
hydrographs  provide  extremely  useful  information  on 
aquifer  properties  (Davis  and  Dewiest,  1966).  Observa- 
tion wells,  located  adjacent  to  springs,  provide  informa- 
tion on  piezometric  level  versus  the  amount  and  periods 
of  spring  discharge.  In  addition,  well  hydrographs  de- 
pict long-term  deprivation  effects  due  to  decreasing 
amounts  of  rainfall,  and  nearby  production  wells  (Figure 
22)  . 

Fractures  and  solution  openings,  formation  contacts, 
depth  of  well,  and  extent  and  condition  of  casing  in 
wells;  can  be  ascertained  by  the  use  of  television  logs. 
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Other  common  borehole  geophysical  techniques  are:  caliper 
logs,  resistivity  logs,  temperature  logs,  spontaneous  po- 
tential, and  nuclear  logs  (Fetter,  1980) .  Information 
obtained  directly  from  geophysical  well  logging  includes 
porosity,  lithology,  temperature,  and  water  zones. 

Surficial  geophysical  methods  include  seismic, 
resistivity,  and  gravimetric  surveys.  Information  ob- 
tained from  surveys  include:  vertical  variation  in 
lithology,  structural  disposition  of  the  area,  and  lo- 
cation of  faults. 

Petrophysical  investigations  of  rock  samples 
acquired  from  field  recce  entail  permeability  and  porosi- 
ty measurements.  Falling-head  and  constant-head  permea- 
meter  tests  are  conducted  on  unconsolidated  or  loosely 
cemented  sediments  (Davis  and  Dewiest,  1966)  .  High- 
pressure  gas  permeameter  tests  can  measure  sample  per- 
meability of  indurated  sediments.  Klinkenberg  (1941)  ac- 
counted for  slip-flow  phenomena  associated  with  gas  flow 
through  porous  media  by  the  following  equation: 

K,  =  K   (1  +  b/p)  (2.31) 

where  K.  is  the  intrinsic  permeability  of  the  media,  Ko 
is  the  intrinsic  permeability  as  observed  for  an  incom- 
pressible fluid,  p  is  the  average  flowing  pressure,  and  b 
is  a  constant  characteristic  of  both  gas  and  the  medium. 

Typical  methods  to  determine  porosity  of  rock  sam- 
ples in  the  laboratory  include:   the  direct  method,  the 
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gas  expansion  method,  the  mercury  injection  method,  the 
imbibition  method,  and  the  statistical  method  (Collins, 
1961) .  The  most  accurate  method  to  determine  porosity  is 
the  imbibition  technique.  The  gas  expansion  method  based 
on  Boyles  gas  law  can  determine  pore  volume  by  the  fol- 
lowing equation: 

Volume  of  pores  =  V   =  V_  -  V   -  V,  — - —        (2.32) 

P    B    a    b  P  _P 

where   V   =   bulk  volume  of  the  sample, 
B 

V   =   volume  of  sample  chamber, 
a 

V.   =   volume  of  evacuated  chamber, 

D 

P-   =   initial  pressure, 
P„   =   final  pressure  (Collins,  1961). 
This  method  is  not  quite  as  accurate  as  other  techniques, 
yet  it  is  the  most  widely  used  method  in  industry. 

Houstin  et  al.  (1982),  assert  that  the  most  effect- 
tive  means  to  assess  hydraulic  conductivity  of  an  aquifer 
is  to  conduct  field  pump  tests  of  long  duration  with  at 
least  two  observation  wells,  in  conjunction  with  labora- 
tory tests.  Also,  Houstin  et  al.  (1982),  claim  the  best 
estimate  of  insitu  hydraulic  conductivity  is  given  by  the 
following  equation: 

K  ■  (KOUT  •  KIN>1/2  (2-33) 

where  K  is  hydraulic  conductivity  measured  from  draw- 
down tests  (pump  test) ,  and  K  is  hydraulic  conductivity 
measured  from  infiltration  tests  (rate  of  rise  tests) . 
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Extensive  literature  exists  on  well  hydraulics  for 
granular  media.  However,  Gringarten  (1982)  states  that 
there  is  scant  field  evidence  in  the  literature  that  val- 
idates the  double-porosity  approach  in  fractured  media. 
Yet,  Briz-Kishore  et  al.  (1982),  conducted  extensive  pump 
tests  in  fracture  crystalline  rock  in  India;  and  conclud- 
ed that  the  majority  of  the  drawdown  occurred  within  10 
minutes,  and  that  the  drawdown  curves  behaved  very  simi- 
lar to  that  of  porous  media.  In  addition,  Briz-Kishore 
et  al.  (1982),  surmised  that  the  amount  of  discharge  for 
a  given  well  has  no  relation  with  the  number  of  fractures 
but  depends  only  on  the  productivity  of  any  single  frac- 
ture . 

Gringarten  (1982)  conducted  pump  and  build-up  tests 
in  a  limestone  karst  aquifer  in  1976  (Figure  23) .  Obser- 
vation wells,  P.,  P_,  and  P  were  located  at  approximate- 
ly the  same  distance  from  the  pumping  well;  yet  all  three 
had  distinctively  different  responses.  Gringarten  (1982) 
knew  that  observation  wells  P0  and  P_  were  connected  to 
the  pumping  well  by  a  horizontal  fracture.  With  the  use 
of  the  packer  test  and  the  horizontal  model  theory 
(Gringarten  et  al.,  1973)  for  single  fracture  flow,  he 
determined  the  average  hydraulic  conductivity  to  be  .28 
ft/day. 

Less  costly  means  to  determine  hydraulic  conductiv- 
ity include  slug  tests  (Bouwer,  1978)  .   In  essence,  a 
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slug  test  entails  instantaneous  removal  of  water  with  a 
bucket  or  bailer,  and  measuring  the  corresponding  rise  of 
water  with  a  pressure  transducer  versus  time.  Another 
technique  commonly  used  is  a  closed  cylinder  filled  with 
sand  or  other  material,  submerged  in  the  well  until  water 
equilibrium  is  attained,  and  then  quickly  pulling  it  out, 
and  measuring  the  rate  of  rise.  Slug  tests  are  useful 
for  partially  penetrating  wells  in  confined  or  unconfined 
conditions . 

Potentiometric  surveys  of  the  aquifer  enable  the  in- 
vestigator to  quantify  the  recharge  area  responsible  for 
spring  flow.  This  survey  is  normally  conducted  after  the 
investigation  of  the  structural  geology  and  lithology  of 
the  project  area  has  been  determined.  Steel  tape  and 
chalk  are  commonly  used  to  measure  water  depth  from  a 
known  surface  elevation  (Figure  24)  .  Then,  piezometric 
elevations  can  be  determined  at  each  point  of  interest  in 
the  project  area  and  plotted  on  an  overlay.  Finally,  a 
potentiometric  map  can  be  constructed  from  the  plotted 
piezometric  elevations,  depicting  distribution  and  occur- 
rence of  equipotential  lines.  General  flow  directions 
and  water  divides  can  be  reasonably  determined  from  the 
potentiometric  map  (Davis  and  Dewiest,  1966)  . 

Water  quality  correlation  can  be  used  to  identify 
aquifers  responsible  for  spring  flow  (Walton,  1970)  .  The 
Stiff  (1951)  graphical  method  using  four  parallel  axes 
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Figure  23.  Drawdown  and  build-up  curves  for  ob- 
servation wells  in  a  limestone  karst 
aquifer  (After  Gringarten,  19  82) . 


Reel 


Graduated  tape  with  chalk 
rubbed  on  lower  part 


Length  of  wetted  tape 
Lead  weight 


Figure  24. 


A  sketch  depicting  water  table  depth 
measurement  utilizing  a  steel  surveyor's 
tape  that  can  be  read  with  an  accuracy 
of  .005  feet  (After  Davis  and  Dewiest, 
1966)  . 
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and  one  vertical  axis,  may  be  useful  in  making  compari- 
sons of  water  types  (Figure  25) .  Four  cations  (Na  +  K, 
Ca,  Mg  and  Fe)  are  plotted  along  each  axis  to  the  left  of 
zero,  and  four  anions  (CI,  HCC>  ,  SO,,  and  CO..)  are  plot- 
ted to  the  right  of  zero  (Figure  25)  .  Connecting  the 
points  of  plotted  anions  and  cations,  depicts  a  charac- 
teristic shape  for  water  of  a  given  aquifer. 

Storage  coefficient  of  an  aquifer  can  be  obtained 
from  well  hydraulic  methods.  Another  method  to  determine 
storage  coefficient  is  based  on  water  level  fluctuations 
of  a  confined  aquifer  as  a  result  of  atmospheric  pressure 
changes.  The  barometric  efficiency,  BE,  is  given  by  the 
following  equation  (Todd,  1959)  : 

Ah 
BE  =  =  l/(a/nB  +  l)  (2.34) 


Ap  /y 
*a  'w 


where 


Y    =   specific  weight  of  water, 

Ah   =   change  in  water  level, 

Ap   =   change  in  barometric  pressure, 

a    =   compressibility  factor  of  an  aquifer, 

n   =   porosity. 

Also,  the  storage  coefficient  of  a  confined  aquifer  is 

defined  by  the  following  equation: 

S   =  Y  b  (a  +  nB  )  (2.35) 

c     w 

where  b  is  the  aquifer  thickness,  and  $  is  the  compress- 
ibility factor  of  water.   Once  a  barometric  efficiency  is 
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obtained  (values  range  from  .25  to  .75),  the  aquifer 
compressibility,  a,  can  be  determined.  Then,  inserting 
the  known  value  of  a  into  Equation  (2.35),  yields  the 
sought  after  storage  coefficient.  Davis  and  Dewiest 
(1966)  assert  that  barometric  efficiency  is  directly  pro- 
portional to  aquifer  rigidity;  where  Ah  (Equation  2.33) 
is  small  for  rigid  aquifers  of  low  permeability,  and 
large  for  unconsolidated  sediments  of  high  permeability. 

Deprivation  of  spring  discharge  can  be  attributed  to 
decreasing  amounts  of  rainfall  and  nearby  production 
wells.  Harp  et  al .  (1975)  ,  devised  a  methodology  in 
which  percentage  of  deprivation  can  be  determined  for 
each  production  well  relative  to  one  well  in  the  same 
aquifer.  The  authors  used  a  well  hydraulics  equation  for 
a  fully  penetrating  well  in  an  unconfined,  homogeneous 
aquifer  of  infinite  aerial  extent,  to  obtain  total 
drawdown  s(x,y),  due  to  adjacent  production  wells  (x., 
Y±>: 

2          X    n  2  2 

s  (x,y)  =  Z      Q   ln[(x  -  x . ) *  +  (y  -  y.)z]    (2.36) 

2ttK  i=l   1  1  x 

where 

s(x,y)   =   total  drawdown,  in  feet,  at  some  point 

x,y,  due  to  n  production  wells, 


Q.(x.,y.)   =  discharge  in  gallons  per  minute  for  the 
i  '  production  well  at  some  point  x.,y.. 
The  authors  reasoned  that  s  is  a  function  of  discharge  Q, 
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2 
and  the  natural  logarithm  of  r 

s2  =  f  (Q,  In  r2)  (2.37a) 

They  deduced  that  discharge  Q  is  inversely  proportional 

2 
to  the  natural  logorithm  of  r  .   Based  on  the  preceding 

logic,  Harp  et  al.   (1975)  ,  assert  that  the  deprivation 

contribution  for  each  production  well  in  relation  to  a 

well  in  a  common  aquifer  can  be  expressed  as  follows: 

^       Q±/ln  (r2) 
Deprivation  contribution  of  i   well  = (2.37b) 

n         2 
I   (Q./ln(rf)) 

i=l   1     1 
The  authors  concluded  that  this  concept  should  be  appli- 
cable to  confined  aquifers  also. 

In  summary,  springs  represent  natural  discharge 
points  for  a  carbonate  aquifer.  Topographic  relief  pro- 
vides the  hydraulic  gradient,  and  geologic  structure  in- 
fluences the  direction  of  groundwater  to  discharge  points 
(White,  1977).  The  magnitude  and  variability  of  spring 
flow  is  closely  related  to:  (1)  the  hydrogeological 
properties  of  the  aquifer,  (2)  the  piezometric  level,  (3) 
the  area  of  the  catchment  basin,  (3)  the  amount  of  annual 
recharge,  and  (4)  deprivation  effects  of  adjacent  produc- 
tion wells. 

2.5   Previous  Hydrogeologic  Findings  of  Artesian  Spring 
Flow,  Sulphur,  Oklahoma,  Area 

Charles  N.  Gould  (1906),  a  geologist,  forwarded  a 

report  to  the  Secretary  of  the  Interior.   The  contents  of 
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the  report  focused  primarily  on  location,  amount  of  flow, 
status  of  development,  and  character  of  the  springs  in 
Piatt  National  Park  (Travertine  District,  Chickasaw  Na- 
tional Recreation  Area) . 

The  report  revealed  that  a  total  of  33  springs  were 
located;  of  which  6  were  freshwater  springs,  and  the  re- 
mainder were  mineralized.  Gould  (1906)  reported  lo- 
cations of  the  springs  in  relation  to  well  known  land- 
marks or  springs.  Unfortunately,  a  grid  system  was  not 
utilized,  and  the  exact  location  of  some  of  these  springs 
that  have  sinced  ceased  flowing,  is  unknown  to  date. 

Gould  (1906)  estimated  freshwater  discharge  to  be 
approximately  2,000  gallons  and  1,500  gallons  per  minute 
for  Antelope  and  Buffalo,  respectively.  At  the  time  of 
the  report  both  springs  were  not  developed  or  modified. 

Gould  (1906)  postulated  that  the  source  of  water  for 
the  springs  in  the  park  was  from  the  upper  part  of  the 
Simpson.  He  based  this  on  one  log  from  a  well.  He  did 
not  state  the  location  of  this  well.  Also,  he  reported 
that  4  artesian  wells,  north  of  the  park,  in  the  city  of 
Sulphur,  were  flowing  at  an  approximate  rate  of  200,000 
gallons  per  day.  He  did  not  divulge  locations  of  the 
artesian  wells. 

Gould  and  Schoff  (1939)  conducted  a  hydrogeologic 
investigation  of  the  park  during  the  summer  of  1939  when 
Buffalo  and  Antelope  ceased  flowing.   The  authors  con- 
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eluded  that  the  dwindling  spring  flow  in  the  park  was  due 
to  diminished  amounts  of  rainfall,  and  aquifer  depriva- 
tion of  the  free  flowing  wells  in  the  city  of  Sulphur. 

Gould  and  Schoff  (1939)  reported  that  14  springs 
could  not  be  located  since  the  initial  investigation  in 
1906.  Also  noted,  was  that  the  Black  Sulphur,  Pavillion, 
and  Hillside  Springs  combined  flow  for  1939  was  one-fifth 
of  the  reported  value  in  1906. 

Schoff  (1939)  determined  the  annual  recharge  to  be 
approximately  1  foot.  The  source  or  area  of  recharge  for 
the  mineralized  springs  and  the  Sulphur  artesian  basin 
was  ascertained  to  be  approximately  12  square  miles  of 
the  outcropping  Simpson  Group,  in  the  Sulphur  Syncline, 
southeast  of  Sulphur  (Plate  I)  .  Schoff  (1939)  surmised 
that  the  surface  elevation  of  the  outcrop  area  was 
approximately  200  ft.  higher  than  the  surface  elevations 
of  the  mineralized  springs  and  flowing  wells  in  the  city 
of  Sulphur,  which  would  account  for  the  artesian  condi- 
tion. In  addition,  the  authors  postulated  that  the 
permeable  sandstone  beds  provided  the  main  source  of 
water  for  the  mineralized  artesian  springs.  The  authors 
reasoned  that  the  pressurized  water  of  the  Simpson  Group 
flowed  through  the  joints  and  fractures  in  the  overlying 
Viola  and  Pontotoc. 

Schoff  (1939)  claimed  that  the  source  of  water  for 
Antelope  and  Euffalo  Springs,  and  the  city  of  Sulphur 
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wells  was  from  the  outcrop  of  the  Pontotoc,  east  and 
south  of  Sulphur.  He  also  speculated  that  because  the 
catchment  basin  for  the  freshwater  springs  is  small, 
spring  discharge  is  particularly  dependent  on  sustained 
precipitation  events.  In  addition,  Schoff  (1939)  stated 
that  freshwater  for  the  springs  and  the  city  of  Sulphur 
wells  could  originate  from  the  Arbuckle  Group.  Also,  the 
investigators  assert  that  the  Pontotoc  and  the  Simpson 
are  to  a  certain  extent  independent  aquifers. 

Prior  to  and  during  the  investigation,  an  inventory 
of  the  artesian  wells  in  the  city  of  Sulphur  were  under- 
taken by  the  city  officials.  A  list  of  12  wells  were 
submitted  to  the  authors  with  a  total  combined  discharge 
of  17,500  gallons  per  minute  for  the  calendar  year  1939. 
Finally,  the  authors  concluded  that  the  Vendome  Well  (an 
artesian  well)  was  the  primary  reason  for  cessation  of 
spring  flow  or  reduced  spring  discharge.  Gould  and 
Schoff  (1939)  attributed  this  deprivation  to  a  high,  sus- 
tained flow  rate  of  2500  gallons  per  minute  since  1922, 
and  its  close  proximity  to  the  park. 

Hart  (1972)  deduced  from  his  investigations  of  the 
area,  that  the  west-northwest  trending  Sulphur  Fault  is 
the  major  controlling  factor  in  the  quality  and  quantity 
of  water  that  flows  from  the  springs  in  the  park,  and 
wells  in  the  city  of  Sulphur.  He  postulated  that  the 
Sulphur  Fault  passes  somewhere  near  Antelope  and  Buffalo, 
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exact  location  was  unknown. 

Contrary  to  previous  investigators,  Schoff  and  Gould 

(1939)  ;  Hart   (1972)   surmised  that  the  outcrop  of  the 

Arbuckle  Group,  east  of  Sulphur,  is  the  source  of  water 

for  Antelope  and  Buffalo  Springs.   He  substantiated  this 

finding  by  the  following  correlations: 

(1)  water  quality  analysis  (dissolved-solids  con- 
tent of  300  mg/1,  calcium  carbonate  type), 

(2)  potentiometric  head  in  the  outcrop  area  of  the 
Arbuckle  Group  versus  surface  elevations  of 
Buffalo  and  Antelope,  and 

(3)  well  hydrographs  of  the  outcrop  area  of  the 
Arbuckle  Group  versus  spring  hydrographs.   He 
also  determined  that  the  Pontotoc  is  relatively 
impremeable,  and  serves  as  a  confining  layer, 
lying  unconformably  over  older  formations. 

Hart  (1972)  claimed  that  another  recharge  area  was 
responsible  for  spring  flow  in  the  park.  The  author  as- 
serts that  south  and  southeast  of  Sulphur,  the  outcrop 
area  of  the  Simpson  Group  is  the  recharge  area  for  the 
only  "true"  mineral  springs  in  the  park,  Bromide  and  Med- 
icine. Location  of  these  two  springs  is  in  the  far  west 
end  of  the  park  (Plate  III)  .  Discharge  from  these  two 
mineral  springs  was  estimated  to  be  1  gallon  per  minute. 
Hart  (1972)  attributed  the  steady,  low-flow  rate  and 
mineralization  (dissolved-solids  content  of  4,000  mg/1, 
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sulphur  type)  to  the  slow,  but  steady  recharge  rate  of 
the  outcrop  area  of  the  Simpson  Group,  and  the  overlying 
confining  layer  of  the  Pontotoc. 

Hart  (1972)  reported  that  mixed  waters,  consisting 
of  water  from  the  Simpson  and  the  elevated  Arbuckle  were 
responsible  for  the  character  of  the  water  of  the  springs 
(Hillside,  Pavillion  and  Black  Sulphur)  in  the  middle  of 
the  park,  and  north  of  the  park  in  the  flowing  wells  of 
the  city  of  Sulphur.  The  results  of  the  water  quality 
analysis  in  the  form  of  Stiff  diagrams,  graphically  de- 
picted higher  than  normal  amounts  of  sodium  and  chloride 
than  are  normally  present  in  waters  of  the  Arbuckle  Group 
(Figure  26)  .  The  dissolved  solids  content  of  these 
springs  was  reported  to  be  from  530  to  816  mg/1,  with  an 
estimated  combined  flow  at  100  gallon  per  minute.  The 
dissolved  solids  content  of  the  Vendome  Well  was  reported 
to  be  1,200  mg/1,  with  an  estimated  discharge  rate  of 
1,000  to  2,000  gallons  per  minute. 

Hart  (1972)  concluded  that  the  Sulphur  Fault  is  the 
primary  controlling  feature  isolating  the  freshwater 
springs,  Antelope  and  Buffalo,  from  the  remainder  of  the 
springs  in  the  park.  He  attributes  the  occurrence  of 
mixed  waters  in  the  center  of  the  park  and  north  of  the 
park  to  waters  from  the  small  veneer  of  Simpson  sediments 
overlying  the  elevated  Arbuckle  dolomite.  He  determined 
that  the  springs,  Bromide  and  Medicine,  originate  from 
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the  sandstone  of  the  Simpson  Group.  He  substantiated 
this  particular  finding  by  conducting  a  water  quality- 
analysis  of  the  Simpson  Group  at  depths  of  2,525  feet  and 
1,575  feet  in  T.IN.,  R3E.,  Sec.  7  (Figure  26).  Finally, 
Hart  (1972)  concurred  with  previous  investigators,  the 
occurrence  of  artesian  springs  in  the  park  was  the  result 
of  pressurized  water  flowing  upward  through  fractures  and 
joints  in  the  younger,  overlying  formations;  ultimately 
discharging  from  the  Pontotoc  Group. 

Harp  et  al.  (1976),  conducted  an  investigation 
primarily  to  determine  the  deprivation  effects  of 
production  wells  on  spring  flow  in  the  park.  Also,  the 
authors  presented  additional  information  to  clarify  and 
elucidate  upon  Hart's  (1972)  findings. 

Harp  et  al.  (1976),  ascertained  that  two  main 
aquifers,  namely  the  Simpson  and  the  Arbuckle,  were  re- 
sponsible for  spring  discharge  in  the  park.  They  es- 
timated the  recharge  area  for  the  mineralized  springs  to 
be  approximately  20  square  miles  of  the  Simpson  Group 
outcrop  located  south  and  southeast  of  the  park.  The  re- 
charge area  for  Buffalo  and  Antelope  was  estimated  to  be 
40  square  miles  of  the  Arbuckle  Group  outcrop,  east  of 
the  park.  In  addition,  the  authors  correlated  spring 
discharge,  precipitation  events,  water  quality  analysis, 
and  observation  well  hydrographs  (two  observation  wells 
were  completed  in  1972)  to  substantiate  the  source  and 
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occurrence  of  the  freshwater  springs.  The  two  observa- 
tion wells  were  completed  in  the  Arbuckle  Group,  verified 
by  a  well-site  geologist. 

The  authors  provided  information  on  location  (Figure 
27)  ,  water  quality  analysis  (Table  6)  ,  and  discharge 
rates  of  selected  wells  and  springs  in  the  project  area. 
Bromide  and  Medicine  springs  were  not  flowing  at  the  time 
the  authors  conducted  the  investigation.  Harp  et  al., 
estimated  the  combined  discharge  of  Pavillion,  Hillside 
and  Black  Sulphur  Springs  to  be  between  100  and  200 
gallons  per  minute.  The  authors  concurred  with  Hart 
(1972)  ,  that  these  springs  and  the  wells  to  the  north, 
derive  their  water  from  both  the  Simpson  and  Arbuckle 
aquifers . 

Harp  et  al.  (1976),  selected  3  well  fields  and  2 
wells  that  discharge  significant  amounts  of  water 
annually,  that  could  adversely  affect  spring  discharge  in 
the  park  (Figure  28) .  The  authors  presented  pumpage  data 
for  the  three  well  fields,  and  estimated  average  dis- 
charge rates  for  the  two  flowing  wells,  Vendome  and 
Wyandotte . 

To  determine  deprivation,  the  authors  used  the  fol- 
lowing steady-state,  well  hydraulics  expression  for  total 
drawdown  due  to  production  wells  in  an  isotropic,  con- 
fined, homogeneous  aquifer  at  some  point  of  interest: 
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n   Q.       R. 
s  -    Z   -i-  In  (-i)  (2.38) 

i=l  2ttT      r. 

1 

where 

n  =  'n'  number  of  wells, 

s  =  total  drawdown  at  some  point  due  to  'n'  number 

of  wells, 

Q.  =  discharge  rate  at  the  i   well, 

T  =  transmissivity ,  constant  for  a  homogeneous, 

isotropic,  confined  aquifer, 

R.  =  radius  of  influence, 
l 

r.  =  radius  from  the  point  of  interest  to  the  i 

well. 

They  reasoned  that  2irT  is  constant  for  a  homogeneous, 

isotropic  aquifer  of  constant  thickness.   Thus,  equation 

(2.37)  provides  a  simplified  expression  which  relates  the 

effect  of  i    well  to  drawdown  at  the  point  of  interest 

and  assumes  the  form: 

R. 
Qi  In  —  (2.39) 

r . 

l 

Harp  et  al.  (1976),  assert  that  if  r.  >  R.,  then  the  pro- 
duction well  will  exert  no  influence  on  the  point  of  in- 
terest. Hence,  the  proportional  effect  of  each  produc- 
tion well  can  be  evaluated  by  using  the  pumpage  data,  Q, 
the  distance  from  the  i  well  to  the  production  well,  r, 
and  the  radius  of  influence,  R,  for  each  i  well.  The 
authors  state  that  if  R  is  unknown,  then  one  can  assume 
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that  R  is  related  to  Q  by  some  constant  a1  ,  which  is 
characteristic  of  the  aquifer,  such  that 

R  =  a'  Q  (2.40) 

Based  on  the  preceding  logic,  the  investigators  de- 
termined the  percentage  of  deprivation  that  each  produc- 
tion well  exerted  on  the  area  of  Buffalo  and  Antelope 
Springs  (Table  7)  .  The  computated  data  indicated  that 
the  highest  deprivational  effect  in  all  instances,  was 
the  Vendome  Well. 

The  major  recommendation  put  forth  by  Harp  et  al. 
(1976)  ,  was  that  the  flow  of  the  Vendome  Well  and  the 
Wyandotte  Well  be  substantially  decreased  or  plugged. 

Fairchild  et  al.  (1983),  U.S.G.S.  conducted  an  ex- 
tensive hydrologic  investigation  of  the  Arbuckle  Moun- 
tains. Even  though  their  investigations  did  not  focus 
specifically  on  the  Sulphur  area;  their  assessment  of  the 
aquifer  characteristics  of  the  Arbuckle-Simpson  aquifer 
on  a  regional  basis,  is  directly  applicable  to  this 
study. 

Fairchild  et  al.  (1983)  ,  determined  the  average 
storage  coefficient  and  transmissivity  for  the  Arbuckle- 
Simpson  aquifer  in  the  outcrop  area,  to  be   .00  8  and 

2 

15,000  ft  /day,  respectively.   The  authors  assumed  that 

the  aquifer  behaves  in  a  homogeneous  manner  on  a  regional 
basis.  They  determined  the  saturated  thickness  of  the 
aquifer  to  be  approximately  3,500  feet  for  the  outcrop 
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area.  In  addition,  the  authors  estimated  the  recharge  to 
the  aquifer  at  about  4.7  inches  per  year.  Also,  water- 
level  fluctuations  in  the  aquifer  were  reported  to  range 
from  8  to  53  feet  each  year  in  response  to  recharge  from 
rainfall . 

Fairchild  et  al.  (1983),  estimated  that  two-thirds 
of  the  Arbuckle-Simpson  aquifer  consists  of  limestones 
and  dolomites.  The  authors  claim  that  the  occurrence  and 
movement  of  groundwater  in  the  Arbuckle-Simpson  aquifer 
are  strongly  controlled  by  lithology  and  structure.  They 
maintained  that  most  springs  are  near  faults  or  fractures 
that  have  been  enlarged  by  solutioning  of  the  more  dense 
Arbuckle  group.  The  authors  assert  that  the  resultant 
network  of  solution  conduits  occurs  laterally  along  bed- 
ding planes;  and  along  fractures,  vertically  and  hori- 
zontally. Small  sinks  or  swallow  shafts  in  the  valleys 
were  encountered  by  the  investigators  in  the  outcrop 
area.  Also,  Fairchild  et  al.  (1983) ,  studied  caliper 
logs  of  deep  wells  in  the  Arbuckle-Simpson  aquifer.  They 
concluded  that  solution  cavities  exist  to  a  depth  of  at 
least  2,000  feet. 

The  authors  conducted  a  potentiometric  survey 
(Winter  1976,  1977)  of  the  eastern  part  of  the  Arbuckle 
Mountains  to  determine  the  source  and  occurrence  of 
groundwater.  They  concluded  that  generally  the  potentio- 
metric surface  follows  the  topography.   They  also  stated 
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that  the  bulk  of  groundwater  flow  in  the  eastern  portion 
of  the  Arbuckles  is  in  an  eastward  direction  with  a  typi- 
cal hydraulic  gradient  of  20  feet  per  mile. 

Fairchild  et  al.  (1983),  reasoned  that  the  occur- 
rence of  artesian  wells  along  Rock  Creek  in  the  vicinity 
of  Sulphur,  is  due  to  the  erosion  of  the  land  surface  be- 
low the  potentiometric  surface  of  the  confined  Arbuckle- 
Simpson  aquifer.  The  authors  concurred  with  Hart  (1972) 
that  sufficiently  higher  surface  elevations  exist  to  the 
east  in  the  recharge  area,  providing  the  necessary  hy- 
draulic head  to  warrant  artesian  conditions  in  the 
Sulphur  area. 

Fairchild  et  al.  (1983),  conducted  an  analysis  of 
the  ratio  of  recharge  to  average  precipitation  based  on 
regional  techniques  proposed  by  Rorabaugh  (1964,  1966). 
The  authors  reported  that  the  ratio  of  recharge  to  aver- 
age precipitation  for  the  month  of  March  to  be. 32,  and 
for  the  months  of  May  and  June  to  be  .15. 

The  authors  conducted  aquifer  tests  to  determine 
storage  coefficient,  specific  capacity,  and  transmissiv- 
ity  at  point  locations  throughout  the  Arbuckle-Simpson 
aquifer.  The  specific  capacity,  amount  of  constant  dis- 
charge in  gallons  per  minute  divided  by  the  drawdown  in 
feet  for  the  period  of  pumpage,  varied  from  .17  to  4.4 
for  the  Simpson,  and  from  .47  to  104  for  the  Arbuckle. 
The  investigators  reported  that  deep  wells  drilled  in  the 
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Arbuckle  Group  have  the  highest  specific  capacity.  They 
attributed  this  phenomena  to  a  greater  probability  of  en- 
countering interconnected  solution  cavities  and  fractures 
with  increasing  depth.  The  investigators  obtained  infor- 
mation from  land  owners  and  drillers  that  revealed  the 
upper  two  hundred  feet  of  the  Arbuckle  Group  has  a  much 
lower  permeability  than  the  lower  part. 

Estimates  of  transmissivity  were  obtained  from  pump 
tests  and  regional  techniques.  Fairchild  et  al.  (1983), 
used  the  modified  nonequilibrium  equation  by  Ferris 
(1962)  to  determine  transmissivity  T,  from  pump  tests, 

264Q 

T  =   (2.41) 

As 

where  discharge  Q,  is  in  gallons  per  minute,  and  As  is 

change  in  residual  drawdown,  in  feet,  over  one  log  cycle 

of  time.   The  above  equation  assumes  a  fully  penetrating 

well  in  a  homogeneous,  unconfined,  isotropic  aquifer  of 

infinite  areal  extent.   The  investigators  did  realize 

that  the  Arbuckle-Simpson  aquifer  did  not  meet  the  above 

mentioned  criteria,  but  decided  that  such  tests  would 

give  only  a  rough  estimate.   Transmissivity  determined 

2 
from   pump   tests   ranged   from   40   to   2,500   ft  /day. 

Fairchild  et  al.  (1983),  determined  storage  coefficient, 
S  ,  from  pump  tests  by  the  use  of  the  modified  nonequi- 
librium equation  (Ferris,  1962) 
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,3T  t 

S   =  =-^  (2.42) 

c      r2 

where 

t   =  time  intercept,  in  days,  where  the  plotted 
straight  line  intersects  the  zero  drawdown 
axis, 
r  =  distance  from  the  pumping  well  to  the  observa- 
tion well. 

The  values  of  storage  coefficient  were  reported  to  be  in 

-4       -5 
the  range  of  10    to  10 

The  authors  determined  the  average  storage  coeffi- 
cient and  transmissivity  by  regional  techniques  developed 
by  Rorabaugh  (1966),  to  be  .008  and  10,000  ft2/day,  re- 
spectively. 

Another  estimate  of  transmissivity,  T,  was  deter- 
mined by  the  use  of  Darcy's  Law  applied  to  spring  dis- 
charge, 

Q 
T  =  —  (2.43) 

IW 

where 

Q  =  spring  discharge  in  cubic  feet  per  day, 
I  =  hydraulic  gradient,  in  feet  per  feet, 
W  =  width,  in  feet,  of  the  cross  section  through 
which  the  discharge  occurs. 

The  authors  estimated  the  aquifer  transmissivity  to  be 

2 
20,000  ft  /day.   The  transmissivity  of  the  aquifer  was 
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based  on  only  one  spring  in  the  area. 

Fairchild  et  al.  (1983),  devised  a  hydrologic  budget 
for  the  area, 

P'  =  Q»  +  ET  (2.44) 

where 

P'  =  average  annual  precipitation, 

Q'  =  total  average  annual  runoff, 

ET  =  average  annual  evapotranspiration. 
They  estimated  the  annual  average  precipitation  P,  to  be 
38.4  in.  This  value  was  ascertained  from  6  years  of  data 
and  3  rain  gauge  stations.  The  total  average  runoff  Q' , 
was  determined  to  be  7.6  in.  Average  runoff  was  computed 
from  discharge  stations  located  on  Mill  Creek,  Blue  River 
and  Pennington  Creek  at  the  periphery  of  the  outcrop  ar- 
ea. Then,  the  authors  determined  evapotranspiration  to 
be  31  inches  per  year,  or  80  percent  of  the  average  annu- 
al precipitation. 

The  authors  reasoned  that  the  surface  runoff  is 
divided  between  base  flow  and  direct  runoff.  They  ar- 
rived at  an  average  value  of  4.7  in.  per  year  for  base 
flow  or  groundwater  flow.  This  value  was  based  on  low- 
flow  periods  when  surface  runoff  is  nil.  Thus,  ground- 
water flow  accounts  for  12  percent  and  direct  runoff 
accounts  for  8  percent  of  the  average  annual  precipita- 
tion in  the  region. 

Fairchild  et  al.   (1983),  conducted  a  water  quality 
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analysis  of  23  wells,  12  springs  and  12  streams  in  the 
Arbuckle  Mountain  area.  The  authors  reported  an  average 
dissolved  solids  content  of  330  mg/1.  Average  water  tem- 
perature of  all  groundwater  surveyed  was  62°F.  This  av- 
erage temperature  is  almost  identical  to  that  of  the  mean 
annual  air  temperature  of  62.5°.  Also,  samples  taken  at 
depths  to  2,500  feet  showed  little  or  no  change  in  dis- 
solved solids  content. 

The  investigators  pointed  out  that  beds  of  Arbuckle- 
Simpson  aquifer  become  highly  mineralized  when  they  dip 
steeply  beneath  younger  sediments  along  the  periphery  of 
the  outcrop  area.  Fairchild  et  al.  (1983),  cited  work 
done  previously  by  Dott  and  Gunter  (1930)  ,  and  Case 
(1934)  that  dealt  with  the  characteristics  of  waters  of 
the  Ordovician  rocks  that  covered  a  broad  area,  including 
the  Arbuckle  Mountain  area.  The  authors  reported  that 
these  investigators  had  found  that  the  concentration  of 
dissolved  solids  increased  at  a  rate  of  10,000  mg/1  per 
mile  in  the  down-dip  direction. 

Finally,  Fairchild  et  al.  (1983),  tabulated  some 
well  yields  for  the  Simpson  and  Arbuckle  groups.  They 
reported  that  yields  for  the  Simpson  ranges  from  100  to 
400  gallons  per  minute.  Also,  yields  for  the  Arbuckle 
range  from  1  gallon  per  minute  (shallow  wells  in  dense 
dolomite  with  no  interconnected  fractures  or  solution 
cavities)  to  as  high  as  2,500  gallons  per  minute  (deep 
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wells,  in  excess  of  1,000  feet). 

Hensch  (1984)  provided  an  update  to  the  investiga- 
tion conducted  by  Harp  et  al.  (1976).  He  primarily 
addressed  the  stage  versus  flow  relationship  for  the  two 
freshwater  springs,  Buffalo  and  Antelope.  In  addition, 
Hensch  (1984)  was  also  tasked  with  updating  pertinent 
hydrologic  data  for  the  study  area  and  establishing  a 
recession  coefficient  for  the  freshwater  springs. 

Hensch  (1984)  observed  that  Buffalo  Spring  wells  up 
from  a  gravel  bed  into  a  circular  retaining  wall,  25  feet 
in  diameter.  Then,  exits  over  a  4.5  ft.  wide,  rectangu- 
lar-shaped spillway  into  a  creek  bed.  He  also  observed 
additional  water  welling  up  into  the  creek  bed,  immedi- 
ately downstream  from  the  spillway.  According  to  Hensch 
(1984)  flow  in  the  stream  increased  substantially.  Es- 
timating the  discharge  of  springs  in  the  creek  to  be 
three  times  that  of  the  discharge  from  Buffalo  Spring. 

Hensch  (1984)  described  Antelope  Spring  as  a  spring 
eminating  from  the  base  of  the  westward  facing  bluff  into 
a  shallow,  irregular  pool,  approximately  100  square  feet. 
Then,  flow  exits  through  an  opening  in  the  pool  about  5 
feet  wide. 

The  author  stated  that  stage  measurements  at  both 
springs  were  done  by  park  employees  on  an  irregular  ba- 
sis. The  stage  readings,  since  1969,  were  made  on  steel 
pins  embedded  at  the  base  of  each  pool. 
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Hensch   (1984)   ascertained  that  both  springs  exit 

through  a  rectangular  configuration,  similar  to  that  of  a 

broad-crested  weir.   The  author,  together  with  previous 

investigators  (Harp  et  al.  ,  1976)  ,  determined  that  flow 

could  be  roughly  approximated  by  the  following  expression 

Q  =  (C)  (b)  H1'5  (2.45) 

w 

where 

Q  =  flow  in  c.f.s., 

C  =  constant  dependent  on  weir  shape, 

b'  =  width  of  weir, 

H  =  the  depth  of  the  pool  at  a  specified  point. 

Since  C  and  b   are  constants,  the  expressions  can  be  sim- 
w  r 

plified  to  the  following  equation: 

Q  -  K'  H1,5  (2.46) 

where  K'  is  a  constant.  The  investigator  realized  that 
both  springs  do  not  flow  until  a  critical  depth  is 
reached.  Hensch  (1984)  reasoned  that  the  stage  reading 
or  pin  reading  P.R.,  for  each  spring  must  exceed  the  base 
of  the  weir-like  opening  before  flow  begins.  Finally, 
after  field  reconnaissance,  he  formulated  the  following 
flow  equation  for  Buffalo: 

Q  =  10.44  (P.R.  -  1.97)1*5  (2.47) 

and  Antelope; 

Q  =   8.69  (P.R.  -  .6) 1'5  (2.48) 

The  author  attempted  to  determine  the  recession  con- 
stant K  ,  for  the  two  springs.   He  used  a  method  proposed 
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by  Langbein  (1940)  and  is  referenced  by  Davis  and  Dewiest 
(1966).  The  method  entails  plotting  discharge  q  (verti- 
cal scale)  ,  versus  discharge  q  .  ,  where  n  =  t  (incre- 
mental time)  -  1.  Discharge  values  can  be  obtained  from 
the   recession   limb   of   a   stream  hydrograph.    Hensch 

determined  the  K  ,  the  total  flow  recession  constant,  to 

r 

be  approximately  .7  for  both  springs. 

He  mentioned  he  could  not  determine  which  values 

represented  groundwater  flow,  so  he  selected  values,  q 

and  q  ,,,  arbitrarily.   It  is  unclear  as  to  whether  he 
nn+l  J 

understood,  that  all  flow  components  the  spring  hydro- 
graph  represent  groundwater  flow.  The  difference  in  flow 
components  (i.e.,  surface,  interflow,  etc.)  between  the 
stream  hydrograph  and  the  spring  hydrograph  is  termi- 
nology. In  all  likelihood,  he  probably  meant  base  flow, 
in  lieu  of  ground  water  flow.  In  addition,  he  offered  no 
explanation  as  to  why  it  is  necessary  to  determine  a  re- 
cession constant  for  the  two  springs,  other  than  being 
tasked  to  do  so. 

It  is  apparent  from  previous  investigations  in  the 
study  area,  that  the  regularity  of  spring  flow  in  the 
park  is  declining  at  an  alarming  rate.  Two  separate 
aquifers,  the  Simpson  and  the  Arbuckle,  are  responsible 
for  the  waters  of  the  artesian  springs  in  the  park.  The 
geology  of  the  area,  water  quality,  deprivation  effects, 
and  well  and  spring  hydrographs  in  relation  to  precipita- 
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tion  events  have  been  addressed  by  various  investigators 
in  the  past.  It  is  also  understood  that  each  investiga- 
tion provides  additional  input  into  the  analysis  of  the 
artesian  spring  phenomena  in  the  park.  It  is  hoped  that 
this  author  will  provide  additional  insight  to  the 
parameters  that  govern  artesian  spring  flow  in  the  park. 


CHAPTER  III 


METHODOLOGY 


3.1   Geologic  Control  of  Groundwater  Flow 

As  stated  previously,  springs  represent  natural  dis- 
charge points  (sinks)  of  a  carbonate  aquifer.  Determina- 
tion of  source  occurrence  and  movement  of  groundwater  to 
the  springs  requires  an  understanding  of  the  geological 
parameters  that  govern  groundwater  flow  (Table  8) . 

Geologic  maps  of  the  area  were  completed  by  Ham  et 
al.  (1945,  1954),  and  will  be  used  extensively  for  this 
investigation.  Aerial  photographs  (May  and  June,  1963) 
of  the  project  area,  scale  of  approximately  1:20,000, 
were  acquired  from  the  Oklahoma  Geological  Survey  on  a 
loan  basis.  Use  of  the  aerial  photographs  and  a 
stereoscope  will  enhance  preliminary  identification  of 
structures,  formation  contacts,  lineations  (faults  and 
fractures) ,  and  geomorphological  features  of  the  study 
area.  Topographic  maps  (1978),  scale  1:24,000,  will  be 
used  in  conjunction  with  aerial  photographs  to  obtain  lo- 
cations of  findings  to  within  1/16  of  a  mile,  and  provide 
information  on  surface  elevations  to  within  5  feet  of 
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Table  8.   Effect  of  Hydrogeologic  Setting  on 
Carbonate  Aquifers 


Geologic  Element 


Control 


Macro-Structure 
(folds,  faults) 

Topographic  Setting 


Stratigraphic  Sequence 

Mini -Structure 
(joints,  fractures) 

Relief 


Placement   of   carbonate   rock 
units  relative  to  other  rocks 

Placement  of  recharge  and 
discharge  regions 

Thickness  and  chemical  character 
of  aquifer 

Orientation  and  transmissibility 
of  primary  flow  paths 

Defines  hydraulic  gradients 


(after  White,  1977) 
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accuracy. 

Field  reconnaissance  will  commence  after  the  analy- 
sis of  the  aerial  photographs  and  geologic  map  has  been 
completed;  to  substantiate  or  refute  previous  findings. 
Structural  features  and  lithology  will  be  noted  at  se- 
lected sites.  Also,  the  occurrence  of  surface  features, 
such  as  sinks  and  shallow  shafts,  will  be  recorded. 

In  particular,  field  measurements  of  the  fractures 
will  be  obtained  in  close  proximity  to  the  freshwater 
springs,  Buffalo  and  Antelope.  Fracture  measurements 
will  entail:  the  direction  of  traverse,  the  distance  be- 
tween fractures,  fracture  filling,  fracture  aperture  and 
fracture  orientation  and  dip.  Fracture  analysis  and 
other  field  observations  will  be  used  to  approximate  the 
orientation  and  location  of  the  Sulphur  Fault  in  relation 
to  the  freshwater  springs.  In  addition,  the  determina- 
tion of  the  origin  and  occurrence  of  the  freshwater 
springs  will  be  discussed  as  a  result  of  these  findings 
and  previous  investigations  (Hart,  1972). 

A  lineament  map  of  the  Arbuckle  Group  will  be  con- 
structed from  aerial  photographs  that  depicts  fracture 
orientation  in  relation  to  macro-structure,  such  as 
faults.  This  map  will  be  used  in  conjunction  with  field 
measurements,  structural  geology  of  the  area,  and  the 
potentiometric  survey  to  assess  general  groundwater  flow 
direction  in  the  outcrop  area  (Davis  and  Dewiest,  1966)  . 
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A  structural  contour  map  of  the  surface  elevation  of 
the  Arbuckle  Group  of  the  project  will  be  constructed 
from:  (1)  well  completion,  reports  and  well  logs  submit- 
ted to  the  Oklahoma  Corporation  Commission;  (2)  domestic, 
municipal  and  industrial  water  well  reports  submitted  to 
the  Oklahoma  Water  Resources  Board;  and  (3)  thickness  of 
stratigraphic  units  as  determined  by  previous  investiga- 
tors (Ham  et  al.,  1945,  1954,  and  1978).  The  construc- 
tion of  this  map  should  provide  information  for  the  in- 
ference of  the  Sulphur  Fault  and  other  faults  in  the  area 
that  lie  at  depth  below  the  Pontotoc  Group  in  the  Sulphur 
area.  In  addition,  a  relationship  could  develop  between 
potentiometric  levels  and  the  distribution  of  surface 
elevations  of  the  Arbuckle  Group. 

3.2   Determination  of  the  Physical  Properties  of  the 
Simpson-Arbuckle  Aquifer 

Rock  specimens  will  be  obtained  at  selected  sites 
from  the  project  area  during  the  field  reconnaissance 
phase  and  returned  to  the  petrophysics  laboratory  for  de- 
termination of  effective  porosity  and  thin-section  analy- 
sis. 

Thin  sections  of  representative  lithologic  units 
from  the  study  area  will  be  photographed,  utilizing  an 
Olympus  Zoom  Stereo  Microscope,  Model  SZ-TR,  and  a 
Polaroid  Land  camera.  Magnification  of  each  thin  section 
can  be  obtained  by, 
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M  =  EP  x  PEP  x  Z  (3.1) 

where 

M  =  total  magnification, 
EP  =  Eyepiece  magnification, 
PEP  =  Photo  eyepiece  magnification, 

Z  =  Zoom  magnification. 
Thin  section  analysis  should  reveal  the  general  texture 
and  fabric  of  the  Pontotoc,  Simpson  and  Arbuckle  Groups. 

Effective  porosity  measurements  of  selected  core 
samples,  representing  the  Pontotoc,  Simpson,  and  Arbuckle 
Groups,  will  be  ascertained  by  the  use  of  the  Helium 
Porosimeter.  The  porosity  can  be  determined  by  the  use 
of  Boyles'  Law  (Collins,  1961): 

P1V1  =  P2V2  =  constant  (3.2) 

where 

P.  =  Initial  pressure, 

V.  =  Initial  volume, 

Pp  =  Final  pressure, 

V^  =  Final  volume. 
The  following  equations  were  derived  by  the  author  to  ob- 
tain volume  of  the  pores,  V  ,  for  each  core  sample: 

Pvac  VL  +  Pvac  (A  Voh  "  VB  +  Vp'  +  P1V1 

"  P2  (VL  "  VB  +  Vp  +  A  Vch'  <3-3) 

rearranging  equation  (3.3)  and  solving  for  V  ; 

P           (V_    +  A  V  .     -  V_)    +    (P-.V-)    -  P0(VT    -V     +  AV,     +VJ 
v     _       vac       L ch         B  11  2     L         B ch         1 

(P2  "  Pvac)  (3.4) 
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where 

V  =  pore  volume  of  the  core  sample  (c.c), 
P 

V  =  bulk  volume  of  the  core  sample  (c.c.), 

V  =  volume  of  the  lines  (c.c.), 

Li 

V  .  =  volume  of  the  sample  chamber  (c.c), 

V  =  volume  of  disks  in  the  sample  chamber  (c.c), 

AVch  "  Vch  -  VD  (C-C-'' 

V-  =  volume  of  reference  chamber  number  1  (c.c), 

P     =  initial  vacuum  pressure  of  V_  and  [V_,+AV  u~v0] 
vac  LP    ch   B 

(negative  gage  pressure,  p.s.i), 

P.  =  initial  pressure  of  V-  (positive  gage  pressure, 

p.s.i.) , 

P0  =  final  pressure  of  V. ,  VT  and  [V+AV  ,  -V_] 
2  1    L         P    Ch   B 

(positive  or  negative  gage  pressure,  p.s.i.). 

Finally,   effective  porosity,   n,   can  be  determined  by 

(Collins,  1961)  : 

n  =  V/V_  (3.5) 

P   B 

where  V   is  the  volume  of  the  pores,  and  V_  is  the  bulk 

p  *  '  B 

volume  of  the  core  sample. 

Determination  of  the  storage  coefficient  and  the 
transmissivity  of  the  Arbuckle  Group,  in  close  proximity 
to  the  freshwater  springs,  will  be  determined  by  several 
field  methods  and  by  the  use  of  the  groundwater  model 
devised  by  Prickett  and  Lonnquist  (1975).  In  addition,  a 
regional  technique  utilized  by  Fairchild  et  al.  (1983), 
will  be  employed  to  determine  aquifer  transmissivity  in 
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relation  to  spring  discharge. 

An  attempt  will  be  made  to  determine  the  storage  co- 
efficient of  the  Arbuckle  aquifer  by  a  method  derived 
originally  by  Jacob  (1940)  and  reported  by  Todd  (1959)  . 
This  method  is  based  on  water  level  fluctuations  in  an 
elastic,  confined  aquifer  in  response  to  atmospheric 
pressure  changes. 

As  discussed  previously,  the  barometric  efficiency 

BE  is, 

Ah        1 
BE  =  y   =  (3.6) 

W  AP      a 

a     —  +  1 

ne 

where 

Y   =  specific  weight  of  water, 

Ah  =  change  in  water  level, 

AP   =  change  in  barometric  pressure, 

a  =  aquifer  compressibility  factor, 

n  =  porosity  of  the  aquifer, 

6  =  water  compressibility  factor. 

Also,  the  storage  coefficient  S   for  a  confined  aquifer 

can  be  expressed  as 

S^  =  b'Y„  (a  +  n8)  (3.7) 

c      w 

were  b1  is  the  thickness  of  the  confined  aquifer. 

The  procedure  to  obtain  barometric  efficiency  for 
this  investigation  is  as  follows: 
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A  microbarograph  will  be  acquired  from  the 
Meteorology  Department,  University  of  Oklahoma. 
Then,  the  microbarograph  will  be  positioned  inside 
the  Travertine  District  Nature  Center,  located  in 
close  proximity  to  the  observation  wells.  Contin- 
uous readings  of  barometric  pressure  for  an  entire 
month  will  be  recorded  on  the  strip  charts  of  the 
microbarograph.  The  strip  charts  will  then  be  an- 
alyzed for  pressure  changes  (inches  of  Hg)  during 
incremental  periods  of  time  (hours) .  This  then  will 
be  correlated  with  the  strip  chart  data  of  the 
Steven's  Water-Level  Recorders,  positioned  at  the 
East  and  West  Observation  Wells,  for  the  same  incre- 
mental periods  of  time.   Data  will  then  be  plotted, 

AP    (horizontal  axis)   versus  Ah   (vertical  axis) . 
a 

The  slope  of  the  line  will  yield  the  barometric  ef- 
ficiency.  From  the  known  value  of  barometric  effi- 
ciency, the  storage  coefficient  S  ,  can  be  obtained. 
Slug  tests  will  be  accomplished  at  selected  sites  in 
the  project  area.   This  rate  of  rise  technique  to  deter- 
mine hydraulic   conductivity  K,   for  confined  and  un- 
confined  aquifers,  fully  or  partially  penetrated  wells 
can  be  expressed  as  (Bouwer,  1978)  : 

r  2  In  (Re/r  )   1      y 

K  =  — - 5L-  -   In  —  (3.8) 

2L  t     Yt 

e  t 
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where 

Re  =  effective  radial  distance  (connotes  the  same 
meaning  as  effective  radius  in  well  hydraul- 
ics) , 

r   =  distance  from  the  well  center  to  aquifer, 

w  • 

uncased  section, 

r   =  distance  from  the  well  center  to  the  casing, 

c  ' 

L   =  length  of  the  open  section  of  the  well  where 

water  enters, 

y   =  y  at  time  zero, 

y   =  y  at  time  t, 

t  =  time  since  Y  . 

o 

For  flow  in  partially  penetrating  wells  it  was  determined 

experimentally  that  In  (Re/r  )  can  be  determined  from  an 

empirical  equation  (Bouwer,  1978) : 

Re  1 

In  — = (3.9) 

rw      1.1      ^     A  +  B  In  [  (H-Ifr) /rj 

**  (W  +  w 

where  A  and  B  are  constants  that  can  be  obtained  from  a 
nomograph  devised  by  Bouwer  (1978),  if  (L  /r  )  is  known; 
and  H  is  the  saturated  thickness  of  the  aquifer.  If  the 
In  [  (H-L  ) /r  ]  exceeds  the  value  of  6,  a  value  of  6 
should  still  be  used. 

The  procedure  for  the  slug  tests  is  as  follows: 

A  cylinder  filled  with  sand,  capped  at  both 

ends,  will  be  lowered  into  a  well.   A  pressure 
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transducer  will  be  lowered  into  the  well  to  monitor 
water  level  fluctuations.  The  probe  (capped 
cylinder)  will  remain  submerged  until  the  water 
level  is  stabilized,  as  indicated  by  the  pressure 
transducer.  The  pressure  head  above  the  transducer 
is  recorded  at  this  time.  The  probe  will  be 
withdrawn,  the  water  level  Y,  will  drop  momentarily, 
and  begin  to  rise.  The  pressure  transducer  will 
monitor  the  rise  of  head,  and  the  investigator 
records  time  versus  head  readings.  This  data  can 
then  be  plotted  on  semilog  paper,  Y  (pressure  head 
of  water)  vs.  time  t.  The  hydraulic  conductivity 
can  then  be  calculated  by  inserting  known  values 
from  the  nomograph  devised  by  Bouwer  (19  78)  and  from 
the  semilog  plot,  into  Equations  (3.9)  and  (3.8). 
Disadvantages  of  this  technique  are  that  K  is 
measured  over  a  relatively  small  portion  of  the 
aquifer;  and  is  not  adequate  for  determination  of 
storage  coefficient,  S  . 

The  Prickett-Lonnquist  (1975)  model  will  be  used  to 
simulate  drawdown  of  a  recently  implaced  municipal  well 
of  the  city  of  Sulphur  wellfield.  The  data  furnished  by 
the  firm  that  conducted  the  pump  test  was  constant  dis- 
charge, water  temperature,  duration  of  test,  and  total 
drawdown  at  the  end  of  the  test.   The  approximate  range 

of  unknown  values,  S   and  T,  should  be  obtained  with  the 

c 
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use  of  Prickett-Lonnquist  Model  by  trial  and  error  for  a 
time  t.  The  two  observation  wells,  East  and  West,  are 
located  approximately  3,000  and  2,500  feet,  respectively, 
from  the  municipal  well.  Strip  charts  from  the  Steven's 
Water-Level  Recorders  positioned  at  each  observation  well 
will  be  analyzed  for  drawdown,  if  any,  during  the  period 
the  pump  test  was  conducted. 

The  Prickett-Lonnquist  model  is  capable  of  handling 
the  transient  flow  equation  of  a  confined  aquifer  that  is 
either  isotropic  or  anisotropic,  homogeneous  or  heteroge- 
neous. Production  wells  at  specified  nodes  in  the  pro- 
gram are  considered  to  be  fully  penetrating  the  saturated 
thickness  of  the  confined  aquifer.  It  is  recognized  by 
this  author  that  the  municipal  well  only  partially  pene- 
trates the  Upper  Arbuckle  Group,  resulting  in  a  possible 
vertical  flow  component  at  the  base  of  a  well.  It  is  ob- 
vious that  the  values  obtained  from  this  simulation  are 
rough  estimates  at  best.  Also,  Gringarten  (1982) , 
asserts  that  if  a  well  intersects  a  vertical  fracture, 
lateral  hydraulic  conductivity  estimates  could  be  in 
error  by  as  much  as  40  percent. 

The  following  assumptions  will  be  made  in  the  set-up 
of  the  model. 

(1)  The  aquifer  is  relatively  infinite,  homogeneous 
and  isotropic  in  aerial  extent.  Boundary  con- 
ditions will  not  be  considered. 
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(2)  The  likelihood  of  intercepting  an  interconnect- 
ed vertical  fracture  is  small.  Majority  of  so- 
lution activity  occurs  along  bedding  planes. 
Hence,  lateral  flow  predominates  (Boulton  and 
Streltsova,  1977) . 

(3)  The  homogeneous  model  is  adequate  for  drawdown 
analysis  at  the  sink  node. 

(4)  Transmissivity  values  will  be  correlated  to 
drawdown,  if  any,  at  the  observation  wells. 

(5)  Values  obtained  from  this  analysis  is  represen- 
tative of  only  the  saturated  thickness  pene- 
trated by  the  production  well. 

(6)  The  base  of  the  partially  penetrating  well  will 
be  considered  relatively  impermeable  in  the 
vertical  direction;  hence  a  fully  penetrating 
well. 

The  regional  technique  employed  by  Fairchild  et  al. 
(1983),  will  be  used  to  determine  transmissivity  T,  ex- 
pressed in  the  form  of  Darcy's  Law: 


Q 
T  =  —  (3.10) 

IW 


where 


3 
Q  =  spring  discharge  (ft  /day), 

I  =  average  gradient  of  the  aquifer  to  the  spring 

discharge  point  (ft/ ft) , 
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W  =  average  width  of  the  aquifer  contributing  to  the 
spring  discharge. 
Average  values  of  stage  readings  of  Antelope  and  Buffalo 
will  be  obtained  during  the  potentiometric  survey  of  the 
project  area.  The  average  gradient  of  the  aquifer (s) 
contributing  to  the  spring  discharge  will  be  determined 
from  the  survey.  Likewise,  the  average  width  of  the  con- 
tributing aquifer (s)  to  spring  discharge  will  be  ascer- 
tained from  the  survey. 

3.3   Wells  and  Springs 

Spring  response  to  long-term  precipitation  events 
will  be  analyzed  with  the  use  of  a  66-year  precipitation 
histogram  that  also  depicts  periods  of  no  flow  for 
Buffalo  and  Antelope.  From  this  histogram,  an  inch- 
years,  time  bar  chart  will  be  constructed.  This  chart 
will  depict  the  summation  of  annual  precipitation  events 
as  positive  or  negative  for  a  specified  period  of  time. 
Annual  rainfall  that  exceeds  the  mean-annual  rainfall  is 
considered  positive.  Conversely,  annual  rainfall  that  is 
below  the  mean  annual  is  negative.  The  area  (inch-years) 
will  be  summed  for  a  specified  time  period  with  the  mean 
annual  rainfall  as  the  reference  datum.  The  results  will 
be  graphically  shown  on  the  inch-years  bar  chart.  In- 
cluded in  each  bar  will  be  the  periods  of  no  flow.  It  is 
hoped,  that  this  graphical  illustration  will  demonstrate 
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the  deprivation  effects  of  nearby  production  wells  for 
the  66-year  period. 

Well  and  spring  hydrograph  versus  precipitation 
events  will  be  scrutinized.  Information  that  can  be 
obtained  from  this  analysis  is:  (1)  minimum  piezometric 
elevation  of  observation  wells  before  discharge  of 
springs  commence;  (2)  long-term  trend  hydrograph  of 
observation  wells  to  forecast  possible,  spring-flow 
cessation;  (3)  correlates  the  magnitude  of  discharge  of 
the  springs  and  water-level  fluctuations  of  the  obser- 
vation wells  to  precipitation  events  on  a  seasonal  and 
yearly  basis;  and  (4)  hopefully  provide  a  confirmation  of 
the  approximate  size  of  the  recharge  area  as  determined 
from  the  potentiometric  survey  and  the  analysis  of 
geologic  control. 

A  diffusivity  coefficient,  K,  (months  ) ,  of  the 
aquifer (s)  responsible  for  the  freshwater  springs  will  be 
determined  by  methods  reported  by  E.  White  (1974)  and  W. 
White  (1977).  Also,  an  approximate  response  time,  t 
(days)  ,  will  be  ascertained  to  aid  in  a  concise  classi- 
fication of  the  freshwater  springs  (E.  White,  1974). 

The  recession  limb  of  the  spring  hydrograph  can  be 
expressed  as  a  decaying  exponential  function  (W.White, 
1977) , 


Q  =  QQe  Kd  t  (3.11) 


where 
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Q  =  spring  discharge  (cfs)  at  some  time  t  (mos.) , 

Q   =  initial  spring  discharge  located  at  the 

inflection  point  just  below  the  peak  discharge 

of  the  spring  hydrograph  (cfs) , 

K,  =  diffusivity  coefficient  (mos.   ), 

t  =  time  in  months. 

In  addition,  response  time  t  ,  can  be  obtained  by  the 

following  expression  (E.  White,  1974), 

30 
t   =  —  (3.12) 

r    Kd 

where 

t   =  time  of  response  (days) , 

K,  =  diffusivity  constant  (months   ) . 

Another  parameter,  spring  variance  or  variability  of 
springs,  can  be  used  to  generalize  storage  capacity  of  an 
aquifer  responsible  for  spring  discharge  (Davis  and 
Dewiest,  1966) .  The  spring  variance  V* ,  for  a  specified 
period  of  time  can  be  expressed  as: 


Q    -  Q  • 
v;  =  -SS2 EiS  (3.13) 

med 


where 


Q    =  maximum  recorded  spring  discharge  for  the 
specified  time  period, 

Q  .   =  minimum  recorded  spring  discharge  for  the 

mm  c         3 

specified  time  period, 
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Q   ,  =  median  discharge  for  the  prescribed  time 
med 

period. 
Spring  variance  analysis  will  be  conducted  on  Antelope 
and  Buffalo. 

Tables  and  maps  will  depict  location  and  flow  of 
artesian  and  production  wells  in  close  proximity  to  the 
park.  A  hydrologic  budget  will  be  developed  to  demon- 
strate deprivation  effects  of  artesian  and  production 
wells  in  the  city  of  Sulphur  versus  the  total  annual  re- 
charge rate  in  the  recharge  area. 

A  statistical  analysis  will  be  utilized  to  determine 
average  flow  rate  of  Buffalo  and  Antelope  on  a  monthly 
and  yearly  basis.  This  then  will  be  correlated  with  pre- 
cipitation on  a  monthly  and  yearly  basis.  Definite 
trends  in  fluctuations  of  spring  discharge  should  be  es- 
tablished. 

Downhole  TV  logs  of  two  separate  artesian  wells  in 
the  project  area  will  be  viewed.  Both  surveys  were 
completed  by  the  National  Park  Service  in  August  1980. 
Information  that  might  be  divulged  from  these  tapes 
include:  (1)  flow  direction,  (2)  occurrence  of  solution 
cavities  and  fractures  with  depth,  (3)  depth  of  well  and 
casing,  and  (4)  formation  contacts. 

3.4   Determination  of  Recharge  Areas 

Pervious  investigators  Hart  (1972),  and  Harp  et  al. 
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(1976)  ,  identified  two  separate  recharge  areas  responsi- 
ble for  spring  flow  in  the  park.  Location  of  these  re- 
charge areas  were  given  in  general  terms,  no  specific 
locations  were  cited. 

The  aerial  extent  of  this  potentiometric  survey  will 
be  considerably  refined,  because  of  a  previous  potentio- 
metric survey  (Winter  1976,  1977)  conducted  by  Fairchild 
et  al.  (1983),  of  the  Arbuckle  Mountain  region  (Plate 
II)  .  The  purpose  of  this  potentiometric  survey  is  to 
define  the  recharge  areas  for  the  springs  and  establish  a 
groundwater  divide  for  the  Sulphur  area.  This  is  neces- 
sary in  estimating  the  hydrologic  budget  for  the  Sulphur 
area . 

The  procedure  to  obtain  piezometric  elevations  of 
selected  well  sites  throughout  the  project  are  as  fol- 
lows: 

The  survey  will  occur  during  the  winter  months 
of  1985,  when  precipitation  rates  are  low  and  the 
aquifer  system  is  fairly  stable.  The  survey  will  be 
taken  in  the  span  of  two  weeks,  to  minimize  error  of 
water  level  fluctuations.  Surface  elevations  (feet 
above  sea  level)  will  be  obtained  from  topographic 
maps  at  each  well  site,  accuracy  should  be  within  5 
feet.  Depth  to  water  table  at  each  well  site  will 
be  read  with  chalk  and  steel  tape.  The  piezometric 
level  (feet  above  sea  level)  at  each  well  site,  will 
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be  the  difference  between  the  surface  elevation  and 
the  depth  to  water  table  reading.   Total  depth  of 
well  readings  will  not  be  taken  due  to  lack  of  prop- 
er equipment.   However,  inquiry  with  landowners  will 
be  made  to  determine  approximate  depth  of  well  and 
the  general  stratigraphy  the  drillers  encountered. 
Once  the  survey  is  completed,  a  potentiometric  map 
will  be  constructed   from  the  piezometric  elevations. 
General  flow  direction,  recharge  areas  and  water  table 
divides  should  be  reasonably  ascertained.   This  informa- 
tion will  then  be  depicted  on  the  geologic  and  hydro- 
geologic  map  of  the  study  area. 


CHAPTER  IV 


FINDINGS 


4 . 1   Geology 

A  small  outcrop  or  fenster  of  dolomitic  carbonates 
in  the  midst  of  the  Pontotoc  Group  was  discovered  just 
east  of  the  Travertine  District,  Chickasaw  National 
Recreational  Area  (Plate  III)  .  The  outcrop  was  mapped , 
attitudes  (dip  and  strike)  of  beds  were  taken,  and  a  gen- 
eral lithologic  description  was  accomplished  (Appendix 
III) .  Rock  specimens  collected  from  this  area  were  shown 
to  Dr.  Robert  0.  Fay,  Oklahoma  Geological  Survey,  who  has 
extensive  field  experience  in  this  region.  He  identified 
the  rocks  with  reasonable  certainty  as  being  from  the 
Bromide  Massive  and  the  Bromide  Basal,  the  upper-most 
formation  of  the  Simpson  Group. 

Photogeologic  analysis  of  the  study  area  revealed 
the  following  information.  The  heavily  fractured 
Arbuckle  Group  with  numerous  surface  outcrops  can  be 
readily  distinguished  from  the  Pontotoc  and  the  Simpson 
Groups  in  the  study  area.  The  Simpson  is  usually  detect- 
ed where  beds  dip   steeply,   otherwise   soil   cover   is 
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present,  especially  where  the  basal  sands  outcrop  close 
to  the  surface  (Ham,  1945)  .  Lineaments  in  the  Simpson 
are  few  and  well  spaced.  Some  lineaments  were  detected 
in  the  Simpson  Group,  vicinity  of  Hickory.  Small,  dark, 
circular-shaped  areas  were  observed  in  the  stream  valleys 
of  the  Arbuckle  Group,  especially  concentrated  in  Section 
4,  Southwest  1/4  Section,  T1S,  R4E.  The  Pontotoc  Group 
is  heavily  dissected  with  intermittent  streams,  and  ap- 
pears white  in  color  where  it  outcrops  along  the  periph- 
ery of  the  Arbuckle  Group,  east  of  Sulphur.  A  strong 
occurrence  of  lineation  was  noted,  east-southeast  of  the 
park,  in  the  Pontotoc  Group  in  NEj  SEj  of  Sec.  1,  T1S, 
R3E;  and  NWj  SWj  of  Sec.  6,  T1S  R4E. 

A  lineament  map  was  constructed  from  lineations  in 
rock  outcrops,  observed  from  the  black  and  white  aerial 
photographs  of  the  project  area  (Plate  IV) .  The  map  was 
completed  to  facilitate  correlation  between  field  recon- 
naissance of  fractures  and  those  observed  from  aerial 
photos,  and  to  demonstrate  linear  trends  with  respect  to 
macrostructures ,  such  as  faults.  Favorable  correlation 
exists  between  orientation  of  lineaments  and  trends  of 
faults  (Plate  IV) . 

Interestingly  enough,  there  are  few  recognizable 
lineaments  other  than  stream  alignment  along  the  Sulphur 
Fault.  In  fact,  one  cannot  discern  the  exact  location  of 
Sulphur   Fault   from   the   aerial   photos,   unless   rock 
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outcrops  of  the  Simpson  and  Arbuckle  are  present.  Lack 
of  lineation  can  be  attributed  to  accelerated  weathering 
of  the  regolith  because  of  intense  fracturing  in  the 
fault  zone  (Fetter,  1980)  .  In  addition,  inspection  of 
the  aerial  photographs  revealed  that  intermittent  stream 
direction  and  preferred  orientation  of  lineaments  at  var- 
ious locations  are  nearly  identical. 

Fracture  measurements  were  taken  at  Antelope  Spring 
when  it  ceased  flowing  (Table  9)  .  Fractures  are  exposed 
along  a  west-facing  escarpment  of  broken  conglomerate, 
approximately  20  feet  high  and  20  feet  wide.  Later  in 
the  investigation,  spring  flow  was  observed  to  originate 
from  the  base  of  the  west-facing  escarpment  at  the  north 
end. 

Attitude  of  the  Pontotoc  beds  are  nearly  horizontal 
at  or  near  Antelope  Spring.  Orthogonal  fracture  sets 
were  noted.  Major  trends  of  fractures  are  S70°W  and 
N25°W.  Dips  of  fracture  planes  range  from  66°  to  verti- 
cal. Fracture  width  ranges  from  .25  inches  to  3  inches 
at  the  surface  outcrop. 

There  were  no  fracture  measurements  taken  at  Buffalo 
Spring.  Soil  cover  is  well  established  in  and  around  the 
spring  area.  A  man-made,  circular  retaining  wall  con- 
sisting of  native  rock,  encompasses  a  gravel  bed  that 
contains  mainly  lime  cobbles,  and  some  chert  and  quartz 
sand;  from  which  Buffalo  Spring  eminates  and  flows  over  a 
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Table  9.   Fracture  Measurements  of  the  Pontotoc  Group, 
Antelope  Spring* 


FRACTURE 
DIP       WIDTH 
STATION   ORIENTATION   (degrees)   (inches) 


DISTANCE 
BETWEEN 
STATIONS 
(feet) 


TOTAL 

DISTANCE 

(feet) 


1 

S70E 

2 

S27E 

3 

N70E 

TRAVERSE  DIRECTION  IS  SOUTH 
72S        .5  0. 

66S        .5-2.0       10.0 
82N       2.0  13.0 


0. 
10.0 
23.0 


1 

N40E 

81S 

.5 

2 

N70E 

90S 

.25 

3 

S25E 

85S 

.25-. 5 

4 

S52E 

78S 

.5-1.0 

ONE 

READING 

1 

N55E 

84N 

.5-2.0 

TRAVERSE  DIRECTION  IS  S35E 

0. 
2.0 
17.0 
2.0 


0. 

2.0 
19.0 
21.0 


0. 


*  Location  -  SEj  NWj  NWi  NEj  of  Sec.  1,  T1S,  R3E 
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weir-like  opening  (4.5  feet  wide)  into  a  small  unnamed 
creek . 

After  preliminary  photogeologic  analysis,  fracture 
measurements  were  taken  in  the  Pontotoc  Group,  east  of 
the  park  boundary.  The  fractures  that  were  measured  did 
not  coincide  with  the  exact  location  of  those  observed  on 
the  aerial  photos.  Upon  completion  of  the  fracture 
measurement,  closer  inspection  of  the  site  on  aerial 
photos  depicted  subtle  lineament  trends  that  corresponded 
exactly  with  those  measured  (Table  10)  .  Distance  between 
fractures  in  the  Pontotoc  Group  at  this  site,  range  from 
.5  to  36  feet.  Preferred  orientation  of  fractures  are 
west-northwest  (N60W  to  W) . 

Fracture  measurements  were  taken  of  a  durable 
calcite  cemented  sandstone,  approximately  1  to  2  feet 
thick,  exposed  in  the  glass-sand  quarry  (Table  11).  The 
bed  is  located  in  the  friable,  basal  sand  of  the  Oil 
Creek  Formation,  Simpson  Group.  The  bulk  of  the  frac- 
tures trend  west-northwest  (N60W  to  W)  and  southwest 
(S30W  to  S60W) .  The  macrostructure  of  the  area,  two  ma- 
jor faults  just  north  of  the  quarry,  Section  11  and  12  of 
T1S,  R4E,  trend  in  the  same  direction  as  the  preferred 
orientation  of  the  fractures.  The  distance  between  frac- 
tures range  from  0  to  6  feet.  Fracture  width  ranged  from 
.0625  to  1.0  inches.  Attitude  of  the  beds  at  this 
locale,  strike  at  N55W,  with  a  dip  of  7  degrees  to  the 
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Table  10.   Fracture  Measurements  of  the 
Pontotoc  Group* 


STATION 

ORIENTATION 

FRACTURE 
DIP       WIDTH 
(degrees)   (inches) 

DISTANCE 
BETWEEN 
STATIONS 
(feet) 

TOTAL 
DISTANCE 
(feet) 

TRAVERSE 

DIRECTION  IS 

N5E 

1 

N70W 

UNK 

+1.0- 

-2.0 

0. 

0. 

2 

N75W 

it 

II 

2.0 

2.0 

3 

N75W 

it 

II 

2.0 

4.0 

4 

S87W 

it 

It 

1.5 

5.5 

5 

N72W 

it 

It 

1.5 

7.0 

6 

N60W 

it 

tt 

.5 

7.5 

7 

N60W 

ii 

tl 

2.0 

9.5 

8 

W 

it 

II 

3.0 

12.5 

9 

N88W 

it 

It 

1.5 

14.0 

10 

N64W 

tt 

2.0 

3.0 

17.0 

11 

N75W 

ii 

.5- 

■1.5 

6.0 

23.0 

12 

N65W 

it 

+1.0- 

•2.0 

8.0 

31.0 

13 

N50W 

it 

II 

36.0 

67.0 

14 

N75W 

ti 

II 

12.0 

79.0 

15 

N80W 

ii 

II 

4.0 

83.0 

16 

S89W 

it 

II 

26.0 

109.0 

17 

N88W 

it 

It 

10.0 

119.0 

18 

N87W 

it 

II 

6.0 

125.0 

19 

N72W 

ti 

tl 

10.0 

135.0 

*  Location  -  SEj  NEi  NEj  of  Sec.  1,  T1S,  R3E 
Fracture  is  filled  with  soil  and  vegetation. 
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Table  11.   Fracture  Measurements  of  an  Indurated 

Sandstone  of  the  Simpson  Group  Glass  Sand 
Quarry* 


STATION    ORIENTATION 


DIP 
(degrees) 


FRACTURE 
WIDTH 
(inches) 


DISTANCE 
BETWEEN 
STATIONS 
(feet) 


TOTAL 
DISTANCE 
(feet) 


TRAVERSE  DIRECTION  IS  N25W 


1 

S55W 

UNK 

.5-1.0 

0. 

0. 

2 

S70W 

ii 

.25 

3.5 

3.5 

3 

S60W 

ii 

.25 

3.5 

7.0 

4 

S62W 

ii 

.125 

3.5 

10.5 

5 

S25W 

it 

.25 

3.0 

13.5 

5a 

S67W 

ii 

.25 

0. 

13.5 

6 

S75W 

•i 

.125 

4.0 

17.5 

6a 

N70W 

ii 

.125 

0. 

17.5 

7 

S70W 

ii 

.125 

4.0 

21.5 

7a 

N30W 

ii 

.25 

0. 

21.5 

8 

N80W 

n 

.25 

6.0 

27.5 

8a 

S55W 

ii 

.25 

0. 

27.5 

9 

N80W 

ii 

.25 

3.0 

30.5 

10 

S45W 

it 

.25 

2.0 

32.5 

11 

N75W 

ii 

.125 

4.5 

37.0 

12 

S70W 

ii 

.125 

4.0 

41.0 

12a 

N40W 

ii 

.125 

0. 

41.0 

13 

S50W 

ii 

.25-. 5 

3.0 

44.0 

13a 

N80W 

ti 

.5 

0. 

44.0 

14 

S50W 

ii 

.5 

1.5 

45.5 

15 

S35W 

it 

.125 

6.0 

51.5 

15a 

N82W 

it 

.125 

0. 

51.5 

16 

N80W 

ii 

.125 

6.0 

57.5 

17 

S55W 

■t 

.25-. 5 

1.5 

59.0 

18 

N80W 

it 

.25 

1.5 

60.5 

19 

N80W 

it 

.75 

6.0 

66.5 

19a 

N45W 

it 

.75 

0. 

66.5 

20 

S35W 

ii 

.25-. 5 

15.0 

81.5 

20a 

N82W 

ti 

.125 

0. 

81.5 

21 

S50W 

ii 

.25-. 5 

3.0 

84.5 

21a 

N30W 

it 

.125 

0. 

84.5 

22 

N75W 

it 

.125 

10.0 

94.5 

22a 

S15W 

•i 

.0625 

0. 

94.5 

23 

N10W 

n 

.125 

6.0 

100.5 

24 

S40W 

it 

.25 

3.0 

103.5 

*  Location  -  SWj  SWj  NW£  NEj  of  Sec.  14,  T1S,  R4E 
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northeast. 

The  field  reconnoiter  of  the  Arbuckle  Group  revealed 
intensely  fractured  dolomites.  A  distance  of  one  inch  or 
less  between  fractures  is  not  uncommon.  Most  of  the 
small  fractures  (less  than  .125  inches  thick)  at  the  sur- 
face are  filled  with  limonite  and  silica.  Fracture  width 
ranges  from  .0625  to  3  inches.  A  fair  amount  of  swallow 
shafts  and  small  sinks  are  located  in  and  along 
intermittent  streams. 

The  primary  directions  of  fracture  orientation  in 
the  Arbuckle  Group  measured  at  3  sites,  (Tables  12 ,  13 
and  14)  correlates  extremely  well  with  the  lineament  map 
(Plate  IV)  .  The  majority  of  the  fracture  measurements 
were  oriented  in  the  north-northwest  (N30W  to  N)  and  the 
south-southwest  (S60W  to  W)  directions.  The  bulk  of  the 
stream  directions  in  the  project  area,  east  of  Sulphur, 
are  aligned  with  the  preferred  orientation  of  the  frac- 
tures. The  beds  of  the  Arbuckle  Group,  east  of  Sulphur, 
gently  dip  to  the  northwest  (Plate  I) . 

Generally,  other  than  the  small  outcrop  that  was 
discovered  east  of  the  park  and  nomenclature  of  the  con- 
fining bed,  the  Pontotoc  group;  the  geologic  structure 
and  the  stratigraphy  of  the  area  defined  by  previous  in- 
vestigations versus  this  investigation,  are  in  good 
agreement. 

A  structural  contour  map  of  the  surface  of  the 
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Table  12.   Fracture  Measurements  of  the  Arbuckle  Group, 
Mill  Creek  and  State  Highv/ay  7* 


STATION    ORIENTATION 


DIP 
(degrees) 


FRACTURE 
WIDTH 
(inches) 


DISTANCE 
BETWEEN 
STATIONS 
(feet) 


TOTAL 
DISTANCE 
(feet) 


TRAVERSE  DIRECTION  IS  N70W 


near 

1 

S60W 

vertical 

1.0 

0. 

0. 

2 

S70W 

ii 

1.5 

.5 

.5 

3 

S70W 

ii 

1.5 

.5 

1.0 

4 

S40W 

ii 

.5 

1.0 

2.0 

5 

N5W 

it 

.5 

2.0 

4.0 

6 

N35W 

it 

1.5 

2.0 

6.0 

7 

S60W 

ti 

.5 

.5 

6.5 

8 

N50W 

ii 

.5 

1.0 

7.5 

8a 

N30W 

H 

.25 

0. 

7.5 

8b 

S70W 

H 

.25 

0. 

7.5 

9 

S5W 

ii 

.75 

0. 

8.5 

9a 

S85W 

it 

.75 

0. 

8.5 

10 

S60W 

it 

.125 

.5 

9.0 

10a 

S10W 

tt 

.125 

0. 

9.0 

11 

S 

it 

1.0 

.5 

9.5 

11a 

S80W 

it 

.5 

0. 

9.5 

lib 

S10W 

it 

.25 

0. 

9.5 

12 

S45W 

H 

2.0 

2.0 

11.5 

13 

S75W 

ii 

1.0 

1.0 

12.5 

14 

N40W 

tt 

1.0 

1.0 

13.5 

14a 

W 

it 

2.0 

0. 

13.5 

15 

S45W 

it 

1.5 

.5 

14.0 

15a 

N15W 

it 

2.0 

0. 

14.0 

16 

S80W 

it 

2.0 

1.0 

15.0 

17 

N70W 

tt 

2.0 

2.0 

17.0 

18 

N50W 

it 

1.0 

1.0 

15.0 

19 

N70W 

it 

2.0 

2.0 

20.0 

20 

S70W 

it 

3.0 

5.0 

25.0 

20a 

N15W 

it 

1.0 

0. 

25.0 

21 

N15W 

ti 

1.0 

1.0 

26.0 

21a 

S60W 

ii 

1.0 

0. 

26.0 

22 

S40W 

ti 

2.0 

4.0 

30.0 

22a 

N20W 

it 

2.0 

0. 

30.0 

*  Location  -  SE£  NEj  SEi  SEj  of  Sec.  32,  TIN,  R4E 
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Table  13.   Fracture  Measurements  of  the  Arbuckle  Group, 
Intersection  of  State  Highways  7  and  12* 


STATION    ORIENTATION 


DIP 
(degrees) 


FRACTURE 
WIDTH 
(inches) 


DISTANCE 
BETWEEN 
STATIONS 
(feet) 


TOTAL 
DISTANCE 
(feet) 


TRAVERSE  DIRECTION  IS  S30W 


near 

1 

S65W 

vertical 

.125 

0. 

0. 

2 

N20W 

H 

.25 

.1 

.1 

3 

S65W 

n 

.125 

.1 

.2 

4 

S70W 

ii 

.125 

1.0 

1.2 

5 

S75W 

ii 

.125 

.1 

1.3 

6 

S60W 

ii 

.25 

.1 

1.4 

7 

S65W 

ii 

.0625 

.6 

2.0 

8 

S60W 

ii 

.0625 

.1 

2.1 

9 

W 

ii 

.0625 

.2 

2.3 

10 

S60W 

ii 

.125 

.1 

2.4 

11 

S65W 

it 

.125 

.2 

2.6 

12 

S55W 

ii 

.125 

.2 

2.8 

13 

S 

ii 

.25 

.5 

3.3 

13a 

S65W 

H 

.125 

0. 

3.3 

14 

N25W 

ii 

.5 

.2 

3.5 

TRAVERSE  DIRECTION  IS  W 


near 

1 

N20W 

vertical 

.125 

0. 

0. 

2 

N40W 

ii 

.25 

.3 

.3 

3 

N30W 

H 

.25 

.3 

.6 

4 

N30W 

ii 

.25 

.6 

1.2 

4a 

N10W 

H 

.5 

0. 

1.2 

4b 

S40W 

H 

1.0 

0. 

1.2 

*  Location  -  NEj  NEj  NWj  NW|  of  Sec.  2,  T1S,  R4E 
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Table  14.   Fracture  Measurements  of  the  Arbuckle  Group, 
Collins'  Ranch* 


STATION 


DIP 
ORIENTATION   (degrees) 


FRACTURE 
WIDTH 
(inches) 


DISTANCE 
BETWEEN 
STATIONS 
(feet) 


TOTAL 
DISTANCE 
(feet) 


TRAVERSE  DIRECTION  IS  WEST 


near 

1 

S35 

vertical 

.125 

0. 

0. 

la 

N40W 

.125 

0. 

0. 

2 

S40W 

.25 

.25 

.25 

2a 

N 

.125 

0. 

.25 

3 

S30W 

.125 

.25 

.50 

4 

S20W 

.0625 

.1 

.60 

5 

S55W 

.125 

.2 

.80 

6 

N70W 

.0625 

.2 

1.0 

6a 

S45W 

.125 

0. 

1.0 

TRAVERSE  DIRECTION  IS  N10W 


near 

1 

S80W 

vertical 

1.5 

0. 

0. 

la 

S55W 

■I 

1.5 

0. 

0. 

2 

S70W 

ii 

1.0 

.8 

.8 

2a 

S5W 

ii 

.25 

0. 

.8 

3 

S45W 

ii 

.25 

.7 

1.5 

*  Location  -  NW*  NEj  NEj  of  Sec.  6,  T1S,  R4E 
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Arbuckle  Group  of  the  Sulphur  area  was  constructed  to 
infer  subsurface  structure  under  the  confining  layers  of 
the  Pontotoc  (Figure  29)  .  Inferrence  of  the  Sulphur 
Fault  is  based  on:  spring  occurrence,  observation  well 
lithology,  comparatively  large  discharge  of  the  Vendome 
Well,  general  structural  trends  of  the  area,  and  eleva- 
tion tops  of  the  Arbuckle  Group  from  well  completion  re- 
ports and  well  logs.  South  Hickory  Fault  was  inferred 
with  reasonable  assurance  to  Sec.  20  of  TIN,  R3E  (Figure 
29)  .  This  was  based  on  a  reported  fault  in  the  SW§  of 
Sec.  20  of  the  basal  Simpson  and  Upper  Arbuckle  from  a 
well  completion  report.  An  unnamed  fault,  now  called 
South  Sulphur  Fault  for  this  investigation,  is  inferred 
to  the  Bromide  Spring  area  (Figures  29  and  30) . 

Only  a  small  veneer  of  Simpson  overlain  by  the 
Pontotoc  is  present  on  the  upthrown  side  of  the  Sulphur 
Fault;  except  for  a  thicker  sequence  of  Simpson  in  the 
drag  fold  of  the  down-thrown  side  of  the  inferred  South 
Hickory  Fault,  north  of  Sulphur  (Figure  29) . 

The  thickness  of  the  Pontotoc  increases  in  a  wester- 
ly and  northerly  direction;  ranging  from  less  than  one 
foot  at  the  contact  with  the  Arbuckle  Group  in  Sec.  6  of 
T1S,  R4E,  to  900  feet  thick  in  Sec.  31  of  TIN,  R3E.  In- 
formation on  stratigraphic  thickness  was  obtained  from 
well  completion  reports. 

Finally,  a  geologic  and  hydrogeologic  map  was  con- 
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Figure    30.      Geologic   and  hydrogeologic  map  of   the 
Sulphur,    Oklahoma,   Area. 
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structed  (Figure  30)  .  The  geology  is  modified  after  Ham 
(1954).  Inferrence  of  the  Sulphur  Fault,  with  reasonable 
certainty,  is  depicted  in  Figure  30.  Also,  the  recently 
discovered  outcrop  of  the  Simpson,  located  southeast  of 
the  freshwater  springs,  is  shown  on  the  map.  A  geologic 
cross  section  of  the  Sulphur  Syncline  and  the  area  of  the 
springs  is  depicted  also  in  Figure  30. 

4.2   Petrophysical  Results 

Effective  porosity  measurements  with  the  use  of  the 
Helium  Porosimeter  were  conducted  on  selected  rock  speci- 
mens from  the  project  area  (Table  15)  .  The  primary  po- 
rosity of  the  Arbuckle  Group  ranged  from  7.4  to  10.4%. 
The  confining  layer  for  the  Sulphur  artesian  basin,  the 
Pontotoc  conglomerate,  had  extremely  low  porosities, 
ranging  from  .13  to  2.5%.  The  porosities  for  the  Simpson 
Group  ranged  from  2.3  to  11.4%.  Thin  sections  of  select- 
ed rock  specimens  of  the  Arbuckle  Group,  Simpson  Group, 
and  the  Pontotoc  Group  are  depicted  in  Figures  31  through 
34.  Classification  schemes  of  thin  sections  are  in  ac- 
cordance with  Williams  et  al.  (1954). 

Analysis  of  the  specimens  of  the  Arbuckle  Group  re- 
vealed a  mosaic  of  aphanitic  dolomite  anhedra .  Selected 
specimens  from  the  Simpson  Group  included  poorly-sorted, 
calcite-cemented  calcarenites ;  and  calcite-cemented, 
fairly  well  sorted,  medium-grained,  subangular  to  sub- 
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Figure  31 .   Magnification  is  39X;  mosaic  of  aphanitic 
dolomite  anhedra  of  the  West  Spring  Forma- 
tion, Arbuckle  Group.    Specimen  obtained 
from  surface  outcrop  at  C  NE%  NW%  NW%  of 
Sec.  2,  T1S,  R4E. 
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Figure  32 . 


Magnification  is  39X;  poorly-sorted, 
calcite-cemented  calcarenite  of  the 
Oil  Creek  Formation,  Simpson  Group. 

Specimen  obtained  from  quarry  at 
C  NE%  SW%  NE^  Of  Sec.  14,  T1S ,  R4E. 
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Figure  33.   Magnification  is  39X;  calcite-cemented, 
fairly  well  sorted,  subangular  to  sub- 
rounded  arenite  of  the  Basal  Oil  Creek 
Formation,  Simpson  Group.    Specimen 
obtained  from  quarry  at  C  NE%  SW%  NE%  of 
Section  14,  T1S,  R4E. 
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Figure  34  .   Magnification  is  39X;  conglomerate  of  the 
Pontotoc  Group  consisting  of  some  fine  to 
medium  grained  quartz  particles,  and  fine 
grained  to  *pebble-size  lime  clasts,  em- 
bedded in  a  calcite  and  argillaceous 
cement.    Specimen  obtained  from  NW%  NW% 
NW%  NW%  Of  Sec.  2,  T1S,  R3E. 

It  is  not  shown  on  thin  section,  but  quite  apparent 
in  the  core  sample. 
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rounded,  arenites.  Analysis  of  the  specimen  of  the 
Pontotoc  Group  revealed  a  conglomerate  consisting  of  some 
subangular  to  subrounded,  fine  to  medium  grained  quartz 
particles  and  fine  grained  to  pebble-size  lime  clasts 
embedded  in  a  calcite  and  argillaceous  cement. 
Pebble-size  lime  clasts  were  not  observed  in  the  thin 
section,  but  were  quite  evident  from  the  core  sample. 

4.3   Assessment  of  Aquifer  Properties 

A  slug  test  was  performed  on  an  abandoned  wildcat 
well  in  the  Arbuckle  Group,  located  in  C  NEJ  SEj  SE£  of 
Sec.  31,  TIN,  R4E .  Dimensions  of  the  well  were  obtained 
from  a  well  completion  report  dated  21  April  1964.  The 
depth  to  water  table  reading  was  105  feet.  Total  well 
depth  is  237  feet.  The  method  employed  was  the  sudden 
withdrawal  of  a  probe  and  then  measuring  the  rate  of  rise 
with  a  pressure  transducer.  The  recorded  data  is  listed 
in  Table  16.  The  data  was  then  plotted  on  semilog  paper 
and  a  hydraulic  conductivity  for  a  partially  penetrating 
well  in  an  unconfined  aquifer  was  calculated  (Figure  35) . 
The  average  value  of  hydraulic  conductivity  was 
determined  to  be  1.6  ft/day  for  an  uncased,  unmodified, 
open-well  length  of  79  feet  (Figure  35) . 

In  addition,  slug  tests  were  performed  on  the  two 
observation  wells,  East  and  West,  that  were  completed  in 
the  Arbuckle  Group  (Plate  III) .   The  results  proved  nega- 
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Table  16.   Slug  Test  Data  of  Crawford's  Well' 


Time 

Pressure  Head  P 

+Head  of  Water 

Rise 

t  (sec) 

P  (psi) 

H  (ft) 

Y(ft) 

6.78 

23.69 

54.68 

2.56 

15.26 

23.79 

54.91 

2.33 

21.76 

23.85 

55.05 

2.19 

31.65 

23.93 

55.23 

2.01 

40.83 

24.00 

55.39 

1.85 

50.64 

24.07 

55.55 

1.69 

58.61 

24.12 

55.67 

1.57 

65.57 

24.16 

55.76 

1.48 

74.80 

24.21 

55.88 

1.36 

83.60 

24.25 

55.97 

1.27 

92.78 

24.29 

56.06 

1.18 

111.03 

24.36 

56.22 

1.02 

124.48 

24.41 

56.34 

.90 

125.60 

24.45 

56.43 

.81 

149.15 

24.49 

•  56.52 

.72 

160.52 

24.52 

56.39 

.65 

175.35 

24.56 

56.68 

.56 

190.38 

24.59 

56.75 

.49 

201.70 

24.61 

56.80 

.44 

220.73 

24.64 

56.87 

.37 

234.39 

24.66 

56.92 

.32 

245.56 

24.67 

56.94 

.30 

272.99 

24.69 

56.98 

.26 

308.71 

24.71 

57.03 

.21 

337.74 

24.73 

57.08 

.16 

469.83 

24.74 

57.10 

.14 

630.00 

24.75 

57.12 

.12 

819.75 

24.75 

57.12 

.12 

Note:   Initial  or  stabilized  pressure  is  24.80  psi  or 

57.24  feet  of  water  above  the  transducer  prior  to 
withdrawal  of  the  probe. 

*  Location  -  C  NEj  SE£  SEj  of  Sec.  31,  TIN,  R4E 
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tive  for  both  wells.  The  rate  of  rise  of  the  West  Obser- 
vation Well  was  extremely  slow  and  actually  stabilizes  at 
a  depth  of  one  foot  below  the  original  piezometric  level. 

The  rate  of  rise  of  the  East  Observation  Well  was 
also  extremely  slow.  As  a  result,  data  was  determined  to 
be  unreliable.  In  addition,  naturally  occurring  asphalt 
in  this  well  could  conceivably  obstruct  hydraulic  conduc- 
tivity tests. 

The  Prickett-Lonnquist  model  (1975)  was  used  to  sim- 
ulate drawdown  at  a  sink  node  within  a  specified  period 
of  time  (Figure  36).  Unknown  parameters,  storage  coeffi- 
cient and  transmissivity ,  were  entered  on  a  trial  and  er- 
ror basis  until  reasonable  values  were  attained  that  co- 
incided with  drawdown  at  designated  nodes  within  a  spec- 
ified period  of  time.  A  10  x  10  array  with  node  spacing 
of  2,640  ft.  (.5  mile)  was  utilized.  NSTEP  was  chosen  to 
be  10,  which  gives  a  total  time  of  25.92  hours;  closely 
approximating  24  hours. 

The  artesian  well  selected  for  simulation  is  located 
approximately  3,000  feet  from  the  West  Observation  Well 
and  2,500  feet  from  the  East  Observation  Well.  This 
flowing  well  is  a  municipal  water  well  for  the  city  of 
Sulphur,  and  flows  at  approximately  150  gallons  per 
minute  at  a  surface  elevation  of  1,060  feet  above  sea 
level.  Upon  completion  of  the  well,  the  drilling  company 
conducted  a  pump  test  (probably  a  specific  capacity  test) 
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Figure  36. 


Prickett-Lonnquist  Model  (1971)  that 
was  used  to  simulate  drawdown  at  a 
sink  node. 
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Figure  36  .   Prickett-Lonnquist  Model  (1971)  that 
was  used  to  simulate  drawdown  at  a 
sink  node  (continued) . 
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and  submitted  the  appropriate  data  on  a  Multi-Purpose 
Water  Well  Report  to  the  Oklahoma  Water  Resources  Board 

(Figure  37)  .   This  data  was  presumed  to  be  correct  and 
was  used  as  the  basis  for  the  simulation. 

Strip  charts  from  the  Stevens  Water-Level  Recorders 
positioned  at  each  observation  well,  for  the  months  of 
March  and  April,  calendar  year  1983,  were  acquired  from 
the  National  Park  Service.  These  charts  will  be  used  in 
quantifying  values  of  transmissivity  and  storage  coeffi- 
cient in  relation  to  drawdown  around  the  sink  node. 

Results  of  the  simulation  of  drawdown  of  85  feet  in 
24  hours  and  analysis  of  the  strip  charts  are  as  follows 

(Table  17  and  18)  : 

(1)  Drawdown  in  and  around  the  well  bore  is  primar- 
ily a  function  of  the  storage  coefficient  of  an 
aquifer. 

(2)  Storage  coefficient,  S  ,  was  determined  to  be 

2  2 

approximately  1,200  ft   per  .25  mi   or  1.7  x 

io-4. 

(3)  The  drawdown  at  the  West  Observation  Well  as  a 
result  of  the  pump  test  was  nil.  The  drawdown 
at  the  East  Observation  Well  could  not  be  dis- 
cerned because  the  recorder  was  inoperative 
part  of  the  time. 

(4)  Transmissivity,  T,  was  determined  to  be  between 
.1  ft2/hour  (Table  17)  to  1  ft2/hour  (Table  18) 
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Figure    37.      Multi-purpose  water  well   report  on 

Municipal  Well   No.    9,    city   of   Sulphur 
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Average  hydraulic  conductivity  can  be  estimated  for 
the  open  length  of  the  well  bore  (510  ft) ,  by  using  the 

expression  T  =  Kh1.   Values  for  hydraulic  conductivity 

-3  -2 

range  from  4.7  x  10    ft/day  to  4.7  x  10    ft/day.   In 

addition,  specific  capacity  of  the  well  was  determined  to 
be  6.3  gallons  per  minute  per  foot  of  drawdown. 

Questionable  results  were  obtained  for  the  deter- 
mination of  storage  coefficient  of  a  confined  aquifer  due 
to  meteorological  phenomena.  The  strip  charts  from  the 
water-level  recorders  at  each  observation  well  and  the 
microbarograph  (for  the  month  of  January,  1985)  were  an- 
alyzed for  water  level  fluctuations  h  (feet),  and  baro- 
metric pressure  changes  P   (inches  of  Hg)  ,  respectively, 

a. 

for  the  same  time  periods.  Water  level  fluctuations  re- 
corded on  the  strip  chart  (15-19  January,  1985)  for  the 
East  Observation  Well  were  used  (Figure  38)  .  However, 
there  appears  to  be  very  little  correlation  between  water 
level  fluctuations  and  barometric  pressure  changes  for 
the  West  Observation  Well  (Figure  39) .  As  a  result,  data 
from  the  West  well  was  not  used.  Also,  barometric  pres- 
sures were  scrutinized  for  the  time  period  15-19  January, 
1985  (Figure  40)  .  Finally,  the  data  from  both  strip 
charts  were  synthesized  and  listed  in  Table  19.  Values 
from  Table  19,  AP  versus  Ah  (vertical)  axis  were  plotted 
(Figure  41) .  The  plot  is  extremely  scattered,  but  a  best 
•eye1  fit  resulted  in  a  value  of  .47  for  barometric  effi- 
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Figure  38 . 


Strip  chart  from  a  Stevens1 Water  Level 
Recorder  positioned  at  the  East  Obser- 
vation Well,  depicting  continuous  mon- 
itoring of  water  level  fluctuations, 
h  (feet) . 
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Figure  39. 


Strip  chart  from  a  Stevens ■  Water 
Level  Recorder  positioned  at  the 
West  Observation  Well,  depicting 
continuous  monitoring  of  water- 
level  fluctuations,  h  (feet). 
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CHART  NO.  5-1071-SD 

BELFORT    INSTRUMENT   COMPANY 
BALTIMORE  MARYLAND     U  S  A. 


instrument  NO0R2<&r.:£!%£    date/: 

REMARKS      _.">"!  SC"V   . -,<■  '\.'~       ?V        CC. 


.££..       STATION    ><*.  'A7^a  '■     *:..'. 

?0  '-£ 


(continued) 
Figure  40. 


Strip  charts  from  a  mi crobaro graph  positioned 

at  the  Travertine  Nature  Center,  depicting 

continuous  readings  of  barometric  pressure, 

P   (inches  of  H  ) . 
a  g 
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CHART  NO.  5-1071-SD 

BELFORT    INSTRUMENT   COMPANY 
BALTIMORE  MARYLAND.   U  S  A. 

INSTRUMENT   \0leUb9.-  «<?<:  ..<?r DATE   /-/«,  -f~        STATION  AJ*JvO*    Ct*s.X?'.>. 

REMARKS     Z.rAttT..Ht!*'.T...  ?**«.<».»*         .SO.  Wo  


Figure  30.   Strip  charts  from  a  mi crobaro graph  position- 
ed at  the  Travertine  Nature  Center,  depict- 
ing continuous  readings  of  barometric  pressure 

P„  (inches  of  H  ) (continued) . 
a  g 
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ciency,  BE. 

Assumptions  that  were  made  in  determining  the  value 
of  the  storage  coefficient  are:   (1)  saturated  thickness 

of  the  Arbuckle  Group  is  3,500  feet;  (2)  water  compressi- 

-7    2 
bility  8,  is  constant  with  a  value  of  .25  x  10    ft  /lb; 

and  (3)  average  porosity  n,  of  the  Arbuckle  Group  is  .1. 

Equations  (3.6)  and  (3.7)  were  used  in  determining  the 

storage  coefficient  S  ,  resulting  in  a  value  of  1.2  x 

ID'3. 

A  regional  technique,  Equation  (3.10),  was  used  in 
the  determination  of  average  aquifer  transmissivity  T , 
for  combined  discharge  Q_ ,  of  the  freshwater  springs,  and 
the  mixed-water  springs  (Table  20) .  Estimates  of  average 
spring  discharge  coincided  with  the  period  of  the  poten- 
tiometric  survey  (25  January  to  0  8  February,  198  5) . 
Also,  the  average  hydraulic  gradient  I,  and  the  average 
contributing  aquifer  width  W,  were  extracted  from  the 
potentiometric  survey. 

Average  transmissivity  of  the  unconfined  Arbuckle 

Aquifer  responsible  for  combined  discharge  of  the  fresh- 

3    2 
water  springs  is  about  43  x  10   ft  /day.   Also,  the 

average   transmissivity   of   the   confined   Simpson   and 

Arbuckle  Aquifers  responsible  for  the  combined  discharge 

2 
of   mixed-water   springs   is   approximately   3.9   x   10 

ft2 /day. 
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4.4   Wells  and  Springs 

East  and  West  Observation  Wells,  completed  in  the 
Arbuckle  Group,  are  located  in  close  proximity  to  the 
artesian,  freshwater  springs  (Plate  III)  .  The  original 
intent  of  these  wells  were  to  monitor  deprivation  effects 
of  adjacent  production  wells  outside  the  confines  of  the 
park  and  establish  a  relationship  of  spring  discharge 
versus  piezometric  level  in  the  wells  (Moore,  1972). 
Depth  to  water  table  measurements  for  the  beginning  of 
each  month  (since  January  1974)  for  each  observation  well 
are  listed  in  Tables  21  and  22. 

Monthly-average  stage  data  for  Buffalo  and  Antelope 
are  listed  in  Tables  23  and  24,  respectively.  As  de- 
termined by  previous  investigators,  flow  of  both  springs 
can  be  approximated  by  modified  weir  equations.  Monthly 
flow  averages  are  calculated  from  stage  data  (in  Tables 
23  and  24)  and  are  listed  for  both  springs  in  Tables  25 
and  26. 

Hydrographs  of  East  and  West  Observation  Wells,  and 
of  Buffalo  and  Antelope  Springs  are  depicted  in  Figure 
42.  Monthly  precipitation  histogram,  long-term  trend  of 
hydrographs,  and  no-flow  periods  of  Buffalo  and  Antelope 
are  also  included.  Figure  42  graphically  demonstrates 
spring  discharge  as  a  function  of  precipitation  events 
and  water  table  elevations  of  the  East  Observation  Well. 

The  slope  of  the  long-term  trend  of  hydrograph  for 
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Table  21.   Ground  Water  Levels  East  Observation  Well 

Chickasaw  National  Recreation  Area 
Sulphur,  Okla. 
(Approximate  Surface  Elevation  =  1150') 


Depth  to  water,  below  land  surface  datum 
(Beginning  of  Month  Readings) 

Year  Jan    Feb   Mar   Apr  May  June  July  Aug  Sept  Oct  Nov  Dec 

1974*  38.5    44.5   47.0   49.5  51.0  51.0  54.0  58.0  61.0  59.0  56.0  44.0 

1975*  49.0    51.0   42.0   39.0  41.0  40.5  44.5  50.0  56.0  59.5  62.0  64.0 

1976*  66.0   67.5   69.5   68.0  67.0  64.5  65.0  63.5  69.5  71.0  71.0  73.0 

1977*  74.5    76.5   75.5   67.0  64.0+  63.0  64.5  66.0  57.0  59.0  60.5  61.5 

1978*  63.0    66.0   76.0   73.5  71.0  67.5  61.0  64.5  66.0  68.5  71.0  73.0 

1979*  75.0    76.0   77.0   71.0  67.0  63.5  62.0  63.5  66.5  68.5  70.0  71.5 

1980*    —    —    76.5   76.5  81.0  82.5  81.5  83.5  87.0  86.0  86.0  87.0 

1981  86.5*   88.0*  89.0*  86.5*  86.5*  85.5*  75.5*   76.8  78.1  80.2  59.2  58.1 

1982  61.6   59.8    57.2+  58.7  60.5  48.0+  44.2  47.9  53.2  .  57.0  61.2  64.8 

1983  66.8   66.7+   62.0   61.5  60.8  45.2  49.7  56.1  60.9  63.1  65.5  66.3 

1984  70.0   70.9    71.9+  70.0  68.9  70.5  71.4  73.3  75.8  77.7  74.2  69.5 


NOTE:   Flow  at  both  Antelope  and  Buffalo  Springs  cease 
when  this  well  reaches  a  depth  of  between  68  and 
72  feet  below  the  ground  surface. 

*  Data  taken  from  graph  -  original  data  unavailable 

+  Denotes  questionable  value 


(after  Harp,  1984) 
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Table  22.   Ground  Water  Levels  West  Observation  Well 

Chickasaw  National  Recreation  Area 
Sulphur,  Okla. 
(Approximate  Surface  Elevation  =  1080') 


Depth  to  water,  below  land  surface  datum 
(Beginning  of  Month  Readings) 


Year  Jan    Feb   Mar   Apr   May    June   July  Aug  Sept  Oct  Nov  Dec 

1974*  23.0  24.5   23.5   24.0  23.0  25.0  26.5  27.0  28.5  25.5  23.0  23.0 

1975*  24.5  23.5   23.0   20.5  21.0  20.5  23.5  23.5  26.0  27.0  27.0  27.0 

1976*  26.5  27.5   29.0   28.0  28.0  26.5  29.0  29.0  3D. 5  30.5  31.0  31.0 

1977*  31.5  31.0   30.0   29.0  29.0  28.0  28.5  29.5  30.5  30.5  30.5  30.5 

1978*  30.5  30.5   29.5   29.5  29.5  29.0  29.0  29.0  29.0  29.0  29.0  29.0 

1979*  29.5  29.5   29.0   29.0  28.5  25.0  26.0  31.0  31.0  31.5  32.5  33.0 

1980*  —  —    32.0   32.0  33.0  33.0  33.0  34.5  34.5  33.0  34.0  33.5 

1981  32.0*  32.5*  32.5*  31.5*  32.5*  30.8  31.8  33.0  33.0  33.5  28.8  28.6 

1982  32.2  28.5   28.3+  29.4  30.2  26.3  27.3  29.1  29.8  29.0  29.5  30.5 

1983  --  28.7+  28.3   28.1  28.6  27.2  28.1  29.2  30.3  30.6  29.3  29.5 

1984  29.2  29.2   30.1   29.0  28.8  30.9  31.4  32.1  32.4  32.7  30.2  29.6 


*  Data  taken  from  graph  -  original  data  unavailable 
+  Denotes  questionable  data 
(after  Harp,  1984) 
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Table  23.   Stage  Data  Buffalo  Springs 

Chickasaw  National  Recreation  Area 
Sulphur,  Okla. 
(Approximate  Surface  Elevation  =  1078') 


Year   Jan    Feb   Mar   Apr   May    June   July   Aug 


Sept   Oct   Nov   Dec 


1972  -      -  1.95  1.91  1.99 

1973  2.08  2.10  2.15  2.18  2.16  2.13  2.10  2.08 

1974  2.15  2.12  2.08  2.08  2.08  2.09  2.02  2.01 

1975  2.10  2.13  2.15  2.13  2.13  2.10  1.00  2.05 

1976  1.97  1.94  1.96  1.96  1.99  1.98  1.94  1.90 

1977  Dry  Dry  Dry  1.99  2.00  1.99  1.97  1.95 

1978  0.46+  Dry  Dry  1.79  1.93  2.01  2.01  1.95' 

1979  Dry  Dry  1.93  1.97  2.00  2.01  1.02  1.07 

1980  Dry  Dry  Dry  Dry  Dry  Dry  Dry  Dry 

1981  Dry  Dry  Dry  Dry  Dry  Dry  Dry  Dry 

1982  2.23  2.22+  2.22+  2.22  2.26+  2.26  2.21  2.21 

1983  2.16  2.29  2.20  2.20  2.28  2.23  2.21  2.20 

1984  1.78  1.54  1.66  2.06  2.03  1.83  1.17  .05 


1.49  1.06  2.00  - 

2.08  2.12  2.18  2.17 

2.03  2.07  2.13  2.11 

2.05  2.01  1.99  2.00 

1.86  1.56  1.33  0.52 

1.92  1.82  1.54  1.00+ 

1.90  1.76  1.36  Dry 

1.92  1.90  1.59  0.93 

Dry  Dry  Dry  Dry 

Dry  2.20  2.25  2.22 

2.15+  2.13+  2.14+ 

2.16+  2.70  2.10  1.49 

Dry  Dry  1.73  2.18 


+   Denotes  questionable  data 

Monthly  average  as  derived  from  field  data  collected 
during  the  month 


(after  Harp,  1984) 


166 


Table   24.   Stage  Data  Antelope  Springs 

Chickasaw  National  Recreation  Area 
Sulphur,  Okla. 
(Approximate  Surface  Elevation  =  1080') 


Year   Jan    Feb   Mar   Apr   May    June   July   Aug   Sept   Oct   Nov   Dec 


1969  -                 -     -  0.95 

1970  0.87  0.86  0.89  0.92  0.93  0.86 

1971  0.91  0.89  0.87  0.83  0.80+  0.75 

1972  0.88  0.87  0.85+  0.84  0.84  0.78 

1973  0.90  0.97  1.05  1.08  1.07  1.05 

1974  1.06  1.05  1.05  1.06  1.03+  1.01 

1975  1.05  1.11  1.20  1.23  1.23  1.19 

1976  0.76  0.84  0.86  0.86  0.96  0.94 

1977  Dry  Dry  0.84  0.92  0.91  0.90 

1978  Dry  Dry  Dry  0.71  0.80  0.88 

1979  Dry  Dry  0.78  0.78  0.88  0.90 

1980  Dry  Dry  Dry  Dry  Dry  Dry 

1981  Dry  Dry  Dry  Dry  Dry  Dry 

1982  0.94+  0.96+  0.98+   -  0.94  1.02 

1983  0.79  0.87  0.90  0.88  1.00  0.94 

1984  0.65         .50  .56        .73           .69  .63 


0.92  0.88  0.87  0.88  0.88  0.89+ 

0.83  0.79  0.82  0.96  0.94  0.94 

0.69  0.70  0.64  0.70  0.70  0.87 

0.73  0.68  0.62+  0.54  0.87  0.89 

1.01  0.99  1.01  1.02  1.11  1.08 

0.96  0.91  0.92+  1.04  1.03  1.02 

1.00  0.98  0.95  0.95  0.90  0.88 

0.84  0.80  0.73  0.65  Dry  Dry 

0.87  0.83  0.80  0.77  0.70  Dry 

0.84  0.81  0.76  0.68  0.65  Dry 

0.84  0.80  0.75  0.74  0.68+  Dry 

Dry  Dry  Dry  Dry  Dry  Dry 

Dry  Dry  Dry  0.95  0.99  0.99+ 

0.93  1.91+  1.74+  1.29+  0.83+  0.83+ 

0.92  0.84+  0.87+  0.81  0.79  0.78 

.12  Dry  Dry  Dry  .47  .83 


+  Denotes  approximate  value 

Monthly  average  as  derived  from  field  data  collected 
during  the  month 


(after  Harp,  1984) 
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Table  25.   Spring  Flow  in  Cubic  Feet  Per  Second 
Buffalo  Springs 

Chickasaw  National  Recreation  Area 
Sulphur,  Okla. 

For  the  Years  1972  -  1983 


Year 


Jan 


Feb 


Mar   Apr   May 


June 


July   Aug   Sept   Oct 


Nov 


Dec 


1972 

- 

1973 

1.38 

1974 

1.98 

1975 

1.55 

1976 

0.61 

1977 

0.00 

1978 

0.00 

1979 

0.00 

1980 

0.00 

1981 

0.00 

1982 

1.38 

1983 

0.86 

1984 

0.00 

1.55  1.98  2.26  2.07 

1.72  1.38  1.38  1.38 

1.80  1.98  1.80  1.80 

0.43  0.55  0.55  0.73 

0.00  0.00  0.73  0.80 

0.00  0.00  0.00  0.38 

0.00  0.38  0.61  0.80 

0.00  0.00  0.00  0.00 

0.00  0.00  0.00  0.00 

1.30  1.30  1.30  1.63 

1.89  1.15  1.15  1.80 

0.00  0.00  .28  .15 


0.49 

0.28 

1.80 

1.55 

1.46 

0.93 

1.55 

0.73 

0.67 

0.43 

0.73 

0.61 

0.86 

0.86 

0.86 

0.33 

0.00 

0.00 

0.00 

0.00 

1.63 

1.23 

1.38 

1.23 

0.00 

0.00 

0.73  0.00  0.00  0.80  0.00 

1.38  1.38  1.72  2.26  2.16 

0.86  1.00  1.31  1.80  1.63 

1.15  1.15  0.86  0.73  0.80 

0.24  0.08  0.00  0.00  0.00 

0.49  0.33  0.00  0.00  0.00 

0.49  0.24  0.00  0.00  0.00 

0.55  0.33  0.24  0.00  0.00 

0.00  0.00  0.00  0.00  0.00 

0.00  0.00  1.15+  1.55  1.30 

1.23  1.02  0.80  0.67  0.73 

1.15  0.86  .65  0.49  0.0 

0.00  0.00  0.00  0.00  1.00 


All  flows  calculated  using  the  equation:   (for  data  thru 
October  of  1981) 


Q  =  10.44(P.R.-1.82) 


1.5 


After  October  of  1981,  the  equation  was  changed  to: 

Q  =  10.44(P.R.-1.97)1*5 

where:   Q  =  flow  in  cfs, 

P.R.  =  pin  reading  (stage  of  spring,  in  feet 

+  Pin  reset,  approx.  0.15'  lower  than  original  setting 

Monthly  mean  flow  as  derived  from  data  of  Table  22. 

(after  Harp,  1984) 
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Table  26.   Spring  Flow  in  Cubic  Feet  Per  Second 
Antelope  Springs 

Chickasaw  National  Recreation  Area 
Sulphur,  Okla. 

For  the  Years  1969  -  1983 


Year   Jan  Feb   Mar   Apr   May    June  July   Aug  Sept  Oct  Nov  Dec 

1969  -  -                      1.80  1.57    1.29  1.22  1.29  1.29  1.36 

1970  1.22  1.15   1.36   1.57   1.65    1.15  0.96   0.72  0.90  1.88  1.72  1.72 

1971  1.50  1.36   1.22   0.96   0.78   0.50  0.23    0.27  0.07  0.27  0.27  1.22 

1972  1.29  1.22   1.09   1.02   1.02   0.66  0.41   0.20  0.02  0.00  1.22  1.36 

1973  1.43  1.96   2.62   2.89   2.80    2.62  2.28    2.12  2.28  2.37  3.17  2.89 

1974  2.71  2.62   2.62   2.71   2.45    2.28  1.88    1.50  1.50  2.54  2.45  2.37 

1975  2.62  3.17   4.04   4.35   4.35    3.94  2.20    2.04  1.80  1.80  1.43  1.29 

1976  0.56  1.02   1.15   1.15   1.88    1.72  1.02   0.78  0.41  0.10  0.00  0.00 

1977  0.00  0.00   1.02   1.57   1.50    1.43  1.22   0.96  0.78  0.61  0.27  0.00 

1978  0.00  0.00   0.00   0.32   0.78    1.29  1.02   0.84  0.56  0.20  0.10  0.00 

1979  0.00  0.00   0.66   0.66   1.29    1.43  1.02    0.89  0.50  0.46  0.20  0.00 

1980  0.00  0.00   0.00   0.00   0.00   0.00  0.00    0.00  0.00  0.00  0.00  0.00 

1981  0.00  0.00   0.00   0.00   0.00    0.00  0.00   0.00  0.00  1.80  2.12  2.12 

1982  1.72  1.88   2.04   1.88   1.72    2.37  1.65    1.30  1.06  1.01  0.96  0.96 

1983  0.72  1.22   1.43   1.29   2.20    1.72  1.57    1.02  1.22  0.84  0.72  0.66 

1984  0.10  0.00   0.00   0.41   0.24   0.05  0.00   0.00  0.00  0.00  0.00  0.96 

All  Flows  calculated  using  the  equation: 

Q  =  8.69(P.R.-0.6)  * 

where  Q  =  flow  in  cfs 

P.R.  =  pin  reading  (stage  of  spring,  in  feet) 

Monthly  mean  flow  as  derived  from  data  of  Table  23. 
(after  Harp,  1984) 


169 


w 

C         <U    1 

O       Si  u 

•H  T3            3 

4JCCM 

rd  fd  fd 

>         US 

d)  -P         O 

H    (1)   C   M 

d)  a;  o  <w 

IW  -H 

(D      -p  a) 

uo  n!H 

fd  in  >  ,Q 

iw  rH   0)  rd 

H  rH  rH   -P 

3       a) 

10    d>         5-1 

5-1    U    d> 

rd  -H  -P 

..       u  rd 

d)  en  -p  £ 

+J    rH     Q) 

O   rH     S     O 

zoo+J 

12   N 

<u  x: 

e  -h  -p 

•    O  P*    04 

x:  -h      <u 

CU-P        'O 

fd  fd    • 

}-i  >  >i  tr> 

CJl  5-4  rH    c 

o  Q)  <u  -h 

u    W    >  -P 

«d  .a  -h  o 

>iO  -P  fd 

,c        OH 

-p  d>  -p 

cr>  en  ChU 

• 

fi    d)    W    £3 

C 

•h  IS  <l>  tn 

0 

H           V4 

•H 

cu^      >iP 

w  c    -ja 

id 

id  -p 

> 

T3         <D  ti 

0) 

c  -p  a)  a) 

rH 

d  to<n  c 

a) 

td      -h 

rH    W   O     fd 

0) 

rH            00    ^J 

u 

Q)  <H  O  X5 

rd 

^     O   rH     O 

M-l 

• 

CN 

■^ 

a) 

h 

3 

cp 

•H 

Pm 

170 

East  and  West  Observation  wells  is  increasing  (piezo- 
metric  level  is  decreasing)  at  a  rate  of  .65  ft/year 
(Figure  42)  .  The  depth  to  water  table  is  approximately 
72  feet  in  the  East  Observation  Well  when  springs  com- 
mence or  cease  to  flow  (Figure  42)  .  The  following  ex- 
pression can  be  used  to  forecast  total  cessation  of 
spring  flow  with  reasonable  certainty: 
Y    -  Y 

i-JE S  ,4.i) 

bh 

where   X  =  number  of  years  when  springs  cease  flowing 
from  some  time  T, 

Y  =  72  feet  =  depth  to  water  table  in  an  observa- 

tion well  when  springs  cease  flowing  (feet) , 

Y  =  Y.qft.  =  45  feet  =  depth  to  water  table  ob- 

tained from  long-term  trend  of  hydrograph  of 

an  observation  well  at  some  time  T  (feet), 
b,   =  .65  feet/year  =  deprivation  slope  as 

determined  from  long-term  trend  of  hydrograph 

(ft/year) . 
Thus,  in  approximately  4  0  years  from  the  submittal  of 
this  report,  a  strong  possibility  exists  that  the  springs 
will  cease  flowing  unless  certain  conservation  measures 
are  undertaken.  The  long-term  trend  of  hydrograph  can  be 
attributed  to  deprivational  effects  of  adjacent  sinks  and 
not  to  a  lack  of  precipitation  for  the  time  period  1974 
to  1984  (Figure  42) . 
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Peak  discharge,  minimum  depth  to  water  table,  and 
maximum  precipitation  during  the  calendar  year  occur  in 
either  middle  spring  or  early  fall  (Figure  42)  .  Con- 
versely, minimum  discharge,  maximum  depth  to  water  table, 
and  minimum  precipitation  generally  occur  in  late  summer 
or  early  winter. 

Spring  hydrographs  of  Antelope  and  Buffalo  demon- 
strate that  both  eminate  from  the  same  aquifer.  In  addi- 
tion, the  East  Observation  Well,  1,000  feet  east- 
northeast  of  Antelope  and  1,000  feet  north-northeast  of 
Buffalo,  fluctuates  in  response  to  precipitation  events 
in  a  very  similar  manner  to  that  of  wells  located  in  the 
Arbuckle  Group  outcrop,  1-mile  east  of  the  park  boundary. 
The  piezometric  level  of  Crawford's  Well,  located  in  the 
outcrop  area,  rose  nearly  38.5  feet  for  the  time  period, 
15  October  1984  to  25  January  1985  (Table  27) .  Likewise, 
the  piezometric  level  of  the  East  Observation  Well  rose 
23.3  feet  for  the  same  period  of  time. 

Figure  43  depicts  periods  of  no  flow  of  Antelope  and 
Buffalo  springs  with  respect  to  annual  precipitation 
amounts,  from  1919  to  1984.  Particularly  noteworthy,  for 
the  first  time  in  66  years  of  recorded  observations,  the 
cessation  of  spring  flow  occurred  in  1984  when  3  years  of 
above-average  annual  precipitation  preceded  it  (Figure 
43)  .  This  anomalous  behavior  is  indicative  of  long-term 
deprivation  effects  of  free  flowing  wells  and  production 
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Table  27.   Piezometric  Levels  of  Selected  Wells 
in  the  Arbuckle  Group 


East  Observation        +Crawford's 
Well  Well 

Piezometric  Elevation  Piezometric  Elevation 
Date                 (feet)  (feet) 


10/15/84  1071.7  1073.0 

10/27/84  1079.0  1082.0 

11/10/84  1076.3  1085.7 

11/16/84  —  1083.8 

11/19/84  1078.9 

11/21/84  —  1089.2 

1/25/85  1095.0  1111.5 


+  Location  of  Crawford's  Well  -  C  NEj  SEj  of  Sec.  31, 
TIN,  R4E. 
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wells  that  are  gradually  lowering  the  piezometric  level 
in  the  Sulphur  artesian  basin. 

The  inch-year,  time  bar  chart,  Figure  44,  portrays 
the  summation  of  rainfall  that  falls  above  (positive)  or 
below  (negative)  the  mean  annual  precipitation  of  37.96 
inches  (1919-1984)  for  a  specified  period  (i.e.,  15  or  20 
years)  .  Inside  each  bar  is  a  number  that  indicates  the 
total  no-flow  time  of  the  freshwater  springs.  Depriva- 
tion caused  by  withdrawal  of  water  can  be  firmly  estab- 
lished by  the  use  of  this  bar  chart.  For  example,  a  to- 
tal no-flow  period  of  .58  years  occurred  during  the 
15-year  span  from  1919  to  1934.  The  total  amount  of  inch 
years  is  +12.0  (Figure  44).  However,  for  the  time  peri- 
od, 1964-1984,  the  no-flow  period  is  5.08  years  and  the 
inch  years  is  a  +30.0.  Comparison  of  the  two  time  in- 
tervals reveals  that  the  total  amount  of  cessation  is  ap- 
proximately 10  times  longer  during  the  time  frame  1964- 
1984;  even  though, the  total  inch  years  of  1964  to  1984 
exceeded  the  total  inch  years  of  1919  to  1934. 

Locations  of  flowing  wells,  observation  wells,  and 
springs  are  shown  on  Plate  III.  Pumpage  data  of  major 
production  wells  in  close  vicinity  to  the  park  is  listed 
in  Table  28.  The  investigator  was  unable  to  obtain 
pumpage  data  from  O.G.E.  for  the  time  periods:  1973  to 
1979,  and  1983  to  1984  (Table  27). 

The  major  artesian  well  in  close  proximity  to  the 
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Figure    44. 


Inch-year,  time  bar  chart  that  depicts 

above-  or  below-normal  precipitation 
amounts  for  a  specified  period,  in  re- 
lation to  no-flow  periods  of  Antelope 
and  Buffalo. 
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Table  28.   Well  Pumpage  Data, Area  Well  Fields 

Chickasaw  National  Recreation  Area 
Sulphur,  Okla. 


O.G.  &  E 

City  of 

Year 

Wells 

Sulphur 

1964 

— 

— 

1965 

348.66 

240.14 

1966 

296.52 

364.95 

1967 

291.31 

300.11 

1968 

302.39 

210.83 

1969 

328.46 

211.80 

1970 

440.55 

228.10 

1971 

334.98 

235.26 

1972 

348.66 

271.43 

1973 

— 

361.56 

1974 

/ 

340.45 

1975 

— 

274.71 

1976 

— 

298.73 

1977 

— 

336.00 

1978 

— 

— 

1979 

— 

438.20 

1980 

472.39 

371.47 

1981 

466.97 

224.82 

1982 

480.88 

400.75 

1983 

— 

376.02 

1984 

— 

392.01 

Roos 
Ranch 


169.44 
102.97 


162.93 
173.68 
173.68 
217.34 


(All  data  in  million  gallons) 
(after  Harp,  1984) 
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park,  the  Vendome  Well,  was  estimated  to  flow  at  approxi- 
mately  2  ft  /sec  or  472  x  10  gallons  per  year  (Appendix 
IV) .   Currently,  the  total  combined  amount  of  water  that 

is  extracted  by  O.G.E.,  Sulphur  Well  Field,  and  Vendome 

9 
Well  on  an  annual  basis,  is  approximately  1.4  x  10   gal- 
lons.  The  total  combined  discharge  of  the  remainder  of 
the  flowing  wells  in  the  Sulphur  artesian  basin  is  un- 
known . 

The  major  findings  of  the  two  downhole  TV  surveys 
conducted  during  the  calendar  year  1980  on  selected  wells 
in  the  project  area  are  as  follows: 

(1)  A  downhole  TV  survey  was  completed  on  20  August 
1980  of  the  Vendome  Well,  located  at  NEj  NEj 
NWJ  NEJ  of  Sec.  3  T1S,  R3E;  surface  elevation, 
945  feet.  The  Vendome  Well  is  cased  to  a  depth 
of  300  feet  below  the  surface,  with  an  inner 
diameter  of  .5  feet.  The  overall  condition  of 
the  casing  is  fair,  but  it  is  extremely  cor- 
roded at  depths  between  250  and  300  feet.  A 
fairly  uniform  conglomerate  formation  consist- 
ing of  boulders  and  cobbles  is  exposed  in  the 
open  well  from  300  to  320  feet.  A  large 
cavity,  3  feet  in  thickness,  was  observed 
between  320  and  323  feet.  A  distinctly  dif- 
ferent formation,  massive  in  appearance, 
appears  between  a  depth  of  323  and  324  feet. 
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The  total  depth  of  the  well  is  324  feet.  The 
original  depth  of  the  well,  completed  in  1922, 
was  364  feet. 
(2)  A  downhole  TV  survey  was  conducted  on  13  and 
14  March  1980,  of  the  Horseman  Well,  located 
NEi  NWi  NEi  NWj  of  Sec.  1,  T1S,  R3E ;  approxi- 
mate surface  elevation,  1,075  feet.  The  condi- 
tion of  casing  is  good  and  extends  to  a  depth 
of  101  feet,  with  a  .5  feet  inner  diameter. 
Open  well  bore  extends  from  101  feet  to  632 
feet  below  the  surface.  Water  was  extremely 
murky  to  a  depth  of  380  feet,  and  as  a  result, 
visibility  was  poor.  The  open  well  bore  pene- 
trated only  the  West  Spring  Formation  of  the 
Arbuckle  Group.  Near-vertical  fractures,  solu- 
tion cavities,  and  vugs  were  observed  at  all 
depths  in  the  open  well  bore.  At  times,  the 
well  bore  periphery  resembled  a  swiss  cheese 
texture.  Less  frequently,  smooth  sections  of 
5-20  feet  thick  of  the  open  well  bore  were 
observed.  This  would  indicate  a  change  in 
lithology  within  the  West  Spring  Formation, 
Arbuckle  Group. 
Mean  annual  flow  (cfs)  and  one  standard  deviation  is 
tabulated  for  both  freshwater  springs  in  Table  29.  In- 
cluded in  Table  29  is  the  annual  rainfall  amount.   Also, 
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Table  29.   Mean  Annual  Flow  of  Antelope  and  Buffalo 
Springs,  and  Annual  Precipitation 


Year 

Antelope  Spring 
Mean  Flow  and 
*S.D.  (cfs) 

Buffalo  Spring 
Mean  Flow  and 
S.D.  (cfs) 

Annual 
Precipitation 
(inches) 

1970 

1.33 

+ 

.37 

— 

50.07 

1971 

.72 

+ 

.512 

— 

36.35 

1972 

.79 

+ 

.51 

— 

31.85 

1973 

2.45 

+ 

.48 

1.79 

+ 

.35 

58.93 

1974 

2.3 

+ 

.44 

1.40 

+ 

.35 

43.48 

1975 

2.75 

+ 

1.16 

1.32 

+ 

.47 

38.33 

1976 

.81 

+ 

.63 

.36 

+ 

.28 

28.67 

1977 

.78 

+ 

.6 

.31 

+ 

.34 

31.09 

1978 

.43 

+ 

.46 

.24 

+ 

.34 

34.91 

1979 

.58 

+ 

.49 

.34 

+ 

.31 

35.43 

1980 

0.00 

0.00 

32.98 

1981 

.50 

+ 

.91 

.33 

+ 

.61 

45.77 

1982 

1.54 

+ 

.48 

1.19 

+ 

.32 

41.65 

1983 

1.21 

+ 

.47 

1.05 

+ 

.53 

40.48 

1984 

.15 

+ 

.29 

.12 

+ 

.29 

42.59 

*  S.D.  -  one  standard  deviation 
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mean  precipitation  and  mean  flow  of  the  springs  on  a 

monthly  basis  are  listed  in  Table  30.   Data  entered  in 

Table  30  demonstrates  that  the  maximum  average  discharge 

of  the  springs  occurs  in  May  and  the  minimum  discharge 

occurs  in  August,  for  the  time  period  of  1969  to  1984. 

Spring  variance  V*    (Equation  3.13),  for  Antelope 

Spring  is  430%  for  the  time  period  of  June  1969  to 

December  1984  (Table  31) .   Also,  the  spring  variance  V  , 

a 

for  Buffalo  Spring  is  390%  for  the  time  period  of  June 

1972  to  December  1984  (Table  31) .   The  values  of  spring 

variance  V    for  both  springs  are  fairly  close.   The 
a 

relatively  high  variance  values  indicate  that  the 
Arbuckle  Aquifer  has  a  low  storage  capacity. 

Exhaustion  coefficients  e,  or  diffusivity  constants 
K,  (are  used  interchangeably) ,  and  response  times  t  ,  are 
listed  in  Table  32.  Recession  limbs  selected  for  analy- 
sis are  from  calendar  years  1976  to  1979.  The  exhaustion 
coefficients  ranged  from  .49  to  .60  mos-  and  the 
response  times  varied  from  50  to  61  days. 

Finally,  on-site  reconnoiter  of  Buffalo  Spring,  when 
flow  resumed,  revealed  that  there  is  additional  water 
welling  up  in  the  creek  bed  in  close  proximity  to  Buffalo 
Spring.  Total  discharge  of  water  is  unknown  at  this  vi- 
cinity, but  shouldn't  exceed  the  flow  at  Antelope  by  more 
than  50%.  This  was  based  on  visual  comparison  of  the 
streams  that  flow  from  the  spring  areas. 
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Table  30.   Mean  Precipitation  and  Mean  Flow 
of  the  Springs  on  a  Monthly  Basis 


(1969-1984) 
Mean  Precipitation 

(June  1969-1984) 
Antelope  Spring 

(June  1972-1984 
Buffalo  Spring 

Month 

* 
and   S.D. 

(in.) 

Mean  Flow  and  S.D. 
(cfs) 

Mean  Flow  and  S.D. 
(cfs) 

Jan. 

1.38  ± 

.89 

.92  ± 

.95 

.65 

±  .76 

Feb. 

2.07  ± 

1.10 

1.04  ± 

1.04 

.72 

±  .84 

March 

3.16  ± 

1.85 

1.28  ± 

1.12 

.73 

±  .79 

April 

3.48  ± 

1.45 

1.39  ± 

1.08 

.84 

±  .74 

May 

6.29  ± 

3.65 

.  1.51  ± 

1.16 

.96 

±  .75 

June 

3.83  ± 

1.83 

1.44  ± 

1.07 

.88 

±  .65 

July 

2.38  ± 

1.33 

1.06  ± 

.77 

.64 

±  .49 

August 

2.29  ± 

1.49 

.87  ± 

.67 

.64 

±  .49 

Sept. 

4.29  ± 

3.80 

.77  ± 

.70 

.49 

±  .51 

October   5.98  ± 

4.68 

.95  ± 

.89 

.52 

±  .60 

Nov. 

2.50  ± 

1.79 

1.00  ± 

.98 

.64 

±  .78 

Dec. 

1.79  ± 

1.59 

1.06  ± 

.54 

.59 

±  .75 

S.D.  =  one  standard  deviation 
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Table  31.   Input  Variables  and  Resultant  Values  of  Spring 
Variance  for  Antelope  and  Buffalo  Springs 


Maximum 

Minimum 

*Median 

+Spring 

Time 

Discharge 

Discharge 

Discharge 

Variance 

Spring  Period 

Q   (cfs) 
max 

Q  .  (cfs) 

nun 

Wcfs) 

V  % 
a 

June  1972 

Buffalo   to 

2.26 

0. 

.575 

390 

Dec.  1984 

June  1969 

Antelope  to 

4.35 

0. 

1.013 

430 

Dec.  1984 

*  Median  discharge  was  determined  by  statistical 
analysis,  interval  is  .25  cfs. 


Q     -  Q  • 
+  Spring  Variance  =  V  =  -^^ 2LH1  (100%) 

3. 


Q 


med 
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Table    32.      Data    Input   and   Resultant  Values   of 
Exhaustion  Coefficient   and   Response 
Time   of  Antelope   Spring 

Calendar     Base  Flow     Inflection  Point     Time        *Exhaustion     +Response 

Year  Q    (cfs)  Flow  Q      (cfs)  t(mos)      Coefficient     Time 

o 


1976  .1 

1977  .1 

1978  .1 

1979  .1 


In  (Q  /Q) 

*  Exhaustion  coefficient  =  e  =  K,  =  

d      t 
30 
+  Response  time  =  t   =  — 

e 


e    (mos 

*> 

t 

r 

(day) 

1.0 

4.0 

.60 

50 

1.5 

5.5 

.49 

61 

1.0 

4.0 

.57 

53 

1.0 

4.5 

.51 

59 

184 


4.5   Recharge  Areas  Responsible  for  Spring  Flow  in  the 
Travertine  District,  Chickasaw  National  Recreation 
Area 

As  stated  previously,  a  potentiometric  survey  of  the 
Arbuckle  Mountain  area  was  completed  by  Fairchild  et  al. 
(1983),  in  the  winter  of  1976-1977  (Plate  II).   This  map 
was  used  to  define  the  project  area. 

The  potentiometric  survey  for  this  investigation  was 
conducted  from  25  January  to  8  February  1985.  Results  of 
the  survey  are  listed  in  Table  33.  Depth  to  water  table 
readings  were  rounded  off  to  the  nearest  foot,  and  sur- 
face elevations  of  wells  were  estimated  from  topographic 
maps  (Table  33) . 

A  potentiometric  map,  Figure  45,  was  constructed 
from  piezometric  elevations  of  the  Simpson  and  Arbuckle 
Groups.  Also,  equipotential  lines  were  inferred  from 
surface  elevations  of  flowing  wells  in  the  city  of  Sul- 
phur and  surface  elevations  of  springs  in  the  central  and 
western  portion  of  the  Travertine  District,  Chickasaw  Na- 
tional Recreation  Area. 

Recharge  areas  were  primarily  determined  from  the 
potentiometric  survey  conducted  by  Fairchild  et  al. 
(Winter,  1976-1977) ,  and  the  survey  for  this  study.  Geo- 
logic control,  primarily  fractures  and  faults,  were  con- 
sidered also.  The  areas  of  recharge  are  shown  on  the 
geologic  and  hydrogeologic  map,  Figure  30. 

The  recharge  area  for  Antelope  and  Buffalo  is  ap- 
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Figure  45.   Potentiometric  map  (Jan.  and  Feb.,  1985) 
of  the  Simpson -Arbuckle  Aquifer,  Sulphur, 
Oklahoma,  area. 
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2 
proximately  5  mi.  ,  located  north  of  the  west-northwest 

trending  Sulphur  Fault  (Figure  30) .   The  primary  recharge 

area  for  the  remainder  of  the  springs  in  the  park  is  ap- 

2 
proximately  2  mi .  ,  located  south  of  the  Sulphur  Fault, 

along  the  periphery  of  the  Sulphur  Syncline  (Figure  30) . 


4.6   Concise  Nomenclature  and  Origin  of  the  Fresh  Water 
Springs,  Buffalo  and  Antelope. 

Buffalo  and  Antelope  are  third-magnitude,  inter- 
mediate-response, freshwater  artesian  springs;  that  de- 
rive their  waters  from  the  Arbuckle  Group. 

The  springs  occur  as  a  result  of  fractures  that  be- 
come solution  cavities  in  the  confining  layer  of  the 
Pontotoc  that  extend  to  the  surface  top  of  the  pres- 
surized water  of  the  Arbuckle  Group,  just  north  of  the 
Sulphur  Fault.  Both  springs  are  topographically  low, 
compared  to  the  surrounding  area.  If  the  Arbuckle  Group 
receives  sufficient  recharge  in  the  outcrop  area,  the 
piezometric  level  will  rise  above  the  surface  elevation 
of  the  springs,  and  discharge  will  result. 

Groundwater  flow  to  the  springs  is  controlled  by  the 
Sulphur  Fault.  The  elevated  Arbuckle  dolomite  on  the  up- 
thrown  side  of  the  fault  is  extensively  fractured  and 
contains  solution  cavities  and  small  caverns. 
Third-magnitude  springs  normally  eminate  from  a  highly 
permeable  aquifer.  Low  hydraulic  gradients  (10  ft.  per 
mile)  in  the  outcrop  area  of  the  Arbuckle  Group,  east  of 
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Sulphur,  indicate  high  regional  transmissivity.  Con- 
versely, high  hydraulic  gradients  (20  to  100  ft.  per 
mile)  of  the  Simpson  Group,  south  of  the  Sulphur  Fault, 
indicate  lower  transmissivities. 

4.7   Hydrologic  Budget  of  the  Sulphur  Area 

The  piezometric  level  of  the  Simpson-Arbuckle 
aquifer  should  remain  constant  if  inflow  equals  outflow. 
However,  the  piezometric  level  of  the  Sulphur  artesian 
basin  is  declining  at  an  average  annual  rate  of  .65  ft 
per  year.   Hence,  outflow  exceeds  inflow  for  the  area. 

As  previously  discussed,  the  total  amount  of  water 

that  is  currently  withdrawn  from  the  Sulphur  area,  is  in 

9 
excess  of  1.4  x  10   gallons  per  year.   The  total  recharge 

area  for  the  Sulphur  artesian  basin  is  approximately  15 

2 
mi   (Figure  30)  .   Assuming  a  long-term  annual  recharge 

rate  of  4.6  inches  per  year  (Fairchild  et  al.,  1983),  the 

total  amount  of  inflow  on  a  yearly  basis  is  approximately 

9 
1.2  x  10   gallons  per  year.   Since  outflow  is  equivalent 

to  withdrawal,  and  withdrawal  exceeds  outflow;  it  becomes 

quite   evident   as   to   why   the   piezometric   level   is 

declining  in  the  Sulphur  area. 


CHAPTER  V 


DISCUSSION  OF  FINDINGS 


5.1   Geologic  Control  of  Groundwater 

Macrostructure  and  fractures  heavily  influence 
groundwater  occurrence  and  flow  in  the  study  area.  In 
particular,  the  Sulphur  Fault  combined  with  the  regional 
dip  to  the  northwest  on  the  upthrown  side  of  the  same 
fault,  greatly  influences  general  groundwater  flow  direc- 
tion to  freshwater  springs. 

The  Sulphur  Fault  invokes  high  permeability  in  the 
fine-grained  Arbuckle  dolomite  on  the  upthrown  side,  and 
creates  an  impermeable  boundary  in  the  Simpson  elastics 
and  calcarenites  (Figures  32,  33)  on  the  downthrown  side. 
This  statement  generally  agrees  with  the  literature 
(Davis  and  DeWiest  1966)  .  Also,  Ham  (1945)  observed 
small  pinnacles  of  orthoquartzite  of  the  usually  friable 
and  permeable  basal  sands  of  the  McLish  and  the  Oil  Creek 
formations,  adjacent  to  the  Sulphur  Fault. 

Fine-grained  dolomites  (Figure  38)  are  more  likely 
to  fracture  than  deform  (Stearns,  1967).  Dips  greater 
than  20°   in  the  Arbuckle  were  not  observed  by  this 
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investigator.  Field  reconnaissance  of  the  fractures  in 
the  Arbuckle  Group  substantiate  brittle  deformation. 
Conversely,  dips  of  the  Simpson  beds  in  the  grabens  of 
the  Sulphur  and  Hickory  Synclines  easily  surpass  20°. 
This  is  indicative  of  the  ductile  and  cataclastic  behav- 
ior of  the  Simpson  Group. 

The  potentiometric  map,  Figure  44,  shows  good  corre- 
lation with  formations  and  structure  (Figures  29  and  30) . 
Fairchild  et  al.  (1983),  concluded  that  the  structural 
geology  of  the  area  heavily  influences  groundwater  flow 
and  occurrence. 

The  relatively  high  regional  transmissivity  of  the 
Arbuckle  Group  can  be  attributed  to  solution  activities 
of  the  bedding  planes  and  fractures.  A  landowner,  Mr. 
Woods,  informed  this  investigator  that  drillers  on  his 
property  encountered  a  20-foot  thick  cavern,  in  the 
Arbuckle  Group,  at  SWi  NWj  SWj  SEj  of  Sec.  32,  TIN,  R4E. 
This  location  is  2  miles  east-northeast  of  the  freshwater 
springs  and  3/4  mile  north  of  the  Sulphur  Fault.  Fairly 
numerous  sinks  and  swallow  shafts  in  and  along  stream 
valleys  of  the  Arbuckle  Group  were  observed.  Rapid  rise 
in  water  levels  to  intense  precipitation  events  during 
the  fall  and  spring  can  be  attributed  to  fractures  and 
swallow  shafts  at  the  surface,  (high  vertical  perme- 
ability) ,  and  the  low  storage  coefficient  of  the  Arbuckle 
Group.   In  addition,  the   relatively  flat  water   table  of 
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the  Arbuckle  Group  (Figure  44)  could  conceivably  suggest 
karst  hydrology  (Fetter,  1980)  . 

Finally,  an  inferrence  of  the  Sulphur  Fault  through 
the  park  and  the  city  of  Sulphur  was  accomplished  (Fig- 
ures 29,  30)  .  Exact  location  is  not  known,  but  the 
inferrence  of  fault  is  quite  close  in  the  author's  opin- 
ion. The  fault's  trend  was  determined  by:  (1)  the 
geomorphology  of  the  park,  (2)  spring  occurrence,  (3)  the 
newly-mapped  formation,  (4)  lithology  of  the  observation 
wells,  (5)  piezometric  levels  and  hydraulic  gradients 
from  the  potentiometric  survey,  (6)  high  discharge  of  the 
Vendome  Well  relative  to  the  springs  in  the  park,  (7)  wa- 
ter quality  analysis,  (8)  downhole  TV  logs  of  the  Vendome 
Well  and  the  Horseman  Well,  and  (9)  surface  elevation 
tops  of  the  Arbuckle  Group  from  well  completion  report 
and  well  logs. 

5.2   Aquifer  Properties 

Simulation  of  drawdown  at  a  sink  node  demonstrated 
that  drawdown  at  the  wellbore  is  a  function  of  storage 
coefficient,   and  to  a   lesser  degree,   transmissivity . 

Transmissivity  values  obtained  from  the  simulation,  2.4 

2  2 

ft  /day  to  24  ft  /day,  are  fairly  low;  yet  drawdown  at 

the  West  Observation  Well,  approximately  .5  mile  from  the 
production  well,  experienced  negligible  drawdown.  How- 
ever, caution  must  be  exercised,  as  pump  tests  conducted 
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by  Gringarten  (1982)  in  a  karst  aquifer  demonstrated 
anisotropic  behavior.  Hence,  transmissivity  valves  for 
this  simulation  are  fairly  reliable  between  the  observa- 
tion well  and  the  production  well,  but  not  suitable  for 
transmissivity  in  other  directions. 

Also,  Fairchild  et  al.  (1983),  conducted  pump  tests 
in  the  Arbuckle  aquifer  and  observed  no  drawdown  in  ob- 
servation wells  at  distances  of  a  1/2  mile  from  the  sink 

-4 
point.   The  value  of  storage  coefficient,  1.7  x  10   ,  as 

determined  by  the  simulation,  is  reasonable,  and  compares 

_3 
favorably  with  a  storage  coefficient  of  1.2  x  10    at  the 

East  Observation  Well. 

The  major  limitations  to  this  simulation  are:  (1) 
that  the  Prickett  Lonnquist  model  (1975)  assumes  a  fully 
penetrating  well  in  a  confined  aquifer  at  some  specified 
sink  node;  (2)  insufficient  data  from  the  pump  test, 
(i.e.,  rate  of  drawdown);  and  (3)  homogeneous  aquifer  in 
close  proximity  to  the  well  bore. 

The  hydraulic  conductivity  value  of  1.6  ft.  per  day, 
obtained  from  the  slug  test  on  Crawford's  Well,  is  indic- 
ative of  the  high  transmissivity  at  this  locale.  The  re- 
gional hydraulic  conductivity,  as  determined  by  regional 
techniques  (Fairchild  et  al.,  1983)  of  this  area,  is  one 
order  of  magnitude  larger  than  that  determined  by  the 
slug  test.  This  is  generally  the  case,  as  hydraulic 
conductivity  determined  from  slug  tests  are  generally  one 
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order  of  magnitude  less  than  hydraulic  conductivity  de- 
termined by  pump  tests  that  stress  the  aquifer  on  a  much 

* 
larger  scale.    In  addition,  this  slug  test  was  performed 

in  an  area  that  was  determined  to  be  responsible  for  re- 
charge to  the  freshwater  springs  (Figure  30) . 

The  storage  coefficient  value  of  approximately  1  x 

_3 
10   ,  that  was  determined  for  the  East  Observation  Well, 

is  not  entirely  adequate.  The  lack  of  data  combined  with 
inoperative  water-level  recorders,  hindered  a  close  ap- 
proximation of  the  barometric  efficiency.  Under  the  cir- 
cumstances, a  fairly  reasonable  value  was  determine  be- 
cause of  the  following:  (1)  a  relatively  thin,  confining 
layer  of  120  ft.  thick  lies  unconformably  upon  the 
Arbuckle  Group;  and   (2)   the  East  Observation  Well  is 

located  in  close  proximity  to  the  recharge  area  that  has 

-3 
an  overall  storage  coefficient  of  8  x  10    (Fairchild  et 

al. ,  1983)  . 

The  average  aquifer  transmissivity  value  of  43,000 

2 
ft  /day,  as  determined  by  a  regional  technique  as  a  func- 
tion of  spring  discharge,  is  fairly  accurate  in  terms  of 
order  of  magnitude  (Fairchild,  et  al.,  1983).  The  re- 
gional technique,  a  modified  form  of  Darcy's  Law,  assumes 
that  Laplacian  conditions  prevail;  which  is  probably  not 
the  case  for  the  Arbuckle  aquifer. 


Investigator's  personal  experience  of  data  analysis  of 
slug  tests  and  pump  tests  on  the  same  well. 
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However,  Basmaci  et  al.  (1977),  assert  that  on  a  re- 
gional basis,  karst  aquifers  behave  like  granular  media. 
Then,  Basaci  et  al.  (1977) ,  maintain  that  Darcy's  Law  is 
applicable  with  the  assumption  that  some  error  (12-24%) 
is  committed  by  the  linear  flow  assumption  throughout  the 
fractured  karst  system. 

A  comparison  of  the  aquifer  transmissivity  responsi- 
ble for  the  spring  discharge  of  mixed  waters  in  the  cen- 
tral portion  of  the  park,  to  that  of  the  freshwater 
springs,  is  two  orders  of  magnitude  less.  This  is  con- 
firmed by  comparatively  high  hydraulic  gradients  and  much 
smaller  discharge  amounts.  These  findings  agree  quite 
favorably  with  the  commentary  of  Davis  and  DeWiest  (1966) 
on  spring  discharge  amount  as  a  function  of  aquifer  per- 
meability. 

5.3   Data  Analysis  of  Wells  and  Springs 

Well  hydrographs  firmly  establish  a  relationship  be- 
tween spring  discharge  and  piezometric  elevation  (Figure 
42) .  Variations  of  piezometric  levels  in  the  East  Obser- 
vation Well  are  very  similar  to  water-level  fluctuations 
in  the  outcrop  area  of  the  Arbuckle  Group,  east  of  Sul- 
phur (Table  26).  Previous  investigators,  Hart  (1972)  and 
Harp  et  al.  (1976),  also  concur  with  these  observations. 

The  West  Observation  Well  hydrograph  is  considerably 
dampened  with  respect  to  the  hydrograph  of  the  East 
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Observation  Well  (Figure  42)  .  This  can  be  attributed  to 
decreased  permeability  or  increased  storativity  of  the 
Arbuckle  Group.  The  hydraulic  gradient  increases  to  40 
ft.  per  mile,  down  dip  from  the  East  Observation  Well; 
which  indicates  decreasing  permeability  in  the  direction 
of  the  West  Observation  Well. 

Also,  the  slug  test  that  was  performed  on  the  West 
Observation  Well,  revealed  an  extremely  slow  rate  of  rise 
after  the  probe  was  withdrawn.  This  provides  additional 
substantiation  that  the  permeability  of  the  Arbuckle 
Aquifer  is  considerably  reduced  in  the  vicinity  of  the 
West  Observation  Well. 

The  long-term  trend  of  the  hydrograph  can  be  used  to 
ascertain  deprivation  effects  of  adjacent  production 
wells  (Figure  42) .  The  deprivation  slope  or  the  slope  of 
the  long-term  trend  of  the  hydrograph  was  determined  to 
be  .65  ft/year  for  both  observation  wells  (Figure  42). 
It  was  also  established  by  the  use  of  the  inch-year  time 
bar  chart,  that  the  primary  cause  of  increasing  cessation 
of  spring  flow  is  attributed  to  nearby  production  wells. 
Finally,  assuming  that  the  deprivation  slope  remains  con- 
stant, cessation  of  all  spring  flow  will  occur  in  40 
years  (Equation  4.1). 

The  variability  of  Antelope  and  Buffalo  Springs  is 
high;  indicating  that  both  derive  their  waters  from  the 
Arbuckle  Group  (Table  30)  .   The  source  of  the  freshwater 
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springs  can  be  further  correlated  by  well  hydrographs 
that  vary  as  much  as  50  feet  per  year  in  the  outcrop  area 
of  the  Arbuckle  Group,  east  of  Sulphur  (Fairchild  et  al., 
1983)  .  Conversely,  the  variance  of  water  levels  in  the 
Simpson  Group  is  less  than  10  feet  per  year  (Fairchild  et 
al. ,  1983) . 

The  exhaustion  coefficient  e,  indicates  whether  con- 
duit or  diffuse  flow  predominates  in  a  carbonate  aquifer 
(White  et  al.,  1974).   Another  variable,  response  time 

t  ,  can  also  be  utilized  to  assess  conduit  or  diffuse 
r 

flow.  White  (1977)  asserts  that  conduit  flow  prevails 
when  the  recession  limb  of  a  spring  hydrograph  is  similar 
to  a  recession  limb  of  a  stream  hydrograph.  Visual 
inspection  of  the  recession  limbs  of  the  well  hydrographs 
and  spring  hydrographs,  Figure  42,  indicates  a 
predominately  diffuse-flow  system.  Values  of  response 
time  t  ,  listed  in  Table  30,  most  nearly  coincides  with 
the  values  for  intermediate-response  springs  listed  in 
Table  4  (White,  1974)  . 

Statistical  analysis  revealed  that  the  peak  dis- 
charge of  Buffalo  and  Antelope  Springs  closely  coincides 
with  the  peak  precipitation  events  in  spring  and  fall 
(Table  29).  White  et  al.  (1974)  maintains  that  peak  dis- 
charges of  diffuse-flow  springs  will  coincide  with  the 
mean  spacing  of  maximum  precipitation  events  on  a  season- 
al basis.   White's  preceding  statement  agrees  favorably 
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with  these  findings,  and  lends  further  credibility  that 
the  Arbuckle  Group  is  primarily  a  diffuse-flow  system. 

5.4   Analysis  of  Spring  Flow  Occurrence  and  the  Recharge 
Areas 

The  freshwater,  artesian  springs,  Antelope  and 
Buffalo,  discharge  from  orifices  of  the  Pontotoc  Group. 
The  source  of  water  for  both  springs  is  from  fractures 
and  joints  in  the  Pontotoc  Group  that  extend  downward  to 
the  top  of  the  pressurized  Arbuckle  Aquifer.  The  thick- 
ness of  the  confining  layer,  the  Pontotoc  Group  is  be- 
lieved to  be  less  than  50  feet  in  the  vicinity  of  Buffalo 
and  Antelope  Springs.  This  thickness  was  deduced  from  a 
lithologic  log  of  the  East  Observation  Well  (Moore, 
1972) . 

The  Sulphur  Fault  is  inferred  with  virtual  certainty 
just  south  of  Buffalo  Spring  (Figure  30)  .  This  can  be 
readily  justified  by  the  following:  (1)  lithology  of  the 
East  Observation  Well;  (2)  water  quality  analysis  con- 
ducted by  Hart  (1972)  and  Harp  et  al.  (1976);  (3)  the 
newly  mapped  outcrop  of  the  Simpson;  (4)  spring  discharge 
as  a  function  of  aquifer  transmissivity ;  and  (5)  orien- 
tation of  fracture  sets  in  the  Pontotoc  Group  that  trend 
in  a  westerly  direction,  located  just  east  of  the  park 
boundary  (Table  9)  .  This  is  also  in  agreement  with  Hart 
(1972) ,  who  previously  inferred  the  Sulphur  Fault  to 
Buffalo  Spring. 
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The  freshwater  springs  evolved  as  a  result  of  in- 
tense fracturing  and  subsequent  solution  of  the  brittle 
Arbuckle  dolomites  in  close  proximity  to  the  Sulphur 
Fault.  In  addition,  the  permeable  Arbuckle  Group  on  the 
upthrown  side  of  the  Sulphur  Fault  is  offset  next  to  the 
less  permeable  Simpson  Group  on  the  downthrown  side. 

The  recharge  area  for  Buffalo  and  Antelope  springs 
is  in  the  outcrop  of  the  Arbuckle  Group,  east  of  Sulphur 
(Figure  30)  .  In  addition,  some  recharge  probably  occurs 
from  the  Pontotoc  that  overlays  the  Arbuckle  Group, 
located  between  the  springs  and  the  surface  contact  of 
the  two  geologic  units  to  the  east. 

The  recharge  area  was  determined  by  general  flow  di- 
rections from  the  potentiometric  map  and  structural  geol- 
ogy of  the  project  area.  It  is  believed  that  the  ground- 
water flow  is  highly  anisotropic  along  the  Sulphur  Fault. 
The  dip  of  the  Arbuckle  Group  is  to  the  west  and  north- 
west in  the  study  area,  which  also  impacts  on  groundwater 

flow  direction  to  the  artesian  springs. 

2 
The   recharge   area   is   approximately   5   mi    for 

Antelope   and  Buffalo   Springs.    This   estimate   agrees 

favorably  with  the  nomograph  devised  by  Davis  and  Dewiest 

(1966),  Figure  46;  for  values  of  3  cfs  and  4.6  inches  of 

annual  recharge.   Also,  the  primary  recharge  area  for  the 

remainder  of  the  springs  in  the  park  is  estimated  to  be 

2 
less  than  2  mi  ;  in  the  outcrop  of  the  Simpson  Group 
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along  the  periphery  of  the  Sulphur  Syncline,  south  and 
southeast  of  the  park.  Finally,  the  estimate  of  the  to- 
tal recharge  area  responsible  for  spring  flow  in  the  park 
is  considerably  less  than  originally  believed  by  previous 
investigators,  Hart  (1972)  and  Harp  et  al.  (1976). 


CHAPTER  VI 


SUMMARY  AND  CONCLUSIONS 


The  origin  of  the  freshwater  artesian  springs, 
Buffalo  and  Antelope,  is  attributed  to  the  Sulphur  Fault 
crossing  hard,  brittle  dolomite  of  the  Arbuckle  Group  and 
thereby  creating  zones  of  greater  permeability.  In  addi- 
tion, it  is  believed  that  the  fault  has  offset  the 
fractured  and  permeable  Arbuckle  Group  against  the  less 
permeable  rocks  of  the  Simpson  Group. 

The  freshwater  springs  eminate  from  orifices  of  the 
Pontotoc  Group.  The  springs  derive  their  water  from  the 
confined  Arbuckle  Aquifer  that  flows  upward  by  pressure 
through  fissures  and  joints  in  the  Pontotoc  Group,  ulti- 
mately discharging  at  the  surface. 

Structure  and  lithology  govern  groundwater  occur- 
rence in  the  study  area.  The  Sulphur  Fault  that  trends 
west-northwest  separates  groundwater  flow  to  the  park. 
The  freshwater  artesian  springs,  Antelope  and  Buffalo, 
receive  recharge  in  the  form  of  precipitation  that  falls 
on  the  outcrop  of  the  Arbuckle  Group  east  of  Sulphur, 
then  infiltrates  to  the  water  table,  and  moves  downdip  to 
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the  springs  (Figure  30)  .  The  remainder  of  the  mineral- 
ized spring  in  the  park,  receive  recharge  primarily  from 
the  outcrop  of  the  Simpson  along  the  periphery  of  the 

Sulphur  Syncline  (Figure  30)  .    The  area  of  recharge  is 

2 
approximately  5  mi.   for  Antelope  and  Buffalo  Springs. 

Also,  the  recharge  area  for  the  mineralized  springs  is 

2 
estimated  to  be  2  mi .  . 

Geological  reconnaissance  of  the  Arbuckle  Group, 
east  of  Sulphur,  revealed  intensely  fractured  microcrys- 
talline  dolomite,  with  a  fair  amount  of  sinks  and  swallow 
shafts  in  the  stream  valleys.  The  Arbuckle  Group  is  rel- 
atively flat,  gently  dipping  to  the  northwest.  In  addi- 
tion, the  inferrence  of  the  Sulphur  Fault  into  the  park 
boundary  was  greatly  enhanced  by  the  discovery  of  a  small 
outcrop  of  the  Simpson  Group,  east-southeast  of  Buffalo 
Spring. 

Regional  transmissivity  of  the  Arbuckle  Aquifer  re- 
sponsible for  groundwater  flow  to  Antelope  and  Buffalo 

Springs   was   determined   to   be   approximately   43,000 

2 
ft  /day.   Hydraulic  conductivity  of  1.6  ft/day  was  de- 
termined from  results  of  a  slug  test  on  a  selected  well 

in  the  outcrop  area  of  the  Arbuckle  Group.   The  values  of 

-3 
storage  coefficient  ranged  from  about  1.2  x  10    to  1 . 7  x 

-4 
10    for  the  confined  Arbuckle  Aquifer  in  close  vicinity 

to  the  springs . 

Well   and   spring   hydrographs   provided   additional 
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insight  to  the  subsurface  hydrology  of  the  Sulphur  area. 
The  peak  discharge  of  the  springs  coincided  with  the  max- 
imum precipitation  event,  and  the  maximum  piezometric 
elevation  in  the  observation  wells.  The  slope  of  the 
long-term  trend  of  the  hydrograph  was  ascertained  to  be 
.65  ft/year.  This  slope  indicates  that  the  piezometric 
elevation  of  the  Sulphur  area  is  declining  at  a  yearly 
rate  of  .65  ft.  If  this  slope  value  remains  constant, 
permanent  cessation  of  flow  of  the  freshwater  springs 
will  occur  in  40  years. 

Variance  values  of  the  freshwater  artesian  springs 
are  in  excess  of  390%  for  at  least  a  period  of  12  years. 
Thus,  confirming  that  both  springs  derive  their  waters 
from  the  Arbuckle  Group,  an  aquifer  of  low  storativity. 
Also  statistical  analysis  revealed  that  the  peak  dis- 
charge occurs  in  May  and  the  minimum  discharge  in  August. 
The  exhaustion  coefficient  values  ranged  from  .49  to  .6 
mos.  for  Antelope  Spring.  These  values  are  typical  for 
intermediate  response,  diffuse-flow  springs  that  origi- 
nate from  carbonate  aquifers. 

The  annual  hydrologic  budget  for  the  Sulphur  area  is 

2 
based  on  15  mi.   of  recharge  area  with  an  average  annual 

recharge  rate  of  4.6  inches  per  year.   The  total  amount 

of  water  available  for  use  on  a  average  yearly  basis  is 

4  9 

approximately  1.2  x  10  gallons.  Current  analysis  of 
production  wells  and  flowing  wells  in  the  Sulphur  area, 
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q 
reveals  that  withdrawal  of  water  exceeds  1.4  x  10   gal- 


lons per  year. 


NOMENCLATURE 


L  =  length 

T  =  time 

F  =  force 

M  =  mass 

h    =   head  [L] 

2 
A   =   area  [L  ] 

-2, 
P    =   pressure  [FL   ] 

q   =   specific  discharge  [LT   ] 

_3 

Y  =   specific  weight  [FL   ] 

-3 
p    =   density  [ML   ] 

v   =   velocity  [LT   ] 

K    =  hydraulic  conductivity  [LT   ] 

y    =   dynamic  viscosity  [ML   T   ] 

_2 
g    =   gravitational  acceleration  constant  [LT   ] 

h   =   pressure  head  [L] 
P 

z    =   elevation  head  [L] 

3 

V  =   volume  of  the  pores  [L  ] 


P 
d   =   diameter  [L] 


3 
V_   =   bulk  volume  [L  ] 


3  -3 
n    =   porosity  [L  L   ] 

3-1-3  1 
K    =   recession  coefficient  [L  T   L   T  ] 
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NOMENCLATURE     (continued) 


3 

V  =   volume  of  solids  [L  ] 
s 

I*  =  Inflow  [L3] 

O'  =  Outflow  [L3] 

S'  =  Storage  [L3] 

S  =  storage  coefficient  [LL   ] 

5  =  specific  storage  [L   ] 

S 

2  -2 

$    =   force  potential  [L  T   ] 

r  2  -1, 
T    =   transmissivity  [L  T   J 

-2 
P    =   pressure  [FL   ] 

-1  2 

6  =  water  compressibility  [F   L  ] 

-1  2 
a    =   aquifer  compressibility  [F   L  ] 

b   =   aperture  width  of  fracture  [L] 

Kf   =  hydraulic  conductivity  of  a  single  fracture  [LT   ] 

2 
K. .  =   intrinsic  permeability  [L  ] 

q.   =   fracture-flux  vector  per  unit  width  of  conduit 

[lV1] 

3 

V  =   volume  of  sample  chamber  [L  ] 

3 
V,   =  volume  of  evacuated  chamber  [L  ] 
b 

-2  -1  2 
BE   =   barometric  efficiency  [FL   F   L  ] 

r  =  radius  [L] 

x  =  coordinate  in  the  horizontal  axis  [L] 

y  =  coordinate  in  the  vertical  axis  [L] 

R  =  radius  of  influence  [L] 
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NOMENCLATURE  ( con t  inued ) 


_2 
a*   =   constant  aquifer  characteristic  [TL   ] 

Q'  =  total  average  annual  runoff  [L] 

ET  =  average  annual  evapotranspiration  [L] 

P'  =  average  annual  precipitation  [L] 

AP  =  change  in  barometric  pressure  [L] 

b1  =  thickness  of  a  confined  aquifer  [L] 

h'   =   saturated  thickness  of  an  aquifer  [L] 

-3 
E    =   total  energy  per  incremental  volume  [FLL   ] 

t    =   time  [T] 

Q   =   (-)  recharge  or  (+)  withdrawal  [LT   ] 
n 

3  -3 
S    =   specific  yield  [L  L   ] 


2  -1  -1 
K   =   normal  stiffness  factor  [L  T   F   ] 


y 

"n 
S-   =   fracture  spacing  [L] 

I    =   hydraulic  gradient  [LL   ] 

t    =   response  time  [T] 

q2   =   space-averaged  velocity  of  the  fracture  relative 

to  the  matrix  [LT   ] 

r    =  mass  flow  rate  from  block  to  fracture  [MT   ] 

e    =   exhaustion  coefficient  [T   ] 

k,   =   diffusivity  constant  [T   ] 

3-1-3  1 
V'   =   variability  or  spring  variance  [L  T   L   T  ] 

2 
K.   =   intrinsic  permeability  of  the  matrix  [L  ] 

2 

K2   =   intrinsic  permeability  of  the  fracture  [L  ] 

s    =   drawdown  [L] 
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NOMENCLATURE  (continued) 

b   =  width  of  weir  [L] 
w 

H    =   depth  of  pool  at  a  specified  point  [L] 
k'   =  weir  constant,  unitless 
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APPENDIX  I 

CONTINUITY  PRINCIPLE 
(Fetter,  1980) 


Fluid  moves  in  only  one  direction  through  the  ele- 
mental volume.  The  fluid  movement  can  be  subdivided  in 
three  component  flows  aligned  with  the  three  principal 
axes  (refer  to  Figure  13) .   Flux  along  the  x  axis  in  the 

elemental  cube  is  p  q   dydz  (units  in  FTL   ) .   Mass  flux 

w^x   J  ' 

(along   the  x   axis)  out   is   p  q   dydz  +  -r-(p  q  )dxdydz. 

^  Kw^x  J  3x   w^x    J 

The   net,   accumulated  mass   flux   is   influx-out flux, 

-  -r—  (p  q  )  dxdydz .   It   follows  that   the  other   two  net 
3x   Kw^x    J 

fluxes,  y  and  z  components   of  flux,  are  identical  in  de- 
rivation:   (1)  y  axis,  -  -r—  (p  q  )dxdydz;  and  (2)  z  axis, 

-  T-r   (p  q  ) dxdydz.   Combining  all  three  terms  yields  the 

following  equation: 

3  3  9 

-(  —  p  q   +  —  p  q   +  —  p  q  )  dxdydz  (1-1) 

3x   w  x    3y  w  y    3z   w  Z 

The  volume  of  water  in  storage  (fully  saturated)  is 

equal  to  ndxdydz,  and  the  mass  of  water  in  storage  is  p  n 

dxdydz.   The  mass  of  the  matrix  of  the  elemental  volume 

is  (1-n) dxdydz.   Change  of  the  mass  of  water  with  respect 

to  time  is: 
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gM    3   (p  n  dxdydz) 
—  =  —   W  (1-2) 

st   at 

The  fluid  density  and  porosity  will  change  when 
pressure  changes  occur  in  the  elemental  volume.  The  com- 
pressibility of  water,  3,  is  the  rate  of  change  in  den- 
sity with  a  corresponding  change  in  pressure,  P: 

BdP  =  dp  /p  (1-3) 

w  w 

In  addition,  changes  in  volume  occur  when  pressure 

varies.   Due  to  lateral  constraining  forces,  changes  in 

volume  will  primarily  occur  in  the  z  direction.   The 

aquifer  compressibility,  a,  is: 

d(dz) 

adP  =  (1-4) 

dz 

Expansion  and  contraction  of  the  aquifer  changes  po- 
rosity, but  the  volume  of  the  matrix  remains  constant. 
If  deformation  is  in  the  z  axis  only,  then  differentia- 
tion of  the  matrix  volume  with  respect  to  time  yields: 

dzdn  =  (l-n)d(dz)  (1-5) 

and 

(l-n)d(dz) 

dn  =  (1-6) 

dz 

The  pressure,  P,  at  some  point  in  the  aquifer  is 

equal  to  P   +  p  gh,  where  P   is  equal  to  atmospheric 

pressure  (a  constant) ,  and  h  is  the  height  of  the  column 

of  water  above  that  point.   Thus,  dP  =  p  gdh,  and  sub- 

w 

stituting  equations  (1-3)  and  (1-4),  respectively,  into 
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the  preceding  equation  yields  the  following  equations: 


and 


dpw  =  p  B(p  gdh)  (1-7) 

W      WW 


d(dz)  =  -dza(p  gdh)  (1-8) 


Equation   (1-5)   can  be  rearranged  by  substituting  the 

equivalent  of  d(dz)  of  Equation  (1-8): 

dn  =  (n-l)ap  gdh  (1-9) 

The  change  of  mass  with  respect  to  time  (dx  and  dy 

are  constant)  yields  another  version  of  Equation  (1-2) 

(carrying  out  the  differentiation) . 

3M         3 (dz)         9n        3p 

—  =  {p  n  +  p  dz  —  +  ndz  — -}dxdy       (1-10) 

3t      W   3t        W    3t        3t 

Substitution   of   Equations   (1-7)   through   (1-9)   into 

Equation  (1-10)  with  some  algebraic  manipulation  yields: 

3M  3h 

—  =  (ap  g  +  ngp  g)  p  dxdydz  —  (1-11) 
3t               ww        3t 

Set  Equation  (1-2)  equal  to  Equation  (1-11) ;  where  Equa- 
tion (1-2)  is  net  accumulation  of  material,  and  Equation 
(1-11)  is  the  change  in  mass  with  time: 

3(q  )    3(q  )    3  (q  ) 

-  {  —  +  *_  +  — }   p  dxdydz 

3x       3y       3z        W 

(ap    g   +   nBptTg)piT   dxdydz    3h/..  (1-12) 

W  WW  dt 

From  Darcy's   Law: 

3h 
qx    =    -K   —  (1-13) 

X  3x 
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3h 

q      =   -K  — 

Y            ay 

and 

3h 

qz    =   "K  — 
3z 

(1-14) 


(1-15) 


Substituting   Equations    (1-13)    through   (1-15)    into 

Equation  (1-12)  yields: 

2      2      2 
3  h    3  h    3  h  3h 

K  ( — 5-  +  — ^  +  — ~)  =  (ap  g  +  n6p  g)  —        (1-16) 

3xZ    3y     3z         W       W    3t 

which  is  the  governing,  three-dimensional  flow  equation 

for  an  isotropic,   homogeneous  confined  aquifer.    For 

two-dimensional  flow  with  no  vertical  flow  component  and 

the  storage  coefficient,  S   =  m(ap  g  +  nBp  g) ,  and  Trans- 

missivity  T  =  Km  where  m  is  a  constant  saturated  aquifer 

thickness;  Equation  (1-16)  can  be  rearranged  as  follows: 

32h    32h    S    3h 

— 2    +   2  =  —  ~"  (1-17) 

3x     3y     T    3t 

Another  equation  that  is  used  when  steady-state  con- 
ditions prevail  in  a  homogeneous  porous  media  is  the 


LaPlace  equation: 

=  0  (1-18) 


2      2      2 
3  h    3  h    3  h 


or 


2      2      2 
3x     3y     3z 


V2h  =  0  (1-19) 
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APPENDIX  II 
PRICKETT-LONNQUIST  MODEL 

Prickett  and  Lonnquist  (1971)  devised  a  finite  dif- 
ference model  capable  of  handling  transient,  two- 
dimensional  flow  equation  for  homogeneous  or  nonhomoge- 
neous,  isotropic,  non-leaky  confined  aquifer  (Figure 
A.l).  Also,  Prickett  and  Lonnquist  (1971)  outlined  other 
modifications  to  the  model  that  are  capable  of  solving 
phreatic  aquifer  behavior,  leaky  confined  aquifer, 
confined  aquifer  goes  to  phreatic  aquifer,  and  evapo- 
ration of  groundwater. 

The  model  requires  certain  input  parameters,  such  as 
initial  conditions  of  head  at  each  node,  and  boundary 
conditions  (if  applicable)  ,  such  as  no  flow  or  constant 
head.  Additional  input  parameters  include  recharge,  dis- 
charge, storage  factor  (storage  coefficient  x  incremental 
area  (AxAy)) ,  and  transmissivity  in  the  x  and  y  direction 
for  each  node.  Parameter  data  and  data  array  are  also 
required. 

Data  array  is  simply  the  number  of  columns  and  the 
number  of  rows.  The  model  is  designed  for  a  50  x  50  ar- 
ray or  less.   The  parameter  input  includes:   the  number 
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60 
70 
C 


ILLINOIS     STATE    WATER     SURVEY 
BASIC     AOIIIFER     SIMULATION     PROGRAM 


DEF INiriON    OF     VARIAOLES 

Mflll.JI HE40S     4T     START     OF     TIME 

INCREMENT     II,  Jl 
HII.JI HEADS     AT     FNO    OF     TIME 

INCRIMINT     IFI  I 
S >-  1  I  I  . J  1 SIDRAGE     FACTOR     FUR 

ARTESIAN    CON!)  IT  IONS 

IG4L/FT I 
01  I  ,  Jl  — CIINSIANT     WITHIIK4.JAL 

RATES     IGPOI 
Tl  I  ■  J.  1  I AHIIIIFR     TRANSMISSI  VI  TV 

HF  TWFI  N     I  p.l     AVI     I  .  J*  I 

11.41  /II4Y/FT  I 
II  I  .J  .21 Auul^FR    IRANS-IISSIVITY 

HEIwEFN    I ,J    MID    I ♦ I.  J 

I GAL/04Y/F  T  I 
44,HB,CC.D0-CUFEF|C IENTS   IN  WATER 

B4LANCE  E0U4TIONS 

Ml; NO.  OF  ROWS  !N  MODE'. 

NC NO.    OF     COLUMNS     IN    MODEL 

NSTEPS NO.     OF     TIMF     INCREMENTS 

DELTA T|MF     INCREMENTS     ID4YSI 

MM,Sl,0iltTT-DFF4OLT     VALMF  S 

I MODEL     COLUMN    NUMRFR 

J MODEL    ROW    NUMBER 


DIMENSION    HI50.51  l,HOI50,SOI  . 
1SF 1 1  SO ,50  I, 0(50. SO  I .11  SO, SO.? I. 
2BI  SOI  ,r,ISOI,DL  ISO. 50  I 

10RN    OFF     UNDERFLOW     TRAP 

CALL    ERRSFT I20R. 256.-1 .1  I 

DEFINF     INPUT     ANO     OUTPUT     DE  V I C  F     NUMBERS 

INTEGER     OUT 
.   IN. 5 
0UT.6 

READ     PARAMETER    CARD    4N0 
DEF4ULT     VALUE    CARO 

RE40IIN.IOINSTEPS.DEL  T A, ERROR. 
INC .NR ,TT ,S1 .HH.OO 
FORMATI 16. 2F6. 0/216, *F6.0) 

Fill.     ARRAYS     WITH    DEFAULT     VALUES 

00    ?0     l-l.NC 

00  20    J=l .NR 

1 1 l .J, I  |«TI 
Tl l,J.2l=TT 
SF1I I.JI-S1 

HI    I    .J  I.  II H 
HOI  I  ,JI=MH 

01  I    ,11-JU 

Ml  All     NODE     C4RIIS 

RFADI  IN,*0,FND=50  I  I  ,  J.T  I  |,  J,  I  I  , 
lit I  ■  J.2  1  ,SF1 l I . JI.hi I , j I, u i I, j l 
FORMATI  2I3«2F6.0.2F*.0,1F6.0I 
GO    III    30 

START     OF     SIMULATION 

T  IMF-0 

DO    120    ISTf  P'I.N',11  PS 

T  IME-T IME'OELTA 

PREOICT    MEAOS    FOR    NE«T 
TIME     INCREMENT 

DO    70    Ll.NC 

DO    70    J-I.NR 

0"MI I .ji-hoi  |  . j I 

HOI  I , j i.hi | . jl 

F.1.0 

IFIOLI I .JI.EO.O.OIGO   TO    60 

IFI  ISTEP.GT.2  if.o/ch.  I  I  .Jl 

IFIF.GT.5IF.5.0 

IFI  F.LT.O.OIF.0.0 

DLI  I.J1.0 

hi  I  , JI.hi I  ,J|.0»« 


10 

10 

UMBER 

NUMBER 

C 

REFINE    ESTIMATES    OF    HEADS   8Y    IA01    METMOO 

009* 

0001 

c 

0095 

0002 

ITER-0 

0096 

000  J 

80 

E-0.0 

0097 

000k 

ITFR-ITER.l 

0098 

0005 

c 

0099 

0006 

c 

COLUMN    CALCULATIONS 

0100 

000  7 

c 

0101 

OOOR 

00     \90     llaliNC 

0102 

0(11)1 

I'll 

0103 

0010 

IF (MODI I:TFP.1T6R,2).F0.1I     1 aNC-l*! 

010* 

0011 

OH     1 70     J« 1 .NR 

0106 

0012 

c 

0106 

001  1 

c 

CALCULATE     8     ANO    G    ARRAYS 

0107 

noi* 

c 

0108 

0(115 

BB'SFI 1 1 , JI/DELTA 

0109 

on  1 6 

(lO-lilll  1  ,JI«SC1  (  1  .  J  l/DFLTA-qi  |  ,J| 

0110 

Out  I 

AA.0.0 

0111 

001  A 

CC-0.0 

0112 

0019 

IFI J-| 190. 100,90 

0113 

0020 

90 

AA.-TI I.J-1,1 I 

on* 

0021 

BB=BB.T 1 1 ,J-l,l I 

0115 

0022 

100 

IFI J-NRI 1 10,120, 1 10 

0116 

002  J 

110 

CC--TI  I.J.I  1 

0117 

002* 

BB-HS.  I  1   1  ,J,|  1 

0118 

0025 

120 

IFI l-ll 130,1*0,130 

0119 

0026 

no 

BB.HB.TI I-I.J.2) 

0120 

002  7 

00.00. HI I-l,JI»TII-l,J,2  1 

0121 

0028 

1*0 

IFI l-NC) 1 50.160,150 

0122 

0029 

160 

BB-BB.TI 1 ,J,2I 

0123 

0010 

oil. Oil. Ml  l»l,JI»Tt  1  ,J,2) 

012* 

00)1 

160 

U.On-AA.RI J- | | 

0125 

0032 

Bl JI-CC/W 

0126 

0013 

170 

Gl  Jl  =(  UO-AA.i.l  J-  1  1  1  /M 

0127 

003* 

c 

0128 

00)5 

c 

RE-ESTIMATE     MEAOS 

0129 

0016 

c 

0130 

0017 

F'F.ABSIHI | ,NRI-GINRI 1 

0131 

0018 

HI  I .NRI'CINRI 

0132 

0019 

N.NR-1 

0133 

00*0 

IRO 

HA.GI  NI-HINI«M( | ,N<1  ) 

013* 

00*1 

E  .E.A8SIHA-HI  |   .HI  1 

0135 

00*2 

HI  1  ,N).HA 

0136 

00*3 

N-N-l 

0137 

00** 

IFI Ml  190,190,180 

0138 

00*5 

190 

CONTINUE 

0139 

00*6 

c 

01*0 

00*7 

c 

ROW    CALCULATIONS 

01*1 

00*8 

c 

01*2 

00*9 

00    300    JJ*1,NR 

01*3 

0050 

J.JJ 

01** 

0051 

IFI MODI ISTEP»ITER,2).=0.1  1    J-NR-J*! 

01*5 

0052 

00    280    l=l.NC 

01*6 

0053 

IIR.SF1  1  1, JI/DELTA 

01*7 

005* 

DlLHOI  1  ,  JI«SF  1  1  |  ,  Jl/OELTA-QI  1  ,J1 

01*8 

0065 

AA'0.0 

01*9 

0066 

CC'0.0 

0150 

0067 

IFI J- 1)200,2 10, 200 

0151 

0068 

200 

HhsBB.I 1 1 ,J-1 ,11 

0152 

0059 

DO'OD.HI  1  ,J-ll«T  II . J-l, |  1 

0153 

0060 

210 

IF| J-NRI .^O, 210, 220 

015* 

0061 

2  20 

OO'DO.HI 1 ,J«1 )»T(I  ,  J,  1  » 

0155 

0062 

BB=8H»TI 1 ,J,l 1 

0156 

0063 

2'0 

IFI  1-1-12*0,250,2*0 

0137 

006* 

2*0 

BR.BH.TI I-I.J.2I 

0158 

0(1*.  6 

AA.- II 1-1, J, 21 

0159 

0066 

250 

IFI l-NCI?60,270,260 

0160 

0067 

260 

BH.BB. T 1 | ,J,2) 

0161 

0068 

CC'-TI l,J,2) 

0162 

0069 

270 

W.nB-AA»B( l-ll 

0163 

00  70 

Bl  1  1-CC/W 

016* 

0071 

280 

Gl ll-l D0-AA»GII-1)I/W 

0165 

0072 

C 

0166 

0073 

c 

RE-ESTIMATE    HEADS 

0167 

007* 

c 

0168 

0076 

E.E.AHSIHINC.JI-GINCII 

0169 

00  76 

HI NC , JI.GINC 1 

0170 

0077 

N.NC-1 

0171 

0078 

290 

HA.GI Nl -BINI«H|N. 1 .J ) 

0172 

0079 

E.E> ABSIHIN, JI-H4I 

01T3 

0080 

HIN,J)«HA 

017* 

0081 

N-N-l 

0175 

0082 

IFI Nl 300,300.290 

0176 

008) 

300 

CONTINUE 

0I7T 

0084 

IFIE.GT. ERROR)    GO    TO    80 

0178 

0085 

c 

0179 

008* 

c 

PRINT    RESULTS 

0180 

0087 

c 

0181 

0086 

WRITEIOUT,310)T|NE,E.ITER 

0182 

0084 

310 

F0RMATC6H2T  IMf  ..F6.2///.E20.7,  15  1 

0183 

0090 

0ELTA»0£LTA«1.2 

018* 

0091 

00    320    J-l.NR 

018* 

0092 

320 

WRIT(IOUT,330IJ,IHIt  ,jl  ,l*| ,NC  1 

018* 

009) 

330 

FORMAT  1  15,51, IOF 10.*/!  1211  OF  10.*)) 

01ST 

C 

0186 

STOP* 

016V 

ENO 

0190 

Figure  A.l.      Prickett-Lonnquist  Model   for   an   isotropic, 
nonleaky,     confined     aquifer    (After 
Prickett    and   Lonnquist,    1971) . 
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of  increments  NSTEP,  the  time  increments  DELTA,  and  the 
ERROR  value  for  closure.  The  model  is  designed  for  in- 
cremental time  increases  as  depicted  in  Table  A.l.  For 
example,  an  NSTEP  value  of  10  which  is  equivalent  to 
ISTEP  value  of  10  in  Table  A.l;  time  =  ISTEP  =  25.96 
hours  or  25.96  days  for  a  DELTA  value  of  1  hour  or  1  day. 
If  one  chooses  to  model  to  the  exact  time,  line  184  can 
be  deleted.  NSTEP  and  DELTA  would  be  specified,  result- 
ing in  data  output  for  head  at  a  constant  increment  of 
time  to  the  desired  time. 

The  flow  chart  for  the  confined  aquifer  simulation 
model  is  shown  in  Figure  A. 2.  The  program  basically  con- 
sists of  (1)  variable  explanation  in  lines  1  through  33; 
(2)  parameter  data  stored,  and  data  array  filled  with  de- 
fault values  and  node  values  in  lines  34  through  72;  (3) 
computation  of  data  in  lines  73  through  178;  and  (4) 
print  out  of  results  in  lines  180  through  187. 

Finally  the  preceding  discussion  is  just  a  brief 
overview  of  the  Prickett-Lonnquist  model.  For  an 
in-depth  explanation  of  the  simulation  model,  the  reader 
is  referred  to  the  original  article  put  forth  by  Prickett 
and  Lonnquist  (1971) . 
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Table  A.l.   ISTEP  Function  Values 


ISTEP 

ISTEP 

ISTEP 

Function 

ISTEP 

Function 

1 

1.00 

26 

567.37 

2 

2.20 

27 

681.84 

3 

3.64 

28 

819.21 

4 

5.37 

29 

984.05 

5 

7.44 

30 

1181.87 

6 

9.93 

31 

1419.24 

7 

12.92 

32 

1704.08 

8 

16.50 

33 

2045.90 

9 

20.80 

34 

2456.08 

10 

25.96 

35 

2948.29 

11 

32.15 

36 

3538.95 

12 

39.58 

37 

4247.74 

13 

48.50 

38 

5098.28 

14 

59.20 

39 

6118.93 

15 

72.03 

40 

7343.71 

16 

87.44 

41 

8813.45 

17 

105.93 

42 

10577.14 

18 

128.12 

43 

12693.57 

19 

154  74 

44 

15233.28 

20 

18b.69 

45 

18280.93 

21 

225.02 

46 

21938.11 

22 

271.03 

47 

26326.73 

23 

326.23 

48 

31593.07 

24 

392.48 

49 

37912.68 

25 

471.98 

50 

45496.21 

(after  Prickett  and  Lonnquist,  1971) 
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Figure  A. 2 


Simplified  flowchart  of  the  Prickett- 
Lonnquist  Model  (After  Prickett  and 
Lonnquist,  1971) . 
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APPENDIX  III 
LITHOLOGY  AND  STRUCTURE  OF  THE 
SIMPSON  GROUP  OUTCROP 

Station  1.  Location  -  NWj  NWj  NWi  SWi  of  Sec.  6, 
T1S,  R4E;  on  a  ridge  near  a  primitive  trail.  Detected 
basal  conglomerate  of  the  Pontotoc  Group  lying  adjacent 
to  a  carbonate  unit.  Attitude  of  the  carbonate  unit  is, 
N40°W  strike  with  30°  SE  dip.  Broken  specimens  revealed 
a  brown,  coarse-crystalline  dolomitic  lime  with  some 
silica.  Rocks  weathered  to  a  light-buff  color,  and  con- 
tain numerous  microfractures.  The  well-cemented  conglom- 
erate consists  of  lime  cobbles  and  boulders. 

Station  2.  Location  -  SWj  SW£  SWj  NWj  of  Sec.  6, 
T1S,  R4E;  in  a  ravine.  Attitude  is  N60°W  strike  with  16° 
NE  dip.  Brown,  sandy  dolomite  and  lime  with  rounded  and 
subrounded  fine-grained  quartz  and  chert  particles. 
Shaley  limes  outcrop  just  south  of  this  location  on  a 
east-trending  ridge. 

Station  3.  Location  -  NWj  NWi  SWj  NWj  of  Sec.  6, 
T1S,  R4E.  Resistant  formation,  thickness  of  approximate- 
ly 12  feet,  outcrops  in  a  ravine.  The  bottom  half  of  the 
section  consists  of  a  light  tan  (streaks  of  red)  to  a 
light  brown  (small  amounts  of  silica)  dolomitic  limes. 
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The  top  half  of  the  section  contains  light-brown,  crys- 
talline dolomite;  and  pinkish-colored,  silica-cemented, 
fine  to  medium  grained,  subrounded  to  rounded  quartz 
sand.  Bedding  planes  of  dolomitic  limes  are  easily  dis- 
cerned. The  attitude  of  beds  strike  N60°E  with  a  N25°W 
dip.  Two  swallow  shafts  and  one  small  sink  hole  were 
found  in  and  along  this  ravine. 

The  outcrop  was  mapped  with  the  use  of  aerial  photo- 
graphs, and  then  transferred  to  a  topographic  map.  The 
dip  of  the  beds  in  this  outcrop  area  are  extremely  errat- 
ic, suggesting  complex  geologic  structure. 
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APPENDIX  IV 
ESTIMATE  OF  DISCHARGE  OF  THE  VENDOME  WELL 

Determination  of  the  discharge  from  Vendome  Well 
using  the  Mannings  equation  for  open  channel  flow  can  be 


expressed  as  (Morris,  1972): 
Q  =  VA  = 


2/3   1/2 
1.49  R  /J  SX/Z   A 


(IV. 1) 


where 


n' 

Q  =  discharge  from  the  well  into  an  open 

3 
channel  (ft  /sec), 

V  =  average  velocity  (ft/ sec) , 

n'  =  channel  coefficient, 

S  =  gradient  of  channel  (ft/ft) , 

R  =  A/P  =  hydraulic  radius  (ft) , 

2 
A  =  area  of  channel  (ft  ) , 

P  =  perimeter  of  channel  (ft) . 


Given: 


(1)  V   =  average  surface  velocity  =  3.3  ft/sec 

(2)  n'  =  .02,  channel  coefficient  for  concrete 

(3)  A   =  configuration  of  channel  is  rectangular 

=  bh  (44/12  ft)  (2/12  ft)  =  .611  ft2 

(4)  R  =  bh/(b+h)  =  .61  ft2/(48/12)  ft.  =  .153  ft. 
1/2 


Solving  for  S     yields, 
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1/2       n'V 
S1/A    =  ^-pr      =   .155  (IV. 2) 

1.49  R2/J 

Then  Q  can  be  obtained  by  using  Equation  (IV. 1) 

1.49  A  S1/2  A  , 

Q  =  =   2  ft  /sec. 

n' 
or 

g 

Q  =  472  x  10   gallons/year. 


